查看: 97|回复: 1|关注: 0

[已答复] ilaplace 不能在xy平面画图 请大神指教!

[复制链接]

新手

5 麦片

财富积分


050


1

主题

1

帖子

0

最佳答案
发表于 3 天前 | 显示全部楼层 |阅读模式
程序:
syms s t
disp('象函数Lf')
Lf = (s+251.3)/(s*s*s*s+2*160*pi*s*s*s+45600*pi*pi*s*s+20000*160*pi*pi*pi*s+10000*25600*pi*pi*pi*pi)
disp('Lf的原函数f')
f = ilaplace(Lf)
subplot(1,2,1)
ezplot(Lf)
subplot(1,2,2)
ezplot(f)

提示:
>> ilaplace
象函数Lf

Lf =

(s + 2513/10)/(s^4 + 320*pi*s^3 + (3865934085194845*s^2)/8589934592 + (6658547215630757*s)/67108864 + 817126680320561/32768)

Lf的原函数f

f =

8589934592*symsum((root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)*exp(root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)*t))/(2*(3865934085194845*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k) + 4123168604160*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)^2*pi + 17179869184*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)^3 + 426147021800368448)), k, 1, 4) + (10793252814848*symsum(exp(root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)*t)/(2*(17179869184*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630\\\r\n757*s5)/67108864 + 817126680320561/32768, s5, k)^3 + 4123168604160*pi*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)^2 + 3865934085194845*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k) + 426147021800368448)), k, 1, 4))/5

错误使用 ezplot (line 162)
不能在 xy 平面中对 8589934592*symsum((root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 +
(6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)*exp(root(s5^4 + 320*s5^3*pi +
(3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5,
k)*t))/(2*(3865934085194845*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 +
(6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k) + 4123168604160*root(s5^4 + 320*s5^3*pi +
(3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)^2*pi +
17179869184*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 +
817126680320561/32768, s5, k)^3 + 426147021800368448)), k, 1, 4) + (10793252814848*symsum(exp(root(s5^4 + 320*s5^3*pi
+ (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5,
k)*t)/(2*(17179869184*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 +
817126680320561/32768, s5, k)^3 + 4123168604160*pi*root(s5^4 + 320*s5^3*pi + (3865934085194845*s5^2)/8589934592 +
(6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k)^2 + 3865934085194845*root(s5^4 + 320*s5^3*pi +
(3865934085194845*s5^2)/8589934592 + (6658547215630757*s5)/67108864 + 817126680320561/32768, s5, k) +
426147021800368448)), k, 1, 4))/5 绘图。

出错 sym/ezplot (line 61)
   h = ezplot(fhandle(f)); %#ok<EZPLT>

出错 ilaplace (line 9)
ezplot(f)

请问大神们应该怎么修改才能画出f的图像呢 谢谢!

新手

13 麦片

财富积分


050


0

主题

17

帖子

2

最佳答案
发表于 昨天 11:37 | 显示全部楼层
H=simplify(vpa(f))
ezplot(H)

替代原来的ezplot(f)
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表