please please ...... i need help in this programme is there an error

3 次查看（过去 30 天）
sihem chaib
sihem chaib 2014-2-19
评论： sihem chaib 2014-2-21
M=120; EI=20*10^7; ksi=0;%.04; k=0.2*10^8; v=1; n1=5; n=501;
l=(3*pi/2)*(4*EI/k)^0.25; t1=l/v; %%%%%%%%%%%%%%%% Les conditions initiales %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
v_0=0.001; %m/s;vitesse initiale ;T_0=0.0001; %m déplacement initiale; t1=l/v; ww1=0; Q=16*10^4; %N; chargement
for m=1:n1; for i=1:n; xx(i)=l/(n-1); xx1(i)=xx(i)*(i-1);
phi(i,m)=sin(m*pi*(xx1(i))/l);
for j=1:n;
tt(j)=l/v/(n-1);
tt1(j)=tt(j)*(j-1);
w_d(m)=(((((m^4)*(pi^4)*EI)/M/(l^4)))+k/M)^0.5;
%ksi(m)=c/(2*M*w_d);
w_n(m)=w_d(m)*(1-(ksi)^2)^0.5;
B=T_0;
mu(m)=m*pi*v/l;
A1=2*Q/(M*l);
A2(m)=1/(w_n(m)^4+2*(((ksi*w_d(m))^2-w_n(m)^2)*mu(m)^2)+mu(m)^4); %(((w_n(m))^4)+2*((((ksi*w_d)^2)+((w_n(m)^2)))*(mu(m)^2))+mu(m)^4);
A3(m)=2*ksi*w_d(m)*mu(m);
A4(m)=(mu(m)/w_n(m))*(((ksi*w_d(m))^2)-(w_n(m)^2)+mu(m)^2);
A5(m)=((w_n(m))^2-(mu(m)^2));
T11(j,m)=exp(-ksi*w_d(m)*(tt1(j)));
T1(j,m)=A1*A2(m)*(T11(j,m)*(A3(m)*cos(w_n(m)*tt1(j))+A4(m)*sin(w_n(m)*tt1(j)))-A3(m)*cos(mu(m)*tt1(j))+A5(m)*sin(mu(m)*tt1(j)));
ww1(i,j,m)=phi(i,m)*T1(j,m);
end
end
end
for i=1:n; for j=1:n; ww1(i,j)=0; for m=1:n1; ww1(i,j)=ww1(i,j)+ww1(i,j,m);
%ww1(i,j) = sum(ww1(i,j,:));
end
end
end
plot(ww1(:,3),'b');
hold on

请先登录，再回答此问题。

采纳的回答

Jason R
Jason R 2014-2-19
What is the problem? I was able to run your code when I uncommented the variables k, v, etc. Remember that anything after a % will not be ran. So in line 5, T_0 is not getting initialized because it is commented out.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by