1谐波幅度对波形合成的影响.docVIP

1.谐波幅度对波形合成的影响 一次、三次、五次、七次和九次谐波合成方波的源代码如下： function[]=sy1() t=0:0.0001:0.1; f=20; s1=sin(2*pi*f*t); s3=sin(2*pi*3*f*t); s5=sin(2*pi*5*f*t); s7=sin(2*pi*7*f*t); s9=sin(2*pi*9*f*t); x=4/pi*(s1+1/3*s3+1/5*s5+1/7*s7+1/9*s9); plot(t,x,r); 2.相位对波形合成的影响 三次谐波相对于基波有270度相移，再将一次、三次、五次、七次和九次谐波合成方波的源代码如下： function[]=sy2() t=0:0.0001:0.1; f=20; s1=sin(2*pi*f*t); s3=cos(2*pi*3*f*t); s5=sin(2*pi*5*f*t); s7=sin(2*pi*7*f*t); s9=sin(2*pi*9*f*t); x=4/pi*(s1+1/3*s3+1/5*s5+1/7*s7+1/9*s9); plot(t,x,k); 3.信号频谱分析 1）、方波信号、正弦波信号的频谱 a.频率为4KHz占空比为50%的方波信号的频谱。 源代码如下： function[]=s3y_1() f=input(enter f:); T=1/f; tao=50; N=4; n=N*T; t=0:1/(50*f):n; x=(sign((tao/100)*T-mod(t,T))+1)/2;% subplot(2,1,1); plot(t,x,r); axis([min(t) max(t) -1 2]); fs=10000;%设定采样频率 y=fft(x); F=(0:length(y)-1)*fs/length(y); subplot(2,1,2); plot(F,abs(y)) 2）、频率为8KHz占空比为50%的方波信号的频谱。 源代码如下： function[]=s3y_1() f=input(enter f:); T=1/f; tao=50; N=4; n=N*T; t=0:1/(50*f):n; x=(sign((tao/100)*T-mod(t,T))+1)/2; subplot(2,1,1); plot(t,x,r); axis([min(t) max(t) -1 2]); fs=10000; y=fft(x); F=(0:length(y)-1)*fs/length(y); subplot(2,1,2); plot(F,abs(y)) （3）正弦波信号的频谱源代码： function[]=zh() f=20;%input(enter f:); t=-0.1:0.0001:0.1; l=200%input(enter l:); k=0:199; x=sin(2*pi*f*t); subplot(211); plot(t,x,r) xk=fft(x,l); subplot(212); plot(k,abs(xk)); 2、观测实时模拟信号（语音）的频谱，源代码如下： function[]=sy3_2() fs=8000; x1=wavread(D:\sying.wav); sound(x1,10000); y1=fft(x1,1024); f=fs*(0:511)/1024; figure(1) plot(x1) title(原始语音信号); xlabel(time n); ylabel(fuzhi n); figure(2) freqz(x1) title(频率响应图) figure(3) subplot(2,1,1); plot(abs(y1(1:512))) title(原始语音信号FFT频谱) subplot(2,1,2); plot(f,abs(y1(1:512))); title(原始语音信号频谱) xlabel(Hz); ylabel(fuzhi); 4、系统的频率域分析 1）、RC电路，1/RC=1000，用MATLAB绘制出系统的幅度谱和相位谱。 源代码如下： function[]=sy4() RC=0.001; w=0:50:5000; H=(1/RC)./(j*w+1/RC); x=abs(H); y=180*angle(H)/pi; %H=