dsp6.m

LMS的一个好算法

% 姓名：齐家红 学号：200520203020
clc
clear
close all
% M=8;
N=200;M=8;w0=1.2*pi;snr=10;p=1;u1=0.05;
c=p/(10^(snr/10));
v=zeros(N,1);w=zeros(200,8);u11=zeros(1,200);u2=zeros(200,8);e=zeros(100,200);u=zeros(N+7,1);w1=zeros(1,8);
for k=1:1:100
    v=randn(N,1)*c^0.5;
for n=8:1:N+7
    u(n)=exp(j*pi*(n-7)-j*pi)+exp(j*w0*(n-7)-0.79*pi)+v(n-7);
end
u;
   for n2=1:1:200
       u2(n2,:)=[u(n2),u(n2+1),u(n2+2),u(n2+3),u(n2+4),u(n2+5),u(n2+6),u(n2+7)];
       u11(n2)=w(n2,:)*u2(n2,:)';
       e(k,n2)=u(n2+7)-u11(n2);
       w(n2+1,:)=w(n2,:)+u1*u2(n2,:)*e(k,n2);
   end
     w1=w1+w(n2+1,:)/100;
end
   w1;
   Jn2=0;
for k1=1:1:100
    Jn2=Jn2+(abs(e(k1,199)))^2/100;
end
Jn2;
L1=0 ;
 for k3=1:1:200;
     for k2=1:1:100
         L1=L1+abs(e(k2,k3))^2/100;
     end
   L(k3)=L1;
   L1=0;
end
k=1:1:200;
subplot(2,1,1);
plot(L(k));
title('DSP作业第五次图(2) 姓名:齐家红 学号：200520203020');
a=-w1'.';
w=0:0.01:2*pi;
t=zeros(1,length(w));
for l=1:M
    t=t+a(l)*exp(-i*l*w);
end
Par8=c./abs(t+1).^2;
subplot(2,1,2);
plot(w,Par8);