huajian.m

详细介绍了Mathlab的使用方法与技巧

function [systf,syszpk]=huajian(sys)
%传递函数化简函数由3部分组成。[systf,syszpk]=huajian(sys)：主函数，完成化简任务；
%[l,li,r,ri,f,fi,x,xi]=genpai(gen)：将传递函数的分子或分母的根按照零根矩阵、实数根矩阵、复数根矩阵、虚根矩阵予以分类；
%[z,p]=duixiao(sumz,nzn,sump,npn):将传递函数的分子和分母的根按照零根、实数根、复数根、虚根分别调用将相等的根予以对消。
%在huajian函数中组成新的系统传递函数
[num,den]=tfdata(sys,'v');
[z,p,k]=tf2zp(num,den);
nz=length(z);np=length(p);iz=0;ip=0;sumz=0;sump=0;
[zl,zli,zr,zri,zf,zfi,zx,zxi]=genpai(z);
[pl,pli,pr,pri,pf,pfi,px,pxi]=genpai(p);
[z1,p1]=duixiao(zl,zli,pl,pli);
%b=[zr,zri,pr,pri];
[z2,p2]=duixiao(zr,zri,pr,pri);
[z3,p3]=duixiao(zf,zfi,pf,pfi);
[z4,p4]=duixiao(zx,zxi,px,pxi);
z=[z1,z2,z3,z4]';p=[p1,p2,p3,p4];
[num,den]=zp2tf(z,p,k);
systf=tf(num,den);
syszpk=zpk(z,p,k);

     



%[num,den]=zp2tf(z,p,k)
%sys=tf(num,den)