|
程序:
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的图像呢 谢谢!
|
|