x=0:17;
y=x;
a1=[1 0.455 0.346 0.483 0.436 0.5 0.482 0.662];
b1=[0 0.00056 0.00272 0.00003 0.01386 0.00021 0.00062 0.02];
a2=[0 0.0079 0.0154 0.0088 0.0068 0.0089 0.0094 -0.0089];
b2=[0 0.00381 0.00864 0.00546 0.00414 0.00435 0.00544 0.0307];
hold off
for j=1:8
m1(j)=a1(j);
variance1(j)=b1(j);
m2(j)=a2(j);
variance2(j)=b2(j);
for i=1:length(x)
SNR=x(i);
simout=sim('Ka_test');
%simout=sim('Ka');运行simulink中的Ka,上一条指令是运行Ka_test。两个不大相同的仿真框图
y(i)=mean(BitErrorRate);
end
z(j,1:18)=y;
%z(j,:)=y;表示取j行无穷列
end
semilogy(x,z(1,1:18),'k:',x,z(2,1:18),'-r+',x,z(3,1:18),'--yx',x,z(4,1:18),':gs',x,z(5,1:18),'-.bd',x,z(6,1:18),'-k^',x,z(7,1:18),'--mp',x,z(8,1:18),'-.co');
% semilogy(x,z(1,:),'-r+',x,z(2,:),'-yx',x,z(3,:),'-gs');以此类推输出七个
%grid on; 加网格线
set(gca,'xgrid','on');
xlabel('信噪比SNR');
ylabel('误码率Pe');
title('不同天气下幅相均变化时的Pe-SNR仿真');
legend('AWGN','晴天','黑云','小雨','雷雨','小雪','冰霜','中雨');
%多普勒频移影响
x=0:15;
y=x;
FrequencySeparation=24000;
BitRate=10000;
SamplesPerSymbol=2;
Velocity=40;
LightSpeed=3*10^8;
Frequency=825*10^6;
WaveLength=LightSpeed/Frequency;
Fd=Velocity*10^3/3600/WaveLength;
hold off;
for i=1:length(x)
SNR=x(i);
sim('untitled2');
y(i)=mean(BitErrorRate);
End
semilogy(x,y);
grid on;
xlabel('SNR');
ylabel('BER');
title('幅度变化的BER-SNR仿真');
ka.rar_KA_KEQ_Ka Band_Ka频段_matlab simulink
版权申诉
5星 · 超过95%的资源 87 浏览量
2022-09-21
22:05:46
上传
评论
收藏 18KB RAR 举报
alvarocfc
- 粉丝: 109
- 资源: 1万+
评论1