请大家帮忙看一下

孔雀草

for i=1:5;
  x(i)=0.25*i-0.25;
  A=[0,1,0,0;-16,-0.1,5-20.*x(i),0;0,0,0,1;5,0,-25,-0.1]
B=[eig(A)]
C=B(2:3)
b=abs(imag(C))
plot(x(i),b)
end
为什么就是得不出图？

eight

原帖由 孔雀草 于 2007-5-16 10:23 发表

                               
登录/注册后可看大图

for i=1:5;
  x(i)=0.25*i-0.25;
  A=[0,1,0,0;-16,-0.1,5-20.*x(i),0;0,0,0,1;5,0,-25,-0.1]
B=[eig(A)]
C=B(2:3)
b=abs(imag(C))
plot(x(i),b)
end
为什么就是得不出图？

先生成整个序列，再 plot(x,b)

建议多看看基础书

[ 本帖最后由 eight 于 2007-5-16 11:08 编辑 ]

孔雀草

我不知道该怎么生成整个序列！

孔雀草

大家多多帮忙啊，不胜感激！！

xjzuo

象你这种问问题的态度,我虽知道,却不太想回答.

孔雀草

为什么呀，你就给我指点一下把，我看了半天书，就是看不出来该则么办！！

eight

原帖由 孔雀草 于 2007-5-17 10:59 发表

                               
登录/注册后可看大图

为什么呀，你就给我指点一下把，我看了半天书，就是看不出来该则么办！！

没看懂你要做什么：x(i) 是一个标量，b是一个向量，根本不匹配，画不到图

孔雀草

我想画出x与得到的复特征值虚部绝对值之间的关系图！x每取一个值，得到一个不同的矩阵，特征值也会变化！

eight

原帖由 孔雀草 于 2007-5-17 11:59 发表

                               
登录/注册后可看大图

我想画出x与得到的复特征值虚部绝对值之间的关系图！x每取一个值，得到一个不同的矩阵，特征值也会变化！

以下代码可以运行，不过可能不是你要的结果，自己修改一下，学会举一反三：

1. 
   
2. for i=1:5;
   
3.     x=0.25*i-0.25;
   
4.     A=[0,1,0,0;-16,-0.1,5-20.*x,0;0,0,0,1;5,0,-25,-0.1];
   
5.     B=[eig(A)];
   
6.     C=B(2:3);
   
7.     b=abs(imag(C));
   
8.     x = x*ones(1,length(b));
   
9.     plot(x,b);
   
10.     hold on;
    
11. end
    
12. hold off
    

复制代码

孔雀草

多谢大哥