把每个子函数单独存为.m文件就可以了
另外还有个错误，在sdir2cas函数中
如果K*2~=Na，A没有赋值
如果Nb~=0，B没有赋值
输出元素必须赋值
可给A B 赋一个默认值
下面也是巴特沃兹低通滤波器的源码
代码:
wp=30;
ws=50;
rp=1;
rs=30;
[n,Wn]=buttord(wp*2*pi,ws*2*pi,rp,rs,'s');
[b,a]=butter(n,Wn,'s');
[H,Wn]=freqs(b,a);
plot(Wn*rs/(2*pi),abs(H));
grid;
xlabel('频率/Hz');
ylabel('幅值');