R1=1e3; C=10e-6; % R=1K ohms, C=1 uF num = -1; den = [R1*C 1]; T=[0:.001:5]; W=...
R1=1e3; C=10e-6; % R=1K ohms, C=1 uF
num = -1;
den = [R1*C 1];
T=[0:.001:5];
W= logspace(1,3,20)
for i= 1:length(W)
w=W(i);
u=cos(w*T);
y=lsim(H1,u,T);
gain(i)=max(y)/max(u);
plot(T,u,T,y), grid,
title(['x(t) & y(t): w = ' num2str(w)])
shg pause
end semilogx(W,20*log10(gain)),grid,
title('Frequency Response gain H1(w)'); xlabel(‘Freq (r/s)’)
ylabel(‘db(H1)’)
shg
How can this code sweep a sine wave? And measure dB to get a frequency response
Homework Answers
Add Answer to:
R1=1e3; C=10e-6; % R=1K ohms, C=1 uF
num = -1;
den = [R1*C 1];
T=[0:.001:5];
W=...
MATLAB code starts here --------- clear T0=2; w0=2*pi/T0; f0=1/T0; Tmax=4; Nmax=15; %--- i=1; for t=-Tmax: .01:Tmax...
MATLAB
code starts here ---------
clear
T0=2;
w0=2*pi/T0;
f0=1/T0;
Tmax=4;
Nmax=15;
%---
i=1;
for t=-Tmax: .01:Tmax
T(i)=t;
if t>=(T0/2)
while (t>T0/2)
t=t-T0;
end
elseif t<=-(T0/2)
while (t<=-T0/2)
t=t+T0;
end
end
if abs(t)<=(T0/4)
y(i)=1;
else
y(i)=0;
end
i=i+1;
end
plot(T,y),grid, xlabel('Time (sec)'); title('y(t) square wave');
shg
disp('Hit return..');
pause
%---
a0=1/2;
F(1)=0; %dc freq
C(1)=a0;
for n=1:Nmax
a(n)=(2/(n*pi))*sin((n*pi)/2);
b(n)=0;
C(n+1)=sqrt(a(n)^2+b(n)^2);
F(n+1)=n*f0;
end
stem(F,abs,(C)), grid, title(['Line Spectrum: Harmonics = '
num2str(Nmax)]);
xlabel('Freq(Hz)'), ylabel('Cn'), shg
disp('Hit return...');
pause
%---
yest=a0*ones(1,length(T));
for n=1:Nmax
yest=yest+a(n)*cos(2*n*pi*T/T0)+b(n)*sin(2*n*pi*T/T0);...
Can you please help me answer Task 2.b? Please show all work. fs=44100; no_pts=8192; t=([0:no_pts-1]')/fs; y1=sin(2...
Can you please help me answer Task 2.b?
Please show all work.
fs=44100; no_pts=8192;
t=([0:no_pts-1]')/fs;
y1=sin(2*pi*1000*t);
figure;
plot(t,y1);
xlabel('t (second)')
ylabel('y(t)')
axis([0,.004,-1.2,1.2]) % constrain axis so you can actually see
the wave
sound(y1,fs); % play sound using windows driver.
%%
% Check the frequency domain signal. fr is the frequency vector and
f1 is the magnitude of F{y1}.
fr=([0:no_pts-1]')/no_pts*fs; %in Hz
fr=fr(1:no_pts/2); % single-sided spectrum
f1=abs(fft(y1)); % compute fft
f1=f1(1:no_pts/2)/fs;
%%
% F is the continuous time Fourier. (See derivation...
MATLAB
code starts here ---------
clear
T0=2;
w0=2*pi/T0;
f0=1/T0;
Tmax=4;
Nmax=15;
%---
i=1;
for t=-Tmax: .01:Tmax
T(i)=t;
if t>=(T0/2)
while (t>T0/2)
t=t-T0;
end
elseif t<=-(T0/2)
while (t<=-T0/2)
t=t+T0;
end
end
if abs(t)<=(T0/4)
y(i)=1;
else
y(i)=0;
end
i=i+1;
end
plot(T,y),grid, xlabel('Time (sec)'); title('y(t) square wave');
shg
disp('Hit return..');
pause
%---
a0=1/2;
F(1)=0; %dc freq
C(1)=a0;
for n=1:Nmax
a(n)=(2/(n*pi))*sin((n*pi)/2);
b(n)=0;
C(n+1)=sqrt(a(n)^2+b(n)^2);
F(n+1)=n*f0;
end
stem(F,abs,(C)), grid, title(['Line Spectrum: Harmonics = '
num2str(Nmax)]);
xlabel('Freq(Hz)'), ylabel('Cn'), shg
disp('Hit return...');
pause
%---
yest=a0*ones(1,length(T));
for n=1:Nmax
yest=yest+a(n)*cos(2*n*pi*T/T0)+b(n)*sin(2*n*pi*T/T0);...
Can you please help me answer Task 2.b?
Please show all work.
fs=44100; no_pts=8192;
t=([0:no_pts-1]')/fs;
y1=sin(2*pi*1000*t);
figure;
plot(t,y1);
xlabel('t (second)')
ylabel('y(t)')
axis([0,.004,-1.2,1.2]) % constrain axis so you can actually see
the wave
sound(y1,fs); % play sound using windows driver.
%%
% Check the frequency domain signal. fr is the frequency vector and
f1 is the magnitude of F{y1}.
fr=([0:no_pts-1]')/no_pts*fs; %in Hz
fr=fr(1:no_pts/2); % single-sided spectrum
f1=abs(fft(y1)); % compute fft
f1=f1(1:no_pts/2)/fs;
%%
% F is the continuous time Fourier. (See derivation...
Most questions answered within 3 hours.
-
I create an address book where the user enters the name, phone
and email in the...
asked 6 minutes ago -
The production capacity for acrylonitrile
(C3H3N) in the United States exceeds 2
million pounds per year....
asked 14 minutes ago -
explain and comment out your answer
43. How many address lines are required to address a...
asked 20 minutes ago -
A sample of 45 observations is selected from a normal
population. The sample mean is 49,...
asked 35 minutes ago -
A construction company is planning to bid on a building
contract. The bid costs the company...
asked 32 minutes ago -
A firm operating in a purely competitive environment is faced
with a market price of $250....
asked 39 minutes ago -
•Let’s say someone claims the average population size is
600 feet squared and the housing authority...
asked 46 minutes ago -
Cynaide is a deadly poison that blocks the last step in the
electron transport chain of...
asked 50 minutes ago -
Your friend tells you that there is a vending machine on campus
that dispenses M&M packs...
asked 1 hour ago -
What advantages are there to using piperidine rather than
hydroxide as a base?
asked 1 hour ago -
7. The life of a Freeze Breeze electric fan is normally
distributed with a mean 4...
asked 1 hour ago -
1. A 751 mL NaCl solution is diluted to a volume of 1.06 L and a...
asked 1 hour ago