% 姓名：齊家紅 學號：200520203020
clc
clear
close all
% M=8;
N=200;M=8;w0=1.2*pi;snr=10;p=1;nm=1;cm=0.1;u1=0.05;
c=p/(10^(snr/10));
v=zeros(N,1);w=zeros(200,8);u11=zeros(1,200);e=zeros(100,200);u=zeros(N+7,1);w1=zeros(1,8);u2=zeros(8,200);
for k=1:1:100
   k1=zeros(8,1);
   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;
p=eye(M,M)*cm;
   for n2=1:1:200
       u2(:,n2)=[u(n2+7),u(n2+6),u(n2+5),u(n2+4),u(n2+3),u(n2+2),u(n2+1),u(n2)]';
       k1=1/nm*p*u2(:,n2)/(1+1/nm*u2(:,n2)'*p*u2(:,n2));     
       e(k,n2)=u(n2+7)'.'-u2(:,n2)'*w(n2,:)';
       p=1/nm*p-1/nm*k1*u2(:,n2)'*p;
       w(n2+1,:)=w(n2,:)+(k1*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作業第五次圖(3) 姓名:齊家紅  學號：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);