Code Lab 4
Code Lab 4
[x1,n1]=sigshift(x,n,1);
[y,ny] = sigadd(3.4*x,n,6.7*x1,n1);
end
%%
close all
clear
clc
clear
clc
close all
n = 1:40;
a = 2;b = -3;
% Input Signal
x1 = cos(2*pi*0.1*n);
x2 = cos(2*pi*0.4*n);
% Weighted Output
[y1,ny1]=TT(x1,n);
[y2,ny2]=TT(x2,n);
[y,ny] = sigadd(a*y1,ny1,b*y2,ny2);
% Weighted Input
[x3,n3] = sigadd(a*x1,n,b*x2,n);
[yt,nt] = TT(x3,n3) ;
% Difference Signal
[err,ne] = sigadd(yt,nt,-y,ny);
n = 1:40;
D = 20;
% Input Signal
a = 22;b = -2;
x = a*sin(2*pi*0.1*n) + b*cos(2*pi*0.4*n);
[xd,nd]=sigshift(x,n,D);
[yxd,nxd]=TT(xd,nd);
[y,ny]=TT(x,n);
[yd,nyd]=sigshift(y,ny,D);
[err,ne] = sigadd(yd,nyd,-yxd,nxd);
subplot(3,1,1)
stem(nxd,yxd);
ylabel('mplitude');
title('Output Signal to dealy input y=T(x(n-D))' );grid;
subplot(3,1,2)
stem(nyd,yd );
ylabel('Amplitude');
title(['Dealy Output Signal y(n-D)'])
subplot(3,1,3)
stem(ne,err);
xlabel('Time index n'); ylabel('Amplitude');
title('Difference Signal');grid;