%白噪声及有色噪声序列的产生 clear all; close all; L=500; %仿真长度 d=[1 -1.5 0.7 0.1]; c=[1 0.5 0.2]; %D、C多项式的系数(可用roots命令求其根) nd=length(d)-1; nc=length(c)-1; %nd、nc为D、C的阶次 xik=zeros(nc,1); %白噪声初值，相当于ξ(k-1)...ξ(k-nc) ek=zeros(nd,1); %有色噪声初值 xi=randn(L,1); %r

单输入单输出系统噪信比的计算：%SISO系统噪信比的计算 clear all; a=[1 -0.4]; b=[1]; c=[1 -0.3]; %对象参数A、B、C na=length(a)-1; nb=length(b)-1; nc=length(c)-1; %各多项式阶次 n=max(max(na,nb),nc); a0=[a zeros(1,n-na)]; b0=[b zeros(1,n-nb)]; c0=[c zeros(1,n-nc)]; %高次项补0 deltau2=1