MATLAB+Integration 源代码

clear clc c=6; % 子载波个数 bits=4000; % 每个信道的比特数 n=c*bits; % 总的传送比特数 data= 2*round(rand(1,n))-1; %产生信源数据 …

function y=LPF(fc); %实现低通的频域形式 global dt df N t f T Bs=N*df/2; %系统带宽 x1=(Bs-fc/2)*N/(2*Bs); x2=(Bs+fc/2)*N/(2*Bs); y=zeros(…

clear all global dt t f df N T %全局变量 close all N=2^20; dt=0.001; %ms df=1/(N*dt); %KHz T=N*dt; %截短时间 Bs=N*df/2; %系统带宽 …

n=1000; x = randint(n,1); %产生信源. h = modem.pskmod (4); %产生调制句柄 y = modulate(h,x); %对信号进行调制 g = modem.pskdemod (h) …

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…