diff options
Diffstat (limited to '2279/CH5/EX5.7/Ex5_7.sce')
| -rw-r--r-- | 2279/CH5/EX5.7/Ex5_7.sce | 50 |
1 files changed, 50 insertions, 0 deletions
diff --git a/2279/CH5/EX5.7/Ex5_7.sce b/2279/CH5/EX5.7/Ex5_7.sce new file mode 100644 index 000000000..b1a9f8517 --- /dev/null +++ b/2279/CH5/EX5.7/Ex5_7.sce @@ -0,0 +1,50 @@ +//Fourier Series coefficients of half-wave rectifier output +//Assume the period of the signal T=1 +t=-0.5:0.01:1; +T = 1; +for i=1:length(t) + if t(i)<T/2 then + x(i)=sin(2*%pi*t(i)); + else + x(i)=0; + end +end +k=-10:10; +for i=1:length(k) + if k(i)==1 then + ak(i)=1/(4*%i); + elseif k(i)==-1 + ak(i)=-1/(4*%i); + else + ak(i)=(cos(k(i)*%pi/2)*exp(-k(i)*%pi/2*-%i))/(%pi-(%pi*k(i)*k(i))); + end +end + + +disp("The fourier series coefficients are...") +disp(ak) +disp("magnitude of Fourier series coefficient") +disp(abs(ak)) +//PLotting frequency spectrum +subplot(2,1,1) +plot(k,abs(ak),'.'); +xtitle("Magnitude Spectrum","k","|ak|"); +for i=1:length(k) + if k(i)==0 | k(i)==3 | k(i)==-3 | k(i)==-5 |k(i)==5 then + phase(i)=0; + elseif k(i)==-1 then + phase(i)=%pi/2; + elseif k(i)==1 then + phase(i)=-%pi/2; + elseif k(i)==-2 | k(i)==-4 + phase(i)=%pi; + elseif k(i)==2 | k(i)==4 + phase(i)=-%pi; + else + phase(i) = 0; + end +end +subplot(2,1,2) +plot(k,phase,'.'); +xtitle("Phase Spectrum","k","angle(ak)"); +disp(phase)
\ No newline at end of file |