matlab結(jié)業(yè)報告

1、MATLAB課程設計報告 學 院： 電氣信息工程學院 班 級： 通信07-2 姓 名： 班級學號： 指導教師： 王鋼 題目一：利用壓控振蕩器模塊產(chǎn)生一個受10Hz正弦波控制，中心頻率為100HZ，頻偏范圍為50Hz到150Hz的振蕩信號，并用示波器模塊和頻譜儀模塊觀察輸出信號的波形和頻譜。壓控振蕩器模塊系統(tǒng)仿真的參考模型 如圖1-1所示 圖1-1壓控振蕩器模塊系統(tǒng)仿真波形和頻譜圖 如圖1-2、圖1-3所示 圖1-2 圖1-3題目三：雙邊帶抑制載波的調(diào)制系統(tǒng)進行仿真。設傳輸?shù)幕鶐盘枮檎也?，幅度?，頻率范圍為110Hz，載波頻率為100Hz.2、題目三：雙邊帶抑制載波的調(diào)制系統(tǒng)仿真模型 如圖

2、3-1所示 圖3-1雙邊帶抑制載波的調(diào)制系統(tǒng)仿真波形和頻譜圖 如圖3-2、圖3-3所示 圖3-2 圖3-3題目五：設計一個模擬低通濾波器，fp=3400Hz,fs=6000Hz,Rp=3dB，Rs=30dB，分別用巴特沃斯和橢圓濾波器原型，求出其3dB截止頻率和濾波器階數(shù)，傳遞函數(shù)并做出頻譜、相頻特性曲線。1、巴特沃斯低通濾波器Wp=2*pi*3400; %通帶截止角頻率Ws=2*pi*6000; %阻帶截止角頻率Rp=3; %通帶最大衰減Rs=30; %阻帶最小衰減N,Wc = buttord(Wp,Ws,Rp,Rs,'s'); %求巴特沃斯低通濾波器的階數(shù)和3db截止角頻率

3、B,A = butter(N,Wc,'s'); %求傳遞函數(shù)subplot(2,2,1);k=0:511;fk=0:8000/512:8000;wk=2*pi*fk;Hk =freqs(B,A,wk); %頻率響應plot(fk/1000,20*log10(abs(Hk); %頻率響應的幅度的曲線圖 title('巴特沃斯低通濾波器的幅頻特性'); xlabel('頻率/kHz'); ylabel('幅度/dB') axis(0,8,-40,5)巴特沃斯低通濾波器幅頻特性 如圖5-1所示 圖5-12、橢圓低通濾波器Wp=2*pi*

4、3400; %通帶截止角頻率Ws=2*pi*6000; %阻帶截止角頻率Rp=3; %通帶最大衰減Rs=30; %阻帶最小衰減N,Wpo = ellipord(Wp,Ws,Rp,Rs,'s'); %求橢圓低通濾波器的階數(shù)和3db截止角頻率B,A = ellip(N,Rp,Rs,Wpo,'s'); %求傳遞函數(shù)subplot(2,2,1);fk=0:6000/512:6000;wk=2*pi*fk;Hk =freqs(B,A,wk); %頻率響應plot(fk/1000,20*log10(abs(Hk); %頻率響應的幅度的曲線圖 title('橢圓低通濾

5、波器的幅頻特性'); xlabel('頻率/kHz'); ylabel('幅度/dB') axis(0,12,-70,5)橢圓低通濾波器幅頻特性 如圖5-2所示 圖5-2題目六：以數(shù)字調(diào)制的通帶模型建立ASK、2PAM、BPSK以及FSK的調(diào)制模型并觀察調(diào)制輸出波形和頻譜。仿真系統(tǒng)參數(shù)自己設定。4.題目六：ASK模型：clearn=1:8192;m=1:128;x(n)=randint(1,8192,2);x=x(n)'y(n)=zeros(1,8192);z(m)=zeros(1,128);for n=1:8192 for m=1:128 if

6、 n=64*m-63 z(m)=x(n); if m=ceil(n/64) y(64*m-63):(64*m)')=z(m); end end endendn=1:8192;rm2=y(n);x2=rm2; %產(chǎn)生基帶信號n=1:(213);x1=cos(n.*1e9*2*pi/4e9); %載頻1GHzx=x1.*x2; % ASK 頻帶調(diào)制b=blackman(213); % 窗函數(shù)X= b'.*x; % ASK 頻帶調(diào)制加窗x3=ones(1,64) zeros(1,8128); % 獨個基帶信號碼元y1=X(1:(213);y4=x1.*x3;Y1=fft(y1,(21

7、3);magY1=abs(Y1(1:1:(212)+1)/(200); % 求FFTY4=fft(y4,(213);magY4=abs(Y4(1:1:(212)+1)/(37); % 求FFTk1=0:(212);w1=(2*pi/(213)*k1;u=(2*w1/pi)*1e9;figure(1)subplot(2,1,1);plot(u,magY1,'b',u,magY4,'r');gridtitle('ASKr');axis(4e8,1.6e9,0,1.1)X2=b'.*x2; %基帶信號加窗y2=X2(1:(213);Y2=fft

8、(y2,(213);magY2=abs(Y2(1:1:(212)+1)/(200)+eps; % 求FFTk1=0:(212);w1=(2*pi/(213)*k1;u=(2*w1/pi)*1e9;Y3=fft(x3,(213);magY3=abs(Y3(1:1:(212)+1)/(35)+eps;subplot(2,1,2);semilogy(u,magY2,'b',u,magY3,'r');gridtitle('ASKr-modulation');axis(0,1.2e9,3e-2,3)figure(2)subplot(2,1,1);plot(n,x2);title('ASKr');axis(0,640,-0.2,1.2); gridsubplot(2,1,2);plot(n,x);axis(0,640,-1.2,1.2);