求助：工具箱问题

dongchao1207

??? Undefined function or method 'dir2cas' for input
arguments of type 'double'.
程序写好后出现这个问题，是因为我的MATLAB没装相应的工具箱吗？如果是告诉我是哪个工具箱，怎么安装，谢谢:loveliness:

friendchj

在网上找到了dir2cas的源码，感谢原作者

1. function [b0,B,A]=dir2cas(b,a);
   
2. %变直接形式为级联形式
   
3. %[b0,B,A]=dir2cas(b,a)
   
4. %b0=增益系数
   
5. %B=包含各因子系数bk的K行3列矩阵
   
6. %A=包含各因子系数ak的K行3列矩阵
   
7. %a=直接型分子多项式系数
   
8. %b=直接型分母多项式系数
   
9. 
   
10. %计算增益系数
    
11. b0=b(1);b=b/b0;
    
12. a0=a(1);a=a/a0;
    
13. b0=b0/a0;
    
14. %将分子、分母多项式系数的长度补齐进行计算
    
15. M=length(b);N=length(a);
    
16. if N>M
    
17. b=[b zeros(1,N-M)];
    
18. elseif M>N
    
19. a=[a zeros(1,M-N)];N=M;
    
20. else
    
21. NM=0;
    
22. end
    
23. %级联型系数矩阵初始化
    
24. K=floor(N/2);B=zeros(K,3);A=zeros(K,3);
    
25. if K*2==N
    
26. b=[b 0];
    
27. a=[a 0];
    
28. end
    
29. %根据多项式系数利用函数roots求出所有的根
    
30. %利用函数cplxpair进行按实部从小到大的成对排序
    
31. broots=cplxpair(roots(b));
    
32. aroots=cplxpair(roots(a));
    
33. %取出复共轭对的根变换成多项式系数即为所求
    
34. for i=1:2:2*K
    
35. Brow=broots(i:1:i+1,:);
    
36. Brow=real(poly(Brow));
    
37. B(fix(i+1)/2,:)=Brow;
    
38. Arow=aroots(i:1:i+1,:);
    
39. Arow=real(poly(Arow));
    
40. A(fix(i+1)/2,:)=Arow;
    
41. end

复制代码

dongchao1207

谢谢主任:@D 我刚接触MATLAB源码怎么用还是不懂:@L 能再指点下吧:loveliness:

1011820111

放在matlab的work目录下就可以了，.M文件都可以的，这样就可以直接调用了

dongchao1207

:@L 麻烦问下怎么放。。。是编辑成。m文件放哪里就可以了吗？

dongchao1207

我会了，谢谢两位了:kiss:

ChaChing

看看参考下
http://forum.vibunion.com/forum/thread-75501-1-1.html