matlab画响应曲线出错,请大侠们帮忙,程序如下:syms f g k m jsyms s ta=g/(j*s^3+m*s^2+k*s);b=ilaplace(a);j=0.003242;m=0.044;k=769;g=k*pi/9-m;b=1/k*(1-exp(-1/2*t/j*m)*cosh(1/2*t/j*(m^2-4*j*k)^(1/2))-exp(-1/2*t/j*m)*m/(m^2-4*j*k)^(1/2)*s
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/20 10:57:20
![matlab画响应曲线出错,请大侠们帮忙,程序如下:syms f g k m jsyms s ta=g/(j*s^3+m*s^2+k*s);b=ilaplace(a);j=0.003242;m=0.044;k=769;g=k*pi/9-m;b=1/k*(1-exp(-1/2*t/j*m)*cosh(1/2*t/j*(m^2-4*j*k)^(1/2))-exp(-1/2*t/j*m)*m/(m^2-4*j*k)^(1/2)*s](/uploads/image/z/10435308-60-8.jpg?t=matlab%E7%94%BB%E5%93%8D%E5%BA%94%E6%9B%B2%E7%BA%BF%E5%87%BA%E9%94%99%2C%E8%AF%B7%E5%A4%A7%E4%BE%A0%E4%BB%AC%E5%B8%AE%E5%BF%99%2C%E7%A8%8B%E5%BA%8F%E5%A6%82%E4%B8%8B%EF%BC%9Asyms+f+g+k+m+jsyms+s+ta%3Dg%2F%28j%2As%5E3%2Bm%2As%5E2%2Bk%2As%29%3Bb%3Dilaplace%28a%29%3Bj%3D0.003242%3Bm%3D0.044%3Bk%3D769%3Bg%3Dk%2Api%2F9-m%3Bb%3D1%2Fk%2A%281-exp%28-1%2F2%2At%2Fj%2Am%29%2Acosh%281%2F2%2At%2Fj%2A%28m%5E2-4%2Aj%2Ak%29%5E%281%2F2%29%29-exp%28-1%2F2%2At%2Fj%2Am%29%2Am%2F%28m%5E2-4%2Aj%2Ak%29%5E%281%2F2%29%2As)
matlab画响应曲线出错,请大侠们帮忙,程序如下:syms f g k m jsyms s ta=g/(j*s^3+m*s^2+k*s);b=ilaplace(a);j=0.003242;m=0.044;k=769;g=k*pi/9-m;b=1/k*(1-exp(-1/2*t/j*m)*cosh(1/2*t/j*(m^2-4*j*k)^(1/2))-exp(-1/2*t/j*m)*m/(m^2-4*j*k)^(1/2)*s
matlab画响应曲线出错,请大侠们帮忙,
程序如下:
syms f g k m j
syms s t
a=g/(j*s^3+m*s^2+k*s);
b=ilaplace(a);
j=0.003242;
m=0.044;
k=769;
g=k*pi/9-m;
b=1/k*(1-exp(-1/2*t/j*m)*cosh(1/2*t/j*(m^2-4*j*k)^(1/2))-exp(-1/2*t/j*m)*m/(m^2-4*j*k)^(1/2)*sinh(1/2*t/j*(m^2-4*j*k)^(1/2)))*g;
t=(0:1:100);
plot(t,b)
报错:
Error using ==> plot
Conversion to double from sym is not possible.
Error in ==> sh2 at 11
plot(t,b)
matlab画响应曲线出错,请大侠们帮忙,程序如下:syms f g k m jsyms s ta=g/(j*s^3+m*s^2+k*s);b=ilaplace(a);j=0.003242;m=0.044;k=769;g=k*pi/9-m;b=1/k*(1-exp(-1/2*t/j*m)*cosh(1/2*t/j*(m^2-4*j*k)^(1/2))-exp(-1/2*t/j*m)*m/(m^2-4*j*k)^(1/2)*s
syms f g k m j
syms s t
a=g/(j*s^3+m*s^2+k*s);
b=ilaplace(a);
j=0.003242;
m=0.044;
k=769;
g=k*pi/9-m;
b=1/k*(1-exp(-1/2*t/j*m)*cosh(1/2*t/j*(m^2-4*j*k)^(1/2))-exp(-1/2*t/j*m)*m/(m^2-4*j*k)^(1/2)*sinh(1/2*t/j*(m^2-4*j*k)^(1/2)))*g;
t=(0:1:100);
b=eval(b)
figure
plot(t,b)
看看你少了哪一句?
还有一个问题就是你的b是一个虚数,所以画出来的图有问题!
祝你学习愉快!