Matlab+Simulation 源代码

function X=T2F(x) global dt df N t f T %x为时域的取样值矢量，X为x的傅氏变换，X与x长度相同,并为2的整幂。　 H=fft(x); X=[H(N/2+1:N),H(1:N/2)]*dt; end

M = 16; %M进制 x = randint(5000,1,M); y=modulate(modem.qammod(M),x); %调制 ynoisy = awgn(y,15,'measured'); %加噪声 scatterplot(y)…

function x=F2T(X) global dt df t f T N %x为时域的取样值矢量，X为x的傅氏变换，X与x长度相同并为2的整幂 X=[X(N/2+1:N),X(1:N/2)]; x=ifft(X)/dt; end

numbits = 99; %每次迭代的比特数 numchans = 2; %信号的信道数 nsamp = 16; %每个符号的采样率 num=10; numerrs = 0; % Number of bit errors seen so…

%impinvar_bilinear.m %chebyshev1模拟原型滤波器 Omegap=0.2*pi;Rp=1; %设计指标 Omegas=0.3*pi;As=16; [n1,Wn1]=cheb1ord(Omegap,Omegas,…

x = rand(100,1); %输入数据 num = rcosine(1,8,'sqrt'); %滤波器的转移函数 y1 = rcosflt(x,1,8,'filter',num); %在发端对数据进行滤波 z1 = rcosflt(y1,1,8,…

Fd= 1; %输入信号的采样率 Fs = 8; %输出信号的采样率 Delay = 3; %滤波器的群时延 R = 0.5; %滚降因子 [yf,tf] = rcosine(Fd,Fs,'fir',R,Delay); plot(yf) g…

M=512; %滤波器长度 L=M; %每个数据块的长度 N=2*M; %FFT的长度 %% 需要添加u, y,d,e等数据，该程序就能给出类似于11.4节LMS算法的仿真图 w=zeros(1,M); wf=fft([w,z…

len=20000; %假设用前2000个数据作为均衡器的训练数据 Tlen=2000; step=0.001; N=64; %均衡器抽头个数 s=zeros(1,len); s1=zeros(1,len); x=zeros(1,N); w…

u = 0.00005; %迭代步长 n = 64; %WHT-LMS的抽头数必须是2的幂 h = zeros(n ,1); %抽头系数，n为抽头个数 x = zeros(1 ,n); %滤波器抽头输入 %% 注意：读者需要给出输入…