dsp1_1.m

利用R的toepliz性质

%  姓名:齐家红 学号：200520203020
clear
clc
a11=-0.195;a21=-1.595;a31=-1.914;a12=0.95;a22=0.95;a32=0.95;
c1=(1-a12)*((1+a12)^2-(a11)^2)/(1+a12)
c2=(1-a22)*((1+a22)^2-(a21)^2)/(1+a22)
c3=(1-a32)*((1+a32)^2-(a31)^2)/(1+a32)
p11=0.5*(-a11+((a11)^2-4*(a12))^0.5)
p12=0.5*(-a11-((a11)^2-4*(a12))^0.5)
p21 =0.5*(-a21+((a21)^2-4*(a22))^0.5)
p22=0.5*(-a21-((a21)^2-4*(a22))^0.5)
p31=0.5*(-a31+((a31)^2-4*(a32))^0.5)
p32=0.5*(-a21-((a21)^2-4*(a22))^0.5)
r11=-a11/(1+a12)
r12=(a11)^2/(1+(a12))-(a12)
r21=-a21/(1+a22)
r22=(a21)^2/(1+(a22))-(a22)
r31=-a31/(1+a32)
r32=(a31)^2/(1+(a32))-(a32)

for m=3:1:255
    r(1)=0.8179;
    r(2)=0.3546;
    r(m)=1.595*r(m-1)-0.95*r(m-2);
end
r;
w=-pi:pi/100:pi;
r1=zeros(1,255);
for m=1:255
    r1(m)=double(r(256-m));
end
r1;
%利用R的toepliz性质
%a是实数，则r也应该是实数
%x是实信号r(-m)=r(m)
pt1=0.0000000;
pt2=0.0000000;
pt=0.00000000;
for k=1:1:255
  pt1=pt1+r(k)*exp(-i*w*k);
end
pt1;
for l=1:1:255
    pt2=pt2+r1(l)*exp(i*w*l);
end
pt2;
pt=pt1+pt2+1;
subplot(2,1,1);
plot(w,pt);
title('DSP作业第一次图（1）  姓名:齐家红 学号：200520203020');
Hw2=1./(1+a21*exp(-i*w)+a22*exp(-i*2*w));
Pw2=c2.*abs(Hw2).^2;
Pw2;
subplot(2,1,2);
plot(w,Pw2);
figure(2);
subplot(3,1,1);
Hw1=1./(1+a11*exp(-i*w)+a12*exp(-i*2*w));
Pw1=c1.*abs(Hw1).^2;
plot(w,Pw1);
title('DSP作业第一次图（2）  姓名:齐家红 学号：200520203020');
subplot(3,1,2);
Hw2=1./(1+a21*exp(-i*w)+a22*exp(-i*2*w));
Pw2=c2.*abs(Hw2).^2;
plot(w,Pw2);
subplot(3,1,3);
Hw3=1./(1+a31*exp(-i*w)+a32*exp(-i*2*w));
Pw3=c3.*abs(Hw3).^2;
plot(w,Pw3);