matlab通信模擬pdf

⑴ 《通信天線建模與MATLAB模擬分析》中文電子書

http://www.gougou.com/search?search=%E9%80%9A%E4%BF%A1%E5%A4%A9%E7%BA%BF%E5%BB%BA%E6%A8%A1%E4%B8%8EMATLAB%E4%BB%BF%E7%9C%9F%E5%88%86%E6%9E%90&id=1

全部在這里..

⑵ 通信原理課程設計 基於Matlab的通信系統模擬 －模擬調制系統

AM調制的代碼：

% f0; 載波頻率 (Hz)
% fm; 單(多)頻調制頻率點 (Hz)
% B: 實際基帶帶寬 %==================
% MaB; 最大帶寬 (Hz) %==================
% fs; (中頻)采樣頻率 (Hz)
% fs1: 第一采樣頻率
% fs2_expect 預期的(第二)采樣頻率 (Hz)
% T ; 模擬時間 (s)
% beta; 調幅系數
% Kv ; VCO靈敏度 (rad/s/V)
% n0; 雜訊功率譜密度 (W/kHz)
function AM_Simulation(f0,DF,fm,B,MaxB,fs,fs1,fs2_expect,T,beta,K1,K2,Kv ,n0)

%====================================================================================

Ts = 1/fs; % 采樣間隔
N = floor(T/Ts); % 時域采樣總點數 % （基於fs1）
t = 0:Ts:(N-1)*Ts; % 時間采樣

%-------------------------------------------------------------
% 計算以第一采樣頻率fs1欠采以後的載波頻率f01
f01 = f0 - floor(f0/fs1)*fs1;
% if(f01>fs1/2)
% f01 = fs1 - f01;
% end;
%-------------------------------------------------------------

% 計算半帶濾波的級數
K = floor(log2(fs1/fs2_expect)); % fs2_expect為預期的第二采樣頻率
fs2 = fs1/(2^K); % fs2 為實際的第二采樣頻率

%=========================================================================================================================================
% 調制信息

M = length(fm);
am = zeros(1,N);
for m=1:M
am = am + cos(2*pi*fm(m)*t);
end;
dc_offset = max(am)/beta;
a = dc_offset + am ;
%----------------------
am = am/dc_offset;
a = a/dc_offset;%%%%%%
%----------------------

%=========================================================================================================================================
% 調幅

theta0 = 2*pi*rand ; % 載波的隨機初始相位
s = a.*cos( 2*pi * (f0+DF) * t + theta0); % 已調幅信號

%=========================================================================================================================================
% 雜訊信道

noise = GenerateNoise(n0,fs,N,(f0+DF));

r = s + noise;

%-----------------------------------------------------------------------------------------------------------------------------------------
% r_downRate = RateAdjust(r,fs, fs1, 'ItpFlt'); %%%%%% 注意：若對白雜訊（或帶限白雜訊）進行欠采樣則將導致雜訊功率譜密度的提升。
%-----------------------------------------------------------------------------------------------------------------------------------------

%=========================================================================================================================================
% 接收機帶通濾波

fc = MaxB;
% r_BPF = BPF_Filter(r, (f0+DF),fc,fs);
[r_BPF,hBPF] = BPF_Filter(r, f0,fc,fs); % 接收機預先知道f0，但預先不知道DF

%----------------------------------------------------------------------------------------------
%theta_BPF = Phi_BPF(f0,fc,fs, (f0+DF));
theta_BPF = Phi_BPF(hBPF ,fs, (f0+DF)); % 該參數不用於解碼，但在繪制輸入鎖相環的theta1時需要知道該參數

%=========================================================================================================================================
% 對BPF輸出信號進行速率調整

r_BPF_downRate = RateAdjust(r_BPF,fs, fs1, 'ItpFlt');% 選擇fs1需確保對帶限白雜訊進行的不是欠采樣

%---------------------------------------------------------------
Ts1 = 1/fs1; % 采樣間隔
N1 = length(r_BPF_downRate); % 時域采樣總點數 % （基於fs1）
t1 = 0:Ts1:(N1-1)*Ts1; % 時間采樣

%==============================================================================================================================================
% 以預先已知的f01從AM信號中得到I,Q信號

yc = r_BPF_downRate .* (2*cos(2*pi*f01*t1)); %%%%%%
ys = r_BPF_downRate .* (-2*sin(2*pi*f01*t1)); %%%%%%

zI = HalfBandFilters(yc,fs1, K);%%
zQ = HalfBandFilters(ys,fs1, K);%%

[xI, IDelay] = LPF_Filter(zI, B, fs2);
[xQ, QDelay] = LPF_Filter(zQ, B, fs2);

%----------------------------------------------------------
Ts2 = 1/fs2; % 采樣間隔
N2 = length(xI); % 時域采樣總點數 % （基於fs2）
t2 = 0:Ts2:(N2-1)*Ts2; % 時間采樣
%----------------------------------------------------------

thetaIQ = 2*pi*DF*(15*(2^(K)-1)/fs1) + 2*pi*DF* IDelay /fs2; % 該參數不用於解碼，但在繪制輸入鎖相環的theta1時需要知道該參數
%----------------------------------------------------------------------------------------------------------

%==============================================================================================================================================

% 在基帶上用鎖相環對頻差以及初相進行跟蹤以及補償

% g = mean(abs(xI+j*xQ));
% xI = xI/g; % 調整輸入鎖相環的信號幅度, 需改為AGC
% xQ = xQ/g;

% 基帶鎖相
[theta2, Io,Qo, PDo,LFo] = PLL_Baseband(xI,xQ, fs2, K1,K2,Kv, 1); % 1, Carrier

a_demod = Io;
%==============================================================================================================================================
% 去直流
% dc = mean(a_demod(floor(0.2*N2):N2)); % Temp
% a_demod_offsetted = a_demod - dc;

fcL=60;%Hz
a_demod_offsetted = AHPF_Filter(a_demod,fs2, fcL);

%==============================================================================================================================================
% TxRx
figure(1);
subplot(4,2,1);plot(t,a);ylabel('a(t)');
subplot(4,2,3);plot(t,s);ylabel('s_A_M(t)'); axis([min(t) max(t) -inf Inf]);
subplot(4,2,5);plot(t,r);ylabel('r(t)'); axis([min(t) max(t) -inf Inf]);
subplot(4,2,7);plot(t,r_BPF);ylabel('r_B_P_F(t)'); axis([min(t) max(t) -inf Inf]);
xlabel('t (s)')
N_FFT = 2^(floor(log2(fs/B*100)));
subplot(4,2,2);Psd_plot_log_N(a,fs,N_FFT);
subplot(4,2,4);Psd_plot_log_N(s,fs,N_FFT);
subplot(4,2,6);Psd_plot_log_N(r,fs,N_FFT);
subplot(4,2,8);Psd_plot_log_N(r_BPF,fs,N_FFT);

% PLL
figure(2);
plot(t2,PDo);hold on;plot(t2,LFo,'r');hold off;axis([min(t2) max(t2) -inf inf]);grid on;legend('PD\o','LF\o'); hold off;
xlabel('t (s)')

figure(3);
theta1 = theta0 + theta_BPF - thetaIQ + 2*pi * DF * t2; % 基於fs2
theta2 = theta2;
subplot(311);plot(t2,theta1);hold on;plot(t2,theta2,'r');hold off;
legend('\theta_1','\theta_2');axis([min(t2) max(t2) -inf Inf]); grid on; title('Baseband Carrier Synchronization(2)');
theta_e = theta1-theta2;
subplot(312);plot(t2,theta_e);
ylabel('\theta_e(t)');legend('\theta_e');axis([min(t2) max(t2) min(theta_e)-1 max(theta_e)+1]); grid on;
theta_e = angle(exp(j*theta_e));
subplot(313);plot(t2,theta_e);
xlabel('t (s)'),ylabel('\theta_e(t)');legend('\theta_e (-\pi, \pi] ');axis([min(t2) max(t2) min(theta_e)-1 max(theta_e)+1]); grid on;

% TxRx
figure(4);
subplot(611);plot(t,a);xlabel(''),ylabel('a(t)');
subplot(612);plot(t,s);xlabel(''),ylabel('s_A_M(t)');
subplot(613);plot(t,r);xlabel(''),ylabel('r(t)');
subplot(614);plot(t,r_BPF);xlabel(''),ylabel('r_B_P_F(t)');
subplot(615);plot(t2,a_demod);xlabel('t'),ylabel('a_d(t)');
subplot(616);plot(t2,a_demod_offsetted);xlabel('t'),ylabel('m(t)');

% TxRx
figure(5);
subplot(211);plot(t,a);xlabel(''),ylabel('a(t)');
subplot(212);plot(t2,a_demod_offsetted);xlabel('t'),ylabel('m(t)');

figure(4);

%==============================================================================================================================================

⑶ 求 MATLAB/Simulink通信系統建模與模擬實例精講pdf，作者：作 者：邵佳，董辰輝

這本書非常難找，目前能找到的幾本相關圖書如下：

1、MATLAB/Simulink建模與模擬實例精講，張德豐等，機械工業出版社

這本書在新浪ishare有下載，分兩部分，掃描版，含封面421頁，一共53M左右。網路知道裡面關於樓主求助的書有些是通過郵件發的，內容不詳；公開的全部誤為本書。

2、MATLAB通信模擬及應用實例詳解，鄧華，人民郵電出版社

我手裡有這本書的掃描版（PDF帶書簽），427頁，17M左右。網上應該也不難找。

3、MATLAB/SIMULINK 通信系統建模與模擬實例分析，邵玉斌，清華大學出版社

網上流傳的此書不是掃描版，而是提交出版社的書稿。從寫作風格看，我覺得這本書的作者是一個比較嚴肅的學者。大約12M，特別推薦。

4、關於樓主要找的這本書，到目前為止，我只找到11和13章，合計共90頁，約9M。

除非有什麼原因必須使用邵佳的那本書（譬如，學習某課程，老師布置作業是按照該書來的），否則，如果只是想學習相關的知識，個人認為有上面三本書足夠了。而從另一個角度說，邵佳的那本書在網上沒有廣泛流傳，是否也從側面說明該書質量一般？因為，在網上，好書從來都不難找。

作為重點推薦，我先上傳邵玉斌的那本書，樓主看看其它還需要哪本，追問之後我再上傳。