MATLAB数字通信仿真（0/1数据产生，16QAM调制，插值，成型滤波，匹配滤波，采样，16QAM解调，判决，误码率计算）

Simulink是MATLAB最重要的组件之一，它提供一个动态系统建模、仿真和综合分析的集成环境。Simulink已成为国内外高等院校高等数学、数值分析、数字信号处理、自动控制理论以及工程应用等课程的基本教学工具。

《控制系统仿真与CAD》国家精品课课件，其中包括PPT与源代码，PPT中的代码可以直接在PPT中运行，适合自动化领域的学习者学习

本文首先推导了永磁同步电机的数学模型，分析了矢量控制的基本原理，并根据本课题所用电机的实际工作需要选择了最大转矩／电流比控制策略，然后介绍了空间电压矢量脉宽调制(Space Vector Pulse WidthModulation，SVPWM)技术的基本原理与实现方法。在此基础上，本文针对采用传统PI调节的PMSM矢量控制系统在实际应用中控制性能不理想的现象，设计了一种采用改进的电流与速度控制器的矢量控制系统，在MATLAB／Simulink环境中建立了PMSM驱动控制系统的仿真模型并进行了仿真研究，仿真结果分别给出了采用传统PI控制与采用改进电流、速度控制器的PMSM驱动控制系统的输出响应波形，并对二者进行了比较分析，验证了所提出的PMSM矢量控制系统的正确性。

基于MIMO-OFDM可见光通信系统仿真，王朔，张晓光，本文设计了基于MIMO-OFDM可见光通信系统。此外，对可见光通信系统的信道模型进行了详细描述。然后，在相同的信道环境下模拟MIMO-OFDM调

薛定宇 陈阳泉的基于Matlab/Simulink的系统仿真技术与应用（第二版），该书资源（System Simulation Techniques with MATLAB and Simulink ）是英文版的和中文版的第二版是完全对应的。 这里还提供的书中的程序及代码实例，有需要的同学及工程研究人员可以下载参考使用。 只不过它是英文版。

精通MATLAB+Simulink系统仿真(2015代码)

利用proteus进行8通道自动温度检测系统仿真，有完整的程序及仿真电路

Matlab论文：基于Matlab的二进制数字调制系统仿真Matlab论文：基于Matlab的二进制数字调制系统仿真Matlab论文：基于Matlab的二进制数字调制系统仿真Matlab论文：基于Matlab的二进制数字调制系统仿真

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);