数字信号处理仿真作业汇.docVIP

作业1 （1）求 A = [1 -a1 -a2; 0 -(1+a2) 0; 0 -a1 -1]; b = [r0 a1*r0 a2*r0]; x = inv(A)*b; %x = [sigma_v r1 r2]; sigma_v=x(1); r1=x(2); r2=x(3); 只要更改a1和a2的值，得到 a1 a2 P P r(1) r(2) -0.195 0.95 0.0965 0.0975+0.9698i 0.0975-0.9698i 0.1000 -0.9305 -1.595 0.95 0.0323 0.7975+0.5604i 0.7975-0.5604i 0.8179 0.3546 -1.9114 0.95 0.0038 0.9557+0.1914i 0.9557-0.1914i 0.9802 0.9236 产生序列 data_len = 1024; trials = 100; n = 1:data_len; a1 = -1.9114; a2 = 0.95; sigma_v_2 =0.003822; v = sqrt(sigma_v_2) * randn(data_len, 1 , trials); u0 = [0 0]; num = 1; den = [1 a1 a2]; Zi = filtic(num,den,u0); u = filte