Question

1.  close all;

2.  clear all;clc;

3.   

4.   t=1:0.025:5;

5.   desired=5*sin(2*3.*t);

6.   

7.   noise=5*sin(2*50*3.*t);

8.   

9.   refer=5*sin(2*50*3.*t+ 3/20);

10.  

11.  primary=desired+noise;

12.  

13.  subplot(4,1,1);

14.  plot(t,desired);

15.  ylabel('desired');

16.  

17.  

18. subplot(4,1,2);

19.  plot(t,refer);

20.  ylabel('refer');

21.  

22.  

23.  

24.  subplot(4,1,3);

25.  plot(t,primary);

26.  ylabel('primary');

27.  

28.   

29.  

30.  order= 2;

31.  mu=0.005;

32.  n=length(primary)

33.  delayed=zeros(1,order);

34.  adap=zeros(1,order);

35.  cancelled=zeros(1,n);

36.  

37.  for k=1:n,

38.      delayed(1)=refer(k);

39.      y=delayed*adap';

40.      cancelled(k)=primary(k)-y;

41.      adap = adap + 2*mu*cancelled(k) .* delayed;

42.      delayed(2:order)=delayed(1:order-1);

43.  end

44.  

45.  subplot(4,1,4);

46.  plot(t,cancelled);

47.  ylabel('cancelled');