一个很简单的程序，但是编写的很经典，都是用于做论文的

这是我的通信原理作业，不想手画（其实是手残画太丑），所以写了个程序。手画俩小时不如写程序二十分钟。。。可议用于通信原理三种键控的学习

内容详细介绍了基于FPGA的数字模拟调制技术，内有相关代码，其中对相对复杂的MSK调制原理，GMSK的高斯低通滤波器设计做了相关说明，介绍了在VIVADO平台的DDS正弦信号发生器、根据FIR滤波器设计高斯低通等经典IP的使用。适合新手和想做相关知识的大学生下载学习，希望能对你有所帮助。

在MATLAB环境下编程实现的数字调制信号样式自动识别，可以识别出2ASK、4ASK、2PSK、4PSK、2FSK、4FSK和16ＱＡＭ。程序包括４个部分，调制信号源的仿真，信道的模拟、瞬时参数的提取、样式的识别。注意：样式识别中的门限值需要自己统计确定

2011年，华东交通大学，课程设计，模拟仿真，自己做的，学弟学妹们可以参考。基于simulink，按照书上的内容自己做的。文档都写好了

通信原理课程设计—实现ask、fsk、psk、dpsk调制解调，包括matlab源程序，可以运行。

Physical-layer security analysis of PSK quantum-noise randomized cipher in optically amplified links

用于生成4psk调制的小程序，很好用，试试看

通信原理课程设计报告（ASK FSK PSK Matlab仿真--数字调制技术的仿真实现及性能研究）报告一切搞定，只写名字了 程序下载地址：http://download.csdn.net/source/1921154

clear all; close all; fs=8e5;%抽样频率 fm=20e3;%基带频率 n=2*(6*fs/fm); final=(1/fs)*(n-1); fc=2e5; % 载波频率 t=0:1/fs:(final); Fn=fs/2;%耐奎斯特频率 %用正弦波产生方波 %========================================== twopi_fc_t=2*pi*fm*t; A=1; phi=0; x = A * cos(twopi_fc_t + phi); % 方波 am=1; x(x>0)=am; x(x<0)=-1; figure(1) subplot(321); plot(t,x); axis([0 2e-4 -2 2]); title('基带信号'); grid on car=sin(2*pi*fc*t);%载波 ask=x.*car;%载波调制 subplot(322); plot(t,ask); axis([0 200e-6 -2 2]); title('PSK信号'); grid on; %===================================================== vn=0.1; noise=vn*(randn(size(t)));%产生噪音 subplot(323); plot(t,noise); grid on; title('噪音信号'); axis([0 .2e-3 -1 1]); askn=(ask+noise);%调制后加噪 subplot(324); plot(t,askn); axis([0 200e-6 -2 2]); title('加噪后信号'); grid on; %带通滤波 %====================================================================== fBW=40e3; f=[0:3e3:4e5]; w=2*pi*f/fs; z=exp(w*j); BW=2*pi*fBW/fs; a=.8547;%BW=2(1-a)/sqrt(a) p=(j^2*a^2); gain=.135; Hz=gain*(z+1).*(z-1)./(z.^2-(p)); subplot(325); plot(f,abs(Hz)); title('带通滤波器'); grid on; Hz(Hz==0)=10^(8);%avoid log(0) subplot(326); plot(f,20*log10(abs(Hz))); grid on; title('Receiver -3dB Filter Response'); axis([1e5 3e5 -3 1]); %滤波器系数 a=[1 0 0.7305];%[1 0 p] b=[0.135 0 -0.135];%gain*[1 0 -1] faskn=filter(b,a,askn); figure(2) subplot(321); plot(t,faskn); axis([0 100e-6 -2 2]); title('通过带通滤波后输出'); grid on; cm=faskn.*car;%解调 subplot(322); plot(t,cm); axis([0 100e-6 -2 2]); grid on; title('通过相乘器后输出'); %低通滤波器 %================================================================== p=0.72; gain1=0.14;%gain=(1-p)/2 Hz1=gain1*(z+1)./(z-(p)); subplot(323); Hz1(Hz1==0)=10^(-8);%avoid log(0) plot(f,20*log10(abs(Hz1))); grid on; title('LPF -3dB response'); axis([0 5e4 -3 1]); %滤波器系数 a1=[1 -0.72];%(z-(p)) b1=[0.14 0.14];%gain*[1 1] so=filter(b1,a1,cm); so=so*10;%add gain so=so-mean(so);%removes DC component subplot(324);