不知道y的表达式是怎样的 我就只改了括号,都是小问题,仔细点就行了
1:第一行和第五行最后加个分号,不然t和y就会输出在主界面了
2:第三行最后多个括号,第四行最后少个括号
如下:
t=0:0.1:100;
y=15.62*(0.01+2*exp(2.52*t).*((-5.22)*10^(-3)*cos(9.45*t)-1.39*10^(-3)*sin(9.45*t)))
-10.93*(5.22*10^(-3)+2*exp(7.56*t).*(5.58*10^(-5)*cos(28.35*t)-1.48*10^(-5)*sin(28.35*t)))
-5.89*10^(-3)*exp(1.29*t)+5.56*10^(-3)*exp(13.8*t)
+8.69*((-2.56)*10^(-3)*exp(1.29*t)+2.56*10^(-3)*exp(13.8*t));
plot(t,y,'-')
title('位移功率谱密度曲线')
xlabel('频率 rad/s')
ylabel('锚索中间位移功率谱 (m2*s)')
legend('1阶')
原图:
放大图:
t=0:0.1:100;
y=15.62.*(0.01+2.*exp(2.52.*t).*((-5.22)*10^(-3).*cos(9.45.*t)-1.39*10^(-3).*sin(9.45.*t)))...
-10.93.*(5.22.*10^(-3)+2.*exp(7.56.*t).*(5.58.*10^(-5).*cos(28.35.*t)-1.48.*10^(-5).*sin(28.35.*t))...
-5.89.*10^(-3).*exp(1.29.*t)+5.56.*10^(-3).*exp(13.8.*t))...
+8.69.*((-2.56).*10^(-3).*exp(1.29.*t)+2.56.*10^(-3).*exp(13.8.*t));
plot(t,y,'-')
title('位移功率谱密度曲线')
xlabel('频率 rad/s')
ylabel('锚索中间位移功率谱 (m2*s)')
legend('1阶')
按上面语句试试