close?all;
clc;
%%?結(jié)論：頻率顛倒的問題出在采樣率?。。。?br/>%?Matlab自帶音頻，例如：chrip、train，是以Fs?=?8192Hz采樣得到的，其信號(hào)總時(shí)間為T?=?N?/?Fs
%?采樣時(shí)間點(diǎn)為t?=?（0:N-1）/Fs，如果用載波對(duì)信號(hào)直接調(diào)制，需要注意以下事項(xiàng)：
%?1、如果不重新采樣，需要使得載波頻率fc?%?2、如果想用更高頻率的載波，必須提高音頻信號(hào)的采樣率。假設(shè)新的采樣率為Fs1，不能單純令t?=?（0:N-1）/Fs1，
%?就認(rèn)為提高了采樣率。我們還需要在音頻時(shí)長(zhǎng)不變的情況下，對(duì)音頻信號(hào)進(jìn)行插值，以模擬采樣率的提高（會(huì)產(chǎn)生新的高頻分量）
%%?音頻信號(hào)
load?chirp;
%?y：音頻信號(hào)，F(xiàn)s：對(duì)音頻信號(hào)的采樣頻率8192
y?=?y‘;
length（y）;
%%?參數(shù)
fc?=?20e3;
Fst?=?4?*?fc;???????????????%?對(duì)載波的采樣頻率

y?=?resample（yFstFs）;?????%?對(duì)音頻信號(hào)進(jìn)行再采樣
Nt?=?length（y）;
tt?=?（0:Nt-1）?/?Fst;
%%?調(diào)制
disp（‘-----調(diào)制-----‘）;
tic;
mode?=?“fm“;
if?mode?==?“am“
????yi?=?（0.02?+?y）?.*?cos（2?*?pi?*?fc?*?tt）;
elseif?mode?==?“ssb“
????yi?=?0.5?*?y?.*?cos（2?*?pi?*?fc?*?tt）?-?0.5?*?imag（hilbert（y））?.*?sin（2?*?pi?*?fc?*?tt）;
elseif?mode?==?“dsb“
????yi?=?y?.*?cos（2?*?pi?*?fc?*?tt）;
elseif?mode?==?“fm“
????%?有時(shí)候，繪圖觀察調(diào)頻信號(hào)，發(fā)現(xiàn)其幅度并非恒定
????%?原因：1、信號(hào)頻率過高（Kf過大），采樣率過低，導(dǎo)致繪圖時(shí)出現(xiàn)了近似
????%?2、Kf過小，信號(hào)有一部分特征調(diào)制到包絡(luò)上
????
????%?kf設(shè)置過大，下面按照帶寬計(jì)算公式，合理計(jì)算kf
????%?dfm?=?kf?*?max（abs（y））?/?（2?*?pi）
????%?B?=?2?*?（dfm?/?Fmax?+?1）?*?F?=?2?*?（dfm?+?Fmax）
????%?應(yīng)滿足B?/?2?????
????Fmax?=?4800;???????????%?chirp信號(hào)的最高頻率大概是5000HzFmax對(duì)應(yīng)音頻信號(hào)最高頻率
????Bfm?=?min（2?*?fcFst?-?2?*?fc）;
????Bfm?=?Bfm?*?0.9;???????%?留下一定的帶寬余量
????dfm?=?Bfm?/?2?-?Fmax;??%?dfm表示由音頻信號(hào)調(diào)頻產(chǎn)生的最大頻偏
????kf?=?2?*?pi?*?dfm?/?max（abs（y））;
????phi?=?cumsum（y）?/?Fst;
????yi?=?cos（2?*?pi?*?fc?*?tt?+?kf?*?phi）;
else
????error（‘No?such?mode...‘）;
end
toc;
%%?信道
disp（‘-----信道-----‘）;
snr?=?40;
out?=?awgn（yisnr‘measured‘）;

%%?濾波、解調(diào)
%?電臺(tái)信號(hào)仿真需不需要解調(diào)呢？
%?假設(shè)接收機(jī)的采樣頻率為Fsr。
%?若接收機(jī)采樣頻率Fsr?>?Fs1，仿真時(shí)需要插值
%?若接收機(jī)采樣頻率Fsr?disp（‘-----解調(diào)-----‘）;
tic;

Fsr?=?1.2?*?Fst;????????????????????????????%?接收機(jī)采樣頻率