matlab中的程序过不去，我是初学者，请各位大虾帮帮忙！

显示全部楼层 · 2010-12-17 20:19:40

clear all
syms b w w1 A t h r;
w=0.3*pi;A=0.2;h=0.4;
t=linspace(0,50,25000);
xc=dsolve('D2x=-2*b*Dx-w^2*x+h*cos(w1*t)','Dx(0)=0','x(0)=0','t');
b=input('b=');
i=1
for k=0.1:0.01;2.5

r(i)=k;w1=k*w;

xc1=eval(xc);

A(i)=max(real(xc1));

i=i+1;
end
plot(r,A), hold on
%b=0.1或者0.2

千问 · 2010-12-17 20:19:40

clear allw=0.3*pi;A=0.2;h=0.4;t=linspace(0,50,25000);xc=dsolve('D2x=-2*b*Dx-w^2*x+h*cos(w1*t)','Dx(0)=0','x(0)=0','t');b=input('b=');k=0.1:0.01:2.5;r=zeros(1,length(k));w1=zeros(1,length(k));for i=1:length(k)
r(i)=k(i);
w1(i)=k(i)*w;
xc1=eval(xc);
A(i)=max(real(xc1));endplot(r,A), hold

千问 · 2010-12-17 20:19:40

你说的是.m文件吧，直接在m文件的窗口中按F5