Matlab高频雷达海洋回波谱的模拟仿真(6页).doc

资源ID：35409575 资源大小：129.50KB 全文页数：6页
资源格式： DOC 下载积分：15金币

快捷下载

会员登录下载

微信登录下载

三方登录下载：

微信扫一扫登录

下载资源需要15金币

邮箱/手机：
温馨提示：	快捷下载时，用户名和密码都是您填写的邮箱或者手机号，方便查询和重复下载（系统自动生成）。如填写123，账号就是123，密码也是123。
支付方式：
验证码：	换一换

账号：
密码：
验证码：	换一换
当日自动登录忘记密码？

友情提示

1、下载资料失败解决办法

2、PDF文件下载后，可能会被浏览器默认打开，此种情况可以点击浏览器菜单，保存网页到桌面，就可以正常下载了。

3、本站不支持迅雷下载，请使用电脑自带的IE浏览器，或者360浏览器、谷歌浏览器下载即可。

4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩，下载后原文更清晰。

5、试题试卷类文档，如果标题没有明确说明有答案则都视为没有答案，请知晓。

网站客服

侵权投诉

Matlab高频雷达海洋回波谱的模拟仿真(6页).doc

-Matlab高频雷达海洋回波谱的模拟仿真-第 6 页%-一、二阶海杂波谱仿真图-%-线性调频连续波发射信号仿真-%clc;clear all;close all;%-初始参数设置-%T = 0.25; % 脉冲宽度(即脉冲重复周期)TT =64; % 相干积累时间N = TT/T; % 相干积累时间内的脉冲数% Fs = B*2; % 采样频率saNum = 5; % 每个周期的采样点数Ts = T/saNum; % 采样间隔sig_ran = 60; % 目标所在的距离单元Fs = 10e2; %-发射信号仿真-%tt = Ts:Ts:T;t = 0:1/(Fs):TT;%-海杂波信号仿真-%-第一种-% 一阶海杂波f0 = 18e6; % 雷达载频g = 9.8; % 重力加速度c = 3e8;lambda = c/f0; % 电磁波波长alfa = 8.1e-3;beita = 0.74;U = 25; % 待变参数，海面上处的风速kc = g/U2 %截止波束theta = 180; % 雷达波束方向到海浪波列前进方向的转角fi = 0; % 雷达波束方向的方位角% theta+fi %表示海浪波列前进方向的方位角fw = 90 %海面风向的方位角k0 = 2*pi/lambda;fd1 = 0.1012 * (f0/1e6)0.5; % 一阶Bragg峰的多普勒频率,存在洋流时偏移fd2 = -fd1; S1 = 0.005*(-2*k0)-4*exp(-0.74*(kc/(2*k0)2);% 两个Bragg峰的幅度不同(随海面风向的方向角)S2 = 0.005*(-2*k0)-4*exp(-0.74*(kc/(2*k0)2);% 两个Bragg峰的幅度不同(随海面风向的方向角)G1_1 = 4/(3*pi)*(cos(180-fw)*pi/180/2)4;G2_1 = 4/(3*pi)*(cos(0-fw)*pi/180/2)4;Bragg1_1 = 26*k04*pi*S1*G1_1*exp(2*pi*1i*fd1.*t).*exp(-(t).2./(22)*0.4)./(2*pi)0.5*0.4); % 正轴一阶Bragg峰Bragg1_2 = 26*k04*pi*S2*G2_1*exp(2*pi*1i*fd2.*t).*exp(-(t).2./(22)*0.4)./(2*pi)0.5*0.4); % 负轴一阶Bragg峰% Noise = 1.2358e-004*randn(1,length(t);%29也行freq = linspace(-Fs/2,Fs/2,length(t);B_1 = fftshift(abs(fft(Bragg1_1+Bragg1_2);% B_11 = B_1/max(B_1);% B_11 = 20*log10(B_11);% figure(1);% plot(freq,B_11); %一阶bragg峰% xlim(-2 2);% 二阶海杂波f4 = linspace(-1.1*fd1,-2*fd1,20);f3 = linspace(-0.9*fd1,-0.1*fd1,80);f1 = linspace(0.1*fd1,0.9*fd1,80);f2 = linspace(1.1*fd1,2*fd1,20);delta_k = k0/40;delta = 0.011 - i*(0.0121);for n = 1:79 k_n = n*delta_k; for i = 1:length(f1) %m=-1,m=1% if k_n-2*k0< k_p < k_n+2*k0 theta = linspace(0,180,181); k_p = (2*pi*f1(i)+(g*k_n)0.5)2/g;% theta_1 = acos(k_p2-k_n2-4*k02)/4*k_n*k0); theta_2 = asin(k_n*sin(k_n*sin(theta.*pi/180)./k_p)-pi; kafang1 = -i*0.5*(k_n+k_p-k_n*k_p*(2*pi*f1(i)2+(2*pi*fd1)2)/(k_n*k_p)0.5/(2*pi*f1(i)2-(2*pi*fd1)2); kafang2 = 0.5*(k_n*k0*-cos(theta.*pi/180).*(k_p*k0*-cos(theta_2.*pi/180)./k03/delta); kafang = (abs(kafang1+kafang2).2*26*k04*pi*n*(delta_k2); S1 = 0.005*(-k_n).-4*exp(-0.74*(kc/(-k_n)2);% 两个Bragg峰的幅度不同(随海面风向的方向角) S2 = 0.005*(k_p).-4*exp(-0.74*(kc/(k_p)2);% 两个Bragg峰的幅度不同(随海面风向的方向角) G1 = 4/(3*pi)*(cos(theta-fw).*pi/180/2).4; G2 = 4/(3*pi)*(cos(theta_2-fw).*pi/180/2).4; X(i) = trapz(kafang.*G1.*G2.*S1.*S2,theta.*pi/180); theta = linspace(-180,0,181); k_p = (2*pi*f1(i)+(g*k_n)0.5)2/g;% theta_1 = acos(k_p2-k_n2-4*k02)/4*k_n*k0); theta_2 = asin(k_n*sin(k_n*sin(theta.*pi/180)./k_p)+pi; kafang1 = -i*0.5*(k_n+k_p-k_n*k_p*(2*pi*f1(i)2+(2*pi*fd1)2)/(k_n*k_p)0.5/(2*pi*f1(i)2-(2*pi*fd1)2); kafang2 = 0.5*(k_n*k0*-cos(theta.*pi/180).*(k_p*k0*-cos(theta_2.*pi/180)./k03/delta); kafang = (abs(kafang1+kafang2).2*26*k04*pi*n*(delta_k2); S1 = 0.005*(-k_n).-4*exp(-0.74*(kc/(-k_n)2);% 两个Bragg峰的幅度不同(随海面风向的方向角) S2 = 0.005*(k_p).-4*exp(-0.74*(kc/(k_p)2);% 两个Bragg峰的幅度不同(随海面风向的方向角) G1 = 4/(3*pi)*(cos(theta-fw).*pi/180/2).4; G2 = 4/(3*pi)*(cos(theta_2-fw).*pi/180/2).4; Y(i) = trapz(kafang.*G1.*G2.*S1.*S2,theta.*pi/180); %m=1,m=-1 theta = linspace(0,180,181); k_p = (2*pi*f1(i)-(g*k_n)0.5)2/g;% theta_1 = acos(k_p2-k_n2-4*k02)/4*k_n*k0); theta_2 = asin(k_n*sin(k_n*sin(theta.*pi/180)./k_p)-pi; kafang1 = -i*0.5*(k_n+k_p-k_n*k_p*(2*pi*f1(i)2+(2*pi*fd1)2)/(k_n*k_p)0.5/(2*pi*f1(i)2-(2*pi*fd1)2); kafang2 = 0.5*(k_n*k0*-cos(theta.*pi/180).*(k_p*k0*-cos(theta_2.*pi/180)./k03/delta); kafang = (abs(kafang1+kafang2).2*26*k04*pi*n*(delta_k2); S1 = 0.005*(k_n).-4*exp(-0.74*(kc/(k_n)2);% 两个Bragg峰的幅度不同(随海面风向的方向角) S2 = 0.005*(-k_p).-4*exp(-0.74*(kc/(-k_p)2);% 两个Bragg峰的幅度不同(随海面风向的方向角) G1 = 4/(3*pi)*(cos(theta-fw).*pi/180/2).4; G2 = 4/(3*pi)*(cos(theta_2-fw).*pi/180/2).4; A(i) = trapz(kafang.*G1.*G2.*S1.*S2,theta.*pi/180); theta = linspace(-180,0,181); k_p = (2*pi*f1(i)-(g*k_n)0.5)2/g;% theta_1 = acos(k_p2-k_n2-4*k02)/4*k_n*k0); theta_2 = asin(k_n*sin(k_n*sin(theta.*pi/180)./k_p)+pi; kafang1 = -i*0.5*(k_n+k_p-k_n*k_p*(2*pi*f1(i)2+(2*pi*fd1)2)/(k_n*k_p)0.5/(2*pi*f1(i)2-(2*pi*fd1)2); kafang2 = 0.5*(k_n*k0*-cos(theta.*pi/180).*(k_p*k0*-cos(theta_2.*pi/180)./k03/delta); kafang = (abs(kafang1+kafang2).2*26*k04*pi*n*(delta_k2); S1 = 0.005*(k_n).-4*exp(-0.74*(kc/(k_n)2);% 两个Bragg峰的幅度不同(随海面风向的方向角) S2 = 0.005*(-k_p).-4*exp(-0.74*(kc/(-k_p)2);% 两个Bragg峰的幅度不同(随海面风向的方向角) G1 = 4/(3*pi)*(cos(theta-fw).*pi/180/2).4; G2 = 4/(3*pi)*(cos(theta_2-fw).*pi/180/2).4; B(i) = trapz(kafang.*G1.*G2.*S1.*S2,theta.*pi/180); Bragg(i,:) = (X(i)+Y(i)+A(i)+B(i)*exp(2*pi*1i*f1(i)*t); end Bragg1(n,:) = sum(Bragg,1);% endendBB2 = sum(Bragg1,1)./10;B_2 = fftshift(abs(fft(BB2);% BB2 = B_2/max(B_2);% BB2 = 20*log10(BB2);% figure(3)% plot(freq,B_2);% xlim(-2 2);for n = 1:79 k_n = n*delta_k; for i = 1:length(f3) %m=1,m=-1% if k_n-2*k0< k_p < k_n+2*k0 theta = linspace(0,180,181); k_p = (2*pi*f1(i)-(g*k_n)0.5)2/g;% theta_1 = acos(k_p2-k_n2-4*k02)/4*k_n*k0); theta_2 = asin(k_n*sin(k_n*sin(theta.*pi/180)./k_p)-pi; kafang1 = -i*0.5*(k_n+k_p-k_n*k_p*(2*pi*f1(i)2+(2*pi*fd1)2)/(k_n*k_p)0.5/(2*pi*f1(i)2-(2*pi*fd1)2); kafang2 = 0.5*(k_n*k0*-cos(theta.*pi/180).*(k_p*k0*-cos(theta_2.*pi/180)./k03/delta); kafang = (abs(kafang1+kafang2).2*26*k04*pi*n*(delta_k2); S1 = 0.005*(k_n).-4*exp(-0.74*(kc/(k_n)2);% ) S2 = 0.005*(-k_p).-4*exp(-0.74*(kc/(-k_p)2);%) G1 = 4/(3*pi)*(cos(theta-fw).*pi/180/2).4; G2 = 4/(3*pi)*(cos(theta_2-fw).*pi/180/2).4; X(i) = trapz(kafang.*G1.*G2.*S1.*S2,theta.*pi/180); BBB = B_1+B_2 +B_3 ;BBB = BBB/max(BBB);BBB = 20*log10(BBB);plot(freq,BBB);xlabel('多普勒频率（Hz）');ylabel('归一化回波功率（dB）');xlim(-2 2);ylim(-80 0)

注意事项

本文（Matlab高频雷达海洋回波谱的模拟仿真(6页).doc）为本站会员（1595****071）主动上传，淘文阁 - 分享文档赚钱的网站仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对上载内容本身不做任何修改或编辑。若此文所含内容侵犯了您的版权或隐私，请立即通知淘文阁 - 分享文档赚钱的网站（点击联系客服），我们立即给予删除！

温馨提示：如果因为网速或其他原因下载失败请重新下载，重复下载不扣分。