【求助】请大家帮我看看错在哪了，万分感激（非常简单的程序）！

tangxi

大家好，我是Matlab初学者，编的第一个程序，老是出错，郁闷的不行了，请各位帮我看看哪出问题了，谢谢了。

syms w;
ZF=j*(300*w-2700/w); M=300;Rm=180;S=27000;
Zm=(Rm-j*S/w)/(1+(Rm-j*S/w)/ZF)+j*w*M;
for w=0:1:500
    X1=1/(j*w*Zm);
    figure;
    plot(w,abs(X1)); %想画w和X1振幅的关系图，x轴为频率，y轴为X1振幅
end
title('Displacement amplitude X1 of mass M1');
xlabel('w(rad/s)');
ylabel('amplitude(m)');

[ 本帖最后由 tangxi 于 2008-10-6 14:28 编辑 ]

花如月

根据错误提示，一步一步改

sigma665

不知道你说的是什么错误

tangxi

谢谢指点，我也是这么操作的，但后面的错误我真不知道怎么改了

tangxi

比如现在系统提示我“plot(w,abs(X1)); ”错误
可我真不知道哪错了

ChaChing

plot(w,abs(X1))中的X1是symbol object
不能如此画吧

sigma665

第一句去掉
for w=0:1:500
ZF=j*(300*w-2700/w); M=300;Rm=180;S=27000;
Zm=(Rm-j*S/w)/(1+(Rm-j*S/w)/ZF)+j*w*M;

    X1=1/(j*w*Zm);
    figure;
    plot(w,abs(X1)); hold on;%想画w和X1振幅的关系图，x轴为频率，y轴为X1振幅
end

sogooda

小西，你害我matlab不停的画图，最后只好全关了。呵呵。。。
改成下面这样试试，没有错误了，不过好像还是没有绘出图来，楼主检查一下公式吧。

1. 
   
2. M=300;Rm=180;S=27000;
   
3. for w=1:500
   
4.     ZF=j*(300*w-2700/w);
   
5.     Zm=(Rm-j*S/w)/(1+(Rm-j*S/w)/ZF)+j*w*M;
   
6.     X1(w)=1/(j*w*Zm);
   
7. end
   
8.     figure;
   
9.     plot(w,abs(X1)); hold on;%想画w和X1振幅的关系图，x轴为频率，y轴为X1振幅
   

复制代码

[ 本帖最后由 sogooda 于 2008-10-6 21:09 编辑 ]

tangxi

您好，能不能说得再详细一些，我是菜鸟，没能看明白，谢谢了

tangxi

小西的程序让我死机了:@L

tangxi

唉，公式是我从书本上搬过来的，我现在已经晕了:@L

sigma665

他原来就是这么画的:@L

ChaChing

不好意思, 人老了还在使用旧版matlab, 习惯 & 懒了!
symbolic object我使用较少, 还没空消化完手册, 其他人或许可帮忙

没头没脑, 好像被考试的感觉 ! 不确定你到底是要什麽?
仅能稍修一下, 是你要的吗?
w=1:500; ZF=j*(300*w-2700./w); M=300; Rm=180; S=27000;
Zm=(Rm-j*S./w)/(1+(Rm-j*S./w)/ZF)+j*w*M; X1=1./(j*w.*Zm);
plot(w,abs(X1)); title('Displacement amplitude X1 of mass M1');
xlabel('w(rad/s)'); ylabel('amplitude(m)');
还有自己必须时常找Help文件, 个人也是

tangxi

牛人啊，这么简单就搞定了，这就是我要的，谢谢！

ChaChing

我还以为只有我死机! 这件事突显了论坛的一些现象
  1.很多人都复制指令就直接至matlab执行, 未事先看过指令
  2.许多人回覆, 似乎没先搞清楚/执行过问题指令, 甚至回覆的指令自个儿未执行过
这或许只是个人感觉而已, 不针对特定人, 因我也是犯同样的错
说说大家共勉之