cdf2rdf函数问题

kangarooli

我在编程过程中有这样一个语句

[evects,evalue]=eig(jacoA);
[evects,evalue]=cdf2rdf(evects,evalue);

运行后老是提示错误，但我看不出原因，哪位高手给看看，不胜感激
错误提示如下，我的程序名是qiujie

??? Function 'find' is not defined for values of class 'sym'.
Error in ==> find at 39
  [varargout{1:nargout}] = builtin('find', varargin{:});
Error in ==> cdf2rdf at 19
i = find(imag(diag(d))');
Error in ==> qiujie at 69
[evects,evalue]=cdf2rdf(evects,evalue)

ChaChing

个人水平专业有限, 无法从你的帖子中得到你的问题的完整信息, 不能确定判断出错原因:loveliness:
find里头的d是否为sym?

kangarooli

前两个错误提示不是我程序里的，只有第三个才是我程序里的语句，应该说的是cdf2rdf的源程序吧，我就是这不明白，d应该就是指所求的的特征值

rocwoods

呵呵，MATLAB的错误提示机制是由内到外逐步提示，先提示直接发生错误的地方，然后顺藤摸瓜提示外层哪个函数调用引起的这个错误，然后这样一层层外推，从LZ给的错误信息来看，LZ程序qiujie第69行[evects,evalue]=cdf2rdf(evects,evalue)引发了错误，cdf2rdf第19行调用了find，而find不支持符号类型的查找。
综上分析得出的结论：qiujie里面定义了一些符号变量，evects或evalue是符号型的，请改成数值型。

kangarooli

谢谢，看了你的回答很受益，但是我求出的特征值和特征向量都是数值的啊！如下
evects =

[   -.2003-1.566*i,   -.2003+1.566*i,            .2180,                    .7613]
[    .6427+1.521*i,    .6427-1.521*i,           -2.109,                     .7772]
[ .2589+.8416e-1*i, .2589-.8416e-1*i,        -.2140e-2,           -.2919]
[    1.849+3.373*i,    1.849-3.373*i,           -6.379,                     .4653]
evalue =

[            -2.815,                0.,                                0.,                                      0.]
[                0.,      -.9570e-5+.7971*i,                0.,                                       0.]
[                0.,                     0.,                          -9.185,                                   0.]
[                0.,                     0.,                               0.,          -.9570e-5-.7971*i]
cdf2rdf就是使复数特征值以实数块的形式显示，我发现特征值如下显示的话就没问题，可以运行，下面是我自己输的，这是怎么回事呢
evalue =

[         0+.7971*i,                        0.,                              0.,                     0.]
[                0.,                        0-.7971*i,                        0.,                     0.]
[                0.,                                0.,                       -2.815,                   0.]
[                0.,                                0.,                             0.,             -9.185]

ChaChing

...evects =

[   -.2003-1.566*i,   -.2003+1.566*i,            .2180,                    .7613]
[    .6427+1.521*i,    .6427-1.521*i,           -2.109,                     .7772]
[ .2589+.8416e-1*i, .2589-.8416e-1*i,        -.2140e-2,           -.2919]
[    1.849+3.373*i,    1.849-3.373*i,           -6.379,                     .4653]
...

如此显示就是sym
试比较两者差异
syms x; aa=[1,2;3,4], bb=[1,2;x,3]

kangarooli

你说的例子我能明白，但我的那个里面没有系统量！都是数值啊！
那个e也是系统十进制显示，怎么会是sym呢，请指点

ChaChing

用isnumeric检查下吧!

kangarooli

奇怪，返回值怎么会是假呢，不解，都是数怎么还返回零呢，真诚的感谢各位，很长知识，谢谢

ChaChing

不是请楼主详细比较下, 两者输出型式的差异!:loveliness:
syms x; aa=[1,2;3,4], bb=[1,2;x,4]; bb(2,1)=3
注意下, 那些[]符号