|
@ -452,7 +452,7 @@ ltiview(unit_transfer_function(1,1), unit_transfer_function(1,2), ... |
|
|
hold off |
|
|
hold off |
|
|
|
|
|
|
|
|
input_signal_frequency = 2000; % Hz |
|
|
input_signal_frequency = 2000; % Hz |
|
|
input_signal_duty_cycle = 20; |
|
|
input_signal_duty_cycle = 40; |
|
|
|
|
|
|
|
|
number_of_periods_sampled = 120000; |
|
|
number_of_periods_sampled = 120000; |
|
|
sampling_time_seconds = number_of_periods_sampled* ... |
|
|
sampling_time_seconds = number_of_periods_sampled* ... |
|
@ -494,30 +494,6 @@ figure(3) |
|
|
semilogx(frequency_vector,Pyy(1:sampling_length_L/2+1)) |
|
|
semilogx(frequency_vector,Pyy(1:sampling_length_L/2+1)) |
|
|
grid on |
|
|
grid on |
|
|
|
|
|
|
|
|
%{ |
|
|
|
|
|
asdf = 1:length(frequency_vector); |
|
|
|
|
|
for i=1:length(frequency_vector) |
|
|
|
|
|
if i<175000 |
|
|
|
|
|
asdf(i) = 1; |
|
|
|
|
|
elseif mod(i,200)==0 |
|
|
|
|
|
[asdf(i),phase] = ... |
|
|
|
|
|
bode(total_transfer_function, ... |
|
|
|
|
|
frequency_vector(i)*2*pi); |
|
|
|
|
|
else |
|
|
|
|
|
asdf(i) = 0; |
|
|
|
|
|
end |
|
|
|
|
|
%{ |
|
|
|
|
|
if mod(i,10000)==0 |
|
|
|
|
|
disp(i) |
|
|
|
|
|
end |
|
|
|
|
|
%} |
|
|
|
|
|
end |
|
|
|
|
|
|
|
|
|
|
|
figure(4) |
|
|
|
|
|
semilogx(frequency_vector,asdf,'.',frequency_vector,P1,'-') |
|
|
|
|
|
grid on |
|
|
|
|
|
%} |
|
|
|
|
|
|
|
|
|
|
|
% Clears unneeded variable from workspace |
|
|
% Clears unneeded variable from workspace |
|
|
clear low_frequency |
|
|
clear low_frequency |
|
|
clear -regexp _numerator$ |
|
|
clear -regexp _numerator$ |
|
|