正弦波的源程序.docVIP

正弦波的源程序：(一)，用到的函数 1,f2t函数function x=f2t(X)global dt df t f T N%x=f2t(X)%x为时域的取样值矢量%X为x的傅氏变换%X与x长度相同并为2的整幂%本函数需要一个全局变量dt(时域取样间隔)X=[X(N/2+1:N),X(1:N/2)];x=ifft(X)/dt;end 2,t2f函数。 function X=t2f(x)global dt df N t f T%X=t2f(x)%x为时域的取样值矢量%X为x的傅氏变换%X与x长度相同,并为2的整幂。　%本函数需要一个全局变量dt(时域取样间隔)H=fft(x);X=[H(N/2+1:N),H(1:N/2)]*dt;end (二),主程序。 1，%（1）绘出正弦信号波形及频谱global dt df t f N　close allk=input(取样点数＝2^k, k取10左右);if isempty(k), k=10; endf0=input(f0=取1(kz)左右);if isempty(f0), f0=1; endN=2^k;dt=0.01; %msdf=1/(N*dt); %KHzT=N*dt; %截短时间Bs=N*df/2; %系统带宽f=[-Bs+df/2:df:Bs]; %频域横坐标t=[-T/2+dt/2:dt:T/2]; %时域横坐标s=sin(2*pi*f0*t); %输入的正弦信号S=t2f(s); %S是s的傅氏变换a=f2t(S); %a是S的傅氏反变换a=real(a);as=abs(S);subplot(2,1,1) %输出的频谱 plot(f,as,b);gridaxis([-2*f0,+2*f0,min(as),max(as)]) xlabel(f (KHz))ylabel(|S(f)| (V/KHz)) %figure(2)subplot(2,1,2)plot(t,a,black) %输出信号波形画图gridaxis([-2/f0,+2/f0,-1.5,1.5])xlabel(t(ms))ylabel(a(t)(V))gtext(频谱图) 最佳基带系统的源程序：(一)，用到的函数 f2t函数和t2f函数。代码 (二),主程序 globaldt t f df N Tclose allclear Eb_N0 Pek=input(取样点数＝2^k, k取13左右);if isempty(k), k=13; endz=input(每个信号取样点数＝2^z, zk);if isempty(z), z=3; endaa=input(滚降系数=[0.5]);if aa==[],aa=0.5;endN=2^k;L=2^z;M=N/L;Rb=2; %码速率是2Mb/sTs=1/Rb; %码元间隔dt=Ts/L; %时域采样间隔df=1/(N*dt); %频域采样间隔T=N*dt; %截短时间Bs=N*df/2; %系统带宽f=[-Bs+df/2:df:Bs]; %频域横坐标t=[-T/2+dt/2:dt:T/2]; %时域横坐标g=sin(pi*t/Ts).*cos(pi*t*aa/Ts)./[pi*t/Ts.*(1-4*t.^2*aa^2/Ts^2)]; %升余弦脉冲波形GG=t2f(g);GG=abs(GG); %升余弦脉冲的傅式变换GT=sqrt(GG);GR=GT; %最佳系统的发送接收滤波器的傅式变换for l1=1:20;Eb_N0(l1)=(l1-1); %Eb/N0 in dBeb_n0(l1)=10^(Eb_N0(l1)/10);Eb=1;n0=Eb/eb_n0(l1); %信道的噪声谱密度 sita=n0*Bs; %信道中噪声功率n_err=0; %误码计数 for l2=1:20;b=sign(randn(1,M));s=zeros(1,N); %产生冲激序列s(L/2:L:N)=b/dt;SS=t2f(s);S=SS.*GG; %信道的傅式变化a=f2t(S);a=real(a); %不加噪声的输出 n_ch=sqrt(sita)*randn(size(t)); %信道噪声N_CH=t2f(n_ch);nr=real(f2t(N_CH.*GR)); %输出噪声sr=a+nr; %接收信号y=sr(L/2-L/4:L:N-L/4); %取样bb=sign(y); %判决n_err=n_err+length(find(bb~=b));%错误累计endPe(l1)=n_err/(