Question

need matlab code

Exercise: Use MATLAB to generate a low frequency signal: xy(t) = 3 cos(1400) And a higher frequency signal: xz(t) = 2 cos(10,000st) Now add these two signals together to generate the dual-component signal: y(t) = 3 cos(14007a)+ 2 cos(10,000) Points to be addressed: Again, manually calculate the theoretical magnitude spectra of y(t) and show this in your report (use Microsoft Visio or similar for your diagram). Clearly label all axes! Contrast this working to what you see in MATLAB. Calculate the power of this new signal.

Answer 1

Solution:

MATLAB CODE:

clc
% Defining Parameters for first signal
N1 = 1024;        % Number of Samples
Fs1 = 700;         % Frequency of Signal = 700 Hz  
T1 = 1/Fs1;         % Time Period of Signal
J1 = 2;             % Number of Cycles = 2
C1 = J1*T1;            % Stop time of 2 Cycles
t1 = 0:(1/N1)*C1:C1;    % Time Range for 2 cycles
x1 = 3*cos(2*pi*Fs1*t1); % First Sinusoidal signal

% Plotting First Sinusoidal Signal
subplot(3,1,1)
plot(t1,x1);             % Plot x1 vs t1
title('First Sinusoidal Signal');
xlabel('Time')
ylabel('Magnitude')

% Defining Parameters for Second Signal
N2 = 1024;        % Number of Samples
Fs2 = 5000;         % Frequency of Signal = 5000 Hz  
T2 = 1/Fs2;         % Time Period of Signal
J2 = 2;             % Number of Cycles = 2
C2 = J2*T2;            % Stop time of 2 Cycles
t2 = 0:(1/N2)*C2:C2;    % Time Range for 2 cycles
x2 = 2*cos(2*pi*Fs2*t2); % Second Sinusoidal signal

% Plotting Second Sinusoidal Signal
subplot(3,1,2)
plot(t2,x2);             % Plot x2 vs t2
title('Second Sinusoidal Signal');
xlabel('Time')
ylabel('Magnitude')

% Plotting dual Sinusoidal Signal
y = x1+x2;
subplot(3,1,3)
plot(t1,y);             % Plot y vs t1
title('Dual Sinusoidal Signal');
xlabel('Time')
ylabel('Magnitude')

Image Copy:

8 Editor - C:\Users\Welcome\sig.m Qe.m X sig.mx + 1 clc 2 Defining Parameters for first signal 3 Nl = 1024; Number of Samples 4 Fsl = 700; Frequency of Signal = 700 Hz 5 Ti = 1/Fsl; Time Period of Signal 6 Jl = 2; Number of Cycles = 2 7 Cl = Jl*Tl; Stop time of 2 Cycles tl = 0: (1/Nl) *C1:Cl; Time Range for 2 cycles 9 xl = 3* cos(2pi Fsl*ti); First Sinusoidal signal 10 11 Plotting First Sinusoidal Signal 12 - subplot (3,1,1) 13 plot (ti,xl); $ Plot xl vs tl 14 - title('First Sinusoidal Signal'); 15 xlabel('Time') 16 - ylabel ('Magnitude') 17 18 $ Defining Parameters for Second Signal 19 N2 = 1024; Number of Samples 20 Fs2 = 5000; Frequency of Signal = 5000 Hz 21 - T2 = 1/F32; Time Period of Signal 22 J2 = 2; $ Number of Cycles = 2 23 C2 = J2*T2; Stop time of 2 Cycles 24 t2 = 0: (1/N2) *C2:02; Time Range for 2 cycles 25 x2 = 2*cos (2*pi*F32*t2); $ Second Sinusoidal signal ITTI IIIIIII

27 & Plotting Second Sinusoidal Signal 28 - subplot (3,1,2) 29 - plot(t2,x2); & Plot x2 vs t2 30 - title('Second Sinusoidal Signal'); 31 - xlabel('Time') 32 - ylabel('Magnitude') 33 34 Plotting dual Sinusoidal Signal 35 - y = xl+x2; 36 - subplot (3,1,3) 37 - plot(ti,y); Plot y vs ti 38 - title('Dual Sinusoidal Signal'); 39 - xlabel('Time') 40 - ylabel ('Magnitude') Command Window

Output Waveform:

- х Figure 1 File Edit View Insert Tools Desktop Window Help First Sinusoidal Signal Magnitude 0 0.5 2 2.5 1.5 Time Second Sinusoidal Signal 3 x 10-3 Magnitude 1 0.5 1 1.5 2.5 3 3.5 4 2 Time Dual Sinusoidal Signal *10-4 Magnitude 1 0.5 1 2 2.5 3 1.5 Time *10-3