求助 如何在图像中把虚数解标出来

kalfca

b=0;
for j=1:100
    p(j)=0+0.1*j;
  b=p(j);
a=[b-2 b 6-b 4-b];
x=roots(a);
d(j)=x(1)
end
plot(p,d)
hold on

这个程序画的是方程解的图，但是有的解是虚数，matlab好像是直接取实数部分放到图里面去了，请问怎样才能把图里面属于虚数的那些点标成红色

hancy

real(d)代表d的实部，imag(d)是虚部

plot(p,imag(d),'r.')

kalfca

回复 2 # hancy 的帖子

终于画出来了  万分感谢~~

kalfca

回复 2 # hancy 的帖子

才想起来这个图里面的点仅仅只是虚部值的点，不是方程解的点，下面的图就是这个方程的解的图，我希望是将图中那些解的点是虚数用红点标示出来，效果如下，图中的红点是我用图形编辑器画的...我有个思路  就是加个判定语句

hancy

回复 4 # kalfca 的帖子

1. b=0;
   
2. for j=1:100
   
3. p(j)=0+0.1*j;
   
4. b=p(j);
   
5. a=[b-2 b 6-b 4-b];
   
6. x=roots(a);
   
7. d(j)=x(1);
   
8. end
   
9. for i=1:100
   
10. if imag(d(i))~=0
    
11. D(i)=d(i);
    
12. else D(i)=-25;
    
13. end
    
14. end
    
15. plot(p,real(d),p,real(D),'r.')
    
16. axis([0 10 -20 25])
    
17. hold on

复制代码

虽然实现了你的目的，但是过程很繁琐，实际上是根据已知结果对图像的处理
额，我也新手，水平不足。希望有高人提供真正的解决方法！

kalfca

回复 5 # hancy 的帖子

你比我厉害多了 我才是超级新手，为了毕业设计没办法，大学四年没碰过matlab...这个就已经很好很好啦~非常感谢你  我可以向老师交代了~~

kalfca

回复 5 # hancy 的帖子

else D(i)=-25;
end
end
plot(p,real(d),p,real(D),'r.')
axis([0 10 -20 25])
hold on
完了，老大，这有个地方看不懂哇 这个D（i）=-25有什么作用呢 为什么选-25？作用是使点超过图形的表示范围么~~？

hancy

回复 7 # kalfca 的帖子

是啊，下面那个axis（[xmin xmax ymin ymax]）是控制显示范围的
-25在-20下面，不会显示在图上

kalfca

回复 8 # hancy 的帖子

哈哈 这你都想得到

hancy

回复 9 # kalfca 的帖子

这是取巧的办法，往往大型一点的题目就没法用了-_-!!

还是等待高手给个真正的解决办法吧

kalfca

回复 10 # hancy 的帖子

呵呵 好吧 辛苦你了~

meiyongyuandeze

1. clear;
   
2. clc;
   
3. close all
   
4. format long;
   
5. b0=linspace(0,10,100);
   
6. d=arrayfun(@(b)roots([b-2 b 6-b 4-b]),b0,'Uni',0);
   
7. dd=cell2mat(d)

复制代码

用匿名函数+arryfun可以比较简单的完成循环求根，但实在看不出哪里有复根。。。。。所以没画图！
显然你的多项式有一个重根为-1，而且每个b值对应三个根。。。。。

kalfca

回复 12 # meiyongyuandeze 的帖子

这个有复根的啊，具体图形显示不出来我也不知道，我只有带数字进去验算
当b=3的时候 a=(1 3 3 1)
a=[1 3 3 1];
>> roots(a)

ans =

  -1.0000 + 0.0000i
  -1.0000 - 0.0000i
  -1.0000   
这个是结果....

kalfca

回复 12 # meiyongyuandeze 的帖子

只是不知道为什么会出现i  而且i 前面还是0.0000

kalfca

回复 12 # meiyongyuandeze 的帖子

刚才试了下这个程序 有复根 不过很小很小....
dd =

  Columns 1 through 4

  2.000000000000001                       2.053191489361700                       2.112359550561797                       2.178571428571430                     
-0.999999981369628                      -1.000000000000000 + 0.000000016292727i -1.000000000000000 + 0.000000015477330i -0.999999999999999 + 0.000000006252422i
-1.000000018630372                      -1.000000000000000 - 0.000000016292727i -1.000000000000000 - 0.000000015477330i -0.999999999999999 - 0.000000006252422i

这个图怎么画呢 都显示不出来...