1.完整项目描述和程序获取
>面包多安全交易平台:https://mbd.pub/o/bread/Y5qVm5hy
>如果链接失效，可以直接打开本站店铺搜索相关店铺：
>如果链接失效，程序调试报错或者项目合作也可以加微信或者QQ联系。
2.部分仿真图预览
3.算法概述
ITS短波宽带信道模型中的最为主要的一个部分就是功率延迟剖面函数，功率延迟剖面函数的确定是信道仿真的最大难点，这是由于功率延迟剖面函数中的需要参数和变量都需要实际测试才能获得。但是，在实际中，我们可以通过近似的方法来获得功率延迟剖面函数，本文，我们将研究ITS模型中功率延迟剖面函数的等效表示方式。通过相关文献可知，功率延迟谱函数的曲线与某些参数先的Nakagami函数曲线比较接近，所以，我们将通过Nakagami函数来模拟出最小误差的ITS模拟功率延迟剖面函数。
4.部分源码
%ITS宽带短波通信信道建模仿真代码
m = 1.1;
delta = 0.5628;
%随机调制函数的仿真
delta = [0.0466 0.0659 0.0932 0.1318];
f = -200:1:200;
for n = 1:length(f)
for i=1:4
SG(i,n) = 1/sqrt(2*pi*delta(i)) * exp((-f(n)^2/2) * delta(i)^2);
end
end
figure;
plot(f,SG(1,:),'r');hold on;
plot(f,SG(2,:),'b');hold on;
plot(f,SG(3,:),'k');hold on;
plot(f,SG(4,:),'g');hold off;
title('Matlab模拟高斯函数形状');
grid on;
legend('delta = 0.0466','delta = 0.0659','delta = 0.0932','delta = 0.1318');
%ITS仿真模型分析
a = 1; %设置仿真幅度
fc = 300; %载波频率
fs = 3000; %信号的采样频率
thta = pi/6; %信号的初始相位
delay = [0 100]; %多径延迟，有几条多径，就输入几个值
fm = [0 200]; %每个多径的频谱展开
fd = [0 10]; %频率偏移
A = [1 3]; %每条多径的幅度
inter = 30; %瑞丽信道参数
SNR = 2; %信噪比
t = 20*pi/10000:20*pi/10000:20*pi;%仿真时间
st = a*cos(2*pi*fc*t+thta); %原始的发送信号
Ns = length(st);
M = length(A);
%希尔伯特变换
n_delay = ceil(delay*10^-3.*fs);
r1 = zeros(1,Ns+max(n_delay));
%ITS宽带短波通信信道建模仿真代码
m = 1.1;
delta = 0.5628;
C = 1;
tao = 2200/length(r1):2200/length(r1):2200;
tao = tao/1000;
for i=1:length(tao)
P(i) = (1/C)*( (sqrt((2*m-1)*delta^2/(2*m*tao(i)))^(1-2*m))...
* exp( -(2*m-1)/2 - (m*tao(i)^2)/(2*delta^2) ) );
end
%通过信道
for m = 1:M
%加入频偏
r1 = r1.*exp(j*(2*pi*fd(m).*[1:length(r1)]/fs));
%加入多径
r1 = r1 + [zeros(1,n_delay(m)),reylei(st,fm(m),fs,Ns,inter)...
.*A(m),zeros(1,max(n_delay)-n_delay(m))];
%加入功率延迟剖面函数
r2 = r1.*sqrt(P);
end
r = r2;
r = real(r)./sqrt(sum(A.^2));%获得平均值
r = awgn(r,SNR,'measured'); %通过高斯信道
01_159m
