通信原理課程設(shè)計-文庫吧
2024-12-29 19:19 本頁面
【正文】 ,50)。 %注意：MATLAB 中數(shù)組下標(biāo)從 1 開始x(1)=1。subplot(3,1,1)。stem(x)。title(39。單位沖擊信號序列39。)。 b單位階躍序列程序如下：ns=0。nf=10。n0=2。n=[ns:nf]。x=[(nn0)=0]c正弦信號及其頻譜t=1::1。x=5*sin(2*pi*10*t)。N=length(x)。fx=fft(x)。df=100/N。n=0:N/2。f=n*df。subplot(211)。plot(t,x)。grid。 subplot(212)。plot(f,abs(fx(n+1))*2/N)。 grid。d周期鋸齒波sawtooth()Fs=10000。t=0:1/Fs:1。x1=sawtooth(2*pi*50*t,0)。x2=sawtooth(2*pi*50*t,1)。subplot(2,1,1),plot(t,x1),axis([0,1,1])。title(39。鋸齒波139。)。subplot(2,1,2),plot(t,x2),axis([0,1,1])。title(39。鋸齒波239。)f實(shí)指數(shù)序列y(n)=2na=2。x=a.^n。 stem(n,x)。title(39。實(shí)指數(shù)序列39。)e周期方波square()Fs=10000。t=0:1/Fs:1。x1=square(2*pi*50*t,20)。subplot(211),plot(t,x1),title(39。周期方波39。)。axis([0,])。g sin2πf1t*cos2πf2t f1=50Hz f2=2000Hzf1=50。f2=2000。t=0::10。a=sin(2*pi*f1*t)。b=cos(2*pi*f2*t)。y=a.*b。 plot(t,y)。title(39。sin2pif1t*cos2pif2t39。)（2）產(chǎn)生一條2到2之間的Sa(200t)曲線。t=linspace(2,2)。 y=sinc(200*t)。 plot(t,y)。title(39。y=sa(200t)39。)（3）產(chǎn)生下面信號，并繪出頻譜 t 0tt0/4s(t)= t+ t0/4 t0/4t 3t0/4 tt0 3t0/4t t0假設(shè)t0=源程序：clear all。close all。clc。t0=。t=0::。s=t.*(t0amp。t=t0/4)+(t+t0/4).*(t=t0/4amp。t=3*t0/4)+(tt0).*(t=3*t0/4amp。tt0)。plot(t,s) 1．模擬調(diào)制與解調(diào)的仿真（1） 信號的DSB調(diào)制與解調(diào)clear all。close all。clc。t0=。dt=。t=[0:dt:1]。fc=phi0=0。ct=cos(2*pi*fc*t+phi0)。 %載波信號fm=10。fun = @(t)1.*(t=0amp。t=t0/3)+(2).*(t=t0/3amp。t=2*t0/3)+0.*(t=2*t0/3)。mt=fun(t)。ut=*mt.*ct。figure(1)subplot(3,1,1),plot(t,mt),title(39。mt調(diào)制信號39。)subplot(3,1,2),plot(t,ct),title(39。ct載波信號39。)subplot(3,1,3),plot(t,ut),title(39。ut已調(diào)信號39。) [m,n]=size(ut)。ni=*randn(m,n)。u0=ut+ni。figure(2)subplot(3,1,1),plot(t,ni),title(39。ni高斯白噪聲39。)subplot(3,1,2),plot(t,u0),title(39。u0=ut+ni已調(diào)信號+高斯白噪聲39。)w1=2*dt*(fc2*fm)。w2=2*dt*(fc+2*fm)。[b,a]=butter(4,[w1,w2],39。bandpass39。)。u1=filter(b,a,u0)。subplot(3,1,3),plot(t,u1),title(39。u1信號進(jìn)入帶通濾波器39。)figure(3)u2=u1.*ct。subplot(3,1,1),plot(t,u2),title(39。u2與載波相乘后的波形39。)B=2*fm。wn3=2*dt*B。[b,a]=butter(4,wn3,39。low39。)。u3=filter(b,a,u2)。subplot(3,1,2),plot(t,u3),title(39。u3解調(diào)后的信號39。)subplot(3,1,3),plot(t,u3),hold on,title(39。解調(diào)后的信號與原調(diào)制信號比較39。)plot(t,mt,39。r39。)。dt=t(2)t(1)。 % 采樣周期f=1/dt。 % 采樣頻率(Hz)X=fft(ut)。 % 計算x的快速傅立葉變換X
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1