该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
这条程序
function[]=pend() %function增加新的函数
clf% %清除当前图形
figure(1) %生成图形窗口1
clf
%define TF
num=[2.3566 0 0];
den=[1.0000 0.0883 -27.8285 -2.3094 0];
[z1,p1,k1]=tf2zp(num,den) %得到传递函数的零极点
figure(1)
%ask user for controller
numc=input('numc?........');
denc=input('denc?........');
k =input('K?........');
%view compensated system bode
bode(k*conv(numc,num),conv(denc,den))
%view compensated system nyquist
figure(2) %生成图形窗口2
subplot(2,1,1) %生成子图窗
nyquist(k*conv(numc,num),conv(denc,den))
%view compensated CL system impulse response
subplot(2,1,2)
clnum=conv(num,denc); %conv多项式乘法
temp1=k*conv(numc,num);
temp2=conv(denc,den);
clden=polyadd(temp1,temp2);
[z2,p2,k2]=tf2zp(clnum,clden)
impulse(clnum,clden )
打进 numc?........1
denc?........1
K?........1
其它都没什么问题。就是
??? Undefined command/function 'polyadd'.
Error in ==> pend at 25
clden=polyadd(temp1,temp2);
这里编译不了。。请高手指点一�