matlab浮点数求绝对值_MATLAB仿真阵列天线切比雪夫综合法(附代码)
来源:cnblogs
在《自适应天线与相控阵》这门课中,了解到了关于理想低副瓣阵列设计的一些方法,其中切比雪夫等副瓣阵列设计方法是一种基础的方法,故将其设计流程写成maltab程序供以后学习使用。在此分享一下。 此方法全称为道尔夫-切比雪夫综合法,简称为切比雪夫综合法,是一种工程实际中常用的可控制副瓣电平的阵列天线综合方法。切比雪夫阵列的特点是:
(1)等副瓣电平;
(2)在相同副瓣电平和相同阵列长度下主瓣最窄,为最佳阵列;
(3)单元数过多时,阵列两端单元激励幅度跳变大,使馈电困难。一般在雷达系统中,为了使其具有较高的抗干扰、抗反辐射导弹的能力,往往要求雷达天线的副瓣尽量低,而采用道尔夫-切比雪夫综合法以及进一步的泰勒综合法等设计的阵列天线就可以实现低副瓣。最早,道尔夫(C.L.Dolph)利用切比雪夫函数来逼近天线阵列的阵因子函数,得到了这种严谨规范的综合方法。而且,经过前人研究,当天线单元N≤13时,切比雪夫阵列从中间到两端的激励分布是单调减小的;而当N>13时,阵列两端单元的激励开始出现跳变。所以对于大型阵列来说一般不宜采用切比雪夫方法综合阵列。所以下面的Matlab程序正常工作在天线单元数N为3到13这个范围内。关于如何采用切比雪夫多项式去设计阵因子的具体技术步骤,另一篇文章较为详细地介绍了,此处不再赘述,大家可以在文尾或评论区查看。下面是可以综合设计天线单元从3到13单元的切比雪夫综合法的Matlab程序:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 |
|
下面即为一个示例:单元间距d=0.6λ、单元数13、主副瓣电平比26dB、扫描角度80度(相对于单元排布方向)的切比雪夫阵列设计。归一化单元电流幅度比为:0.406 0.432 0.604 0.770 0.908 1.000 0.516 1.000 0.908 0.770 0.604 0.432 0.406
----END上文提到的另一篇文章。
阵列天线综合之切比雪夫低副瓣阵列设计 MATLAB(作者:OLIVERMAHOUT)
相控阵天线中,直线阵列作为重要的一种,有着极为广泛的应用。切比雪夫低副瓣阵列设计是一种典型的设计方法。
切比雪夫方法主要是实现低副瓣、窄波束:
其产生的核心如下:
我的理解:因为能量守恒,所有副瓣都一样的时候,能量会更多的集中在副瓣中,
主瓣最大增益也不会改变,这样就可以使主瓣窄,副瓣电平降低。G=4πS/λ2
结合切比雪夫函数,可以得到:
当具体应用时,解决方案如下:
话不多说,其Matlab中的程序如下:
1
% 2019-11% 切比雪夫低副瓣阵列馈电设计_1.0 (端射阵)
close all;clear% digits(3);
% 参数设置lamda = 1; % 波长d = lamda * 0.6; % d为阵元间距theta0 = (120/180)*pi; % 扫描角度theta = 0: 0.01 : pi; % Θ为方向角u = pi*d*(cos(theta)-cos(theta0))/lamda; %T = Chebyshev; % T为切比雪夫恒等式系数矩阵N = 10; % N为直线阵的阵元数量,M为一侧的单元数(对称)R0dB = 26; % R0dB为副瓣电平
if (mod(N,2)==0)M = N / 2;parity = 0; % parity为奇偶性,0为偶数elseM = (N+1)/2;parity = 1;end
% 导入切比雪夫多项式syms x; T = [1;x;2*x^2-1;4*x^3-3*x;8*x^4-8*x^2+1;16*x^5-20*x^3+5*x;32*x^6-48*x^4+18*x^2-1;64*x^7-112*x^5+56*x^3-7*x;128*x^8-256*x^6+160*x^4-32*x^2+1;256*x^9-576*x^7+432*x^5-120*x^3+9*x;512*x^10-1280*x^8+1120*x^6-400*x^4+50*x^2-1];
% 换算副瓣电平R0R0 = 10 ^ (R0dB / 20);
% 计算x0x0 = ((R0 + sqrt(R0^2 -1))^(1/(N-1)) + (R0 - sqrt(R0^2 -1))^(1/(N-1))) * 1/2;
% 定义馈电幅度矩阵II = sym('I', [1 M]);
% 计算展开的方向图表达式S = T(2) * I(1);
for k = 2 : MS = S + T(2*k) * I(k);end
%collect(S,x)%vpa(S)
S_po = coeffs(S,x); % 含电流的方向图多项式系数T_po = sym2poly(T(N)); % 标准的方向图多项式系数(反向了)T_PO = zeros(1,M); for k = 1 : MT_PO(k) = T_po(2*k-1);S_po(k) = S_po(k)/x0^(2*k-1);end% T_PO% vpa(S_po)
% 系数比较求出电流大小eq = sym('eq',[M 1]); % 系数比较恒等式for k = 1 : Meq(k) = S_po(k) == T_PO(M+1-k);end
vpa(eq)I_st = solve(eq);I_ce = struct2cell(I_st);i = zeros(M,1); % 最终的电流矩阵for k = 1 : Mi(k) = I_ce{k,1};i(k) = i(k);endfor k = 2 : Mi(k) = i(k)/i(1); % 电流归一化endi(1) = 1; ii=[1;0.89;0.706;0.485;0.357]; % 用来检验的数据
% 计算最终的阵因子S_all = zeros(1,length(theta));for k = 1 : MS_all = S_all + i(k)*cos((2*k-1)*u);endSS = S_all;
% 画图 —— 直角坐标系S_max = max(S_all); % 归一化处理S_all = 20*log10(abs(S_all/S_max)); % 取分贝值figure('NumberTitle', 'off', 'Name', 'S Parameter (dB) - Cartesian');theta_ = theta * 180 / pi;plot(theta_,S_all,'k','LineWidth',1.5);grid offxlabel('\theta (°)','FontSize',13);ylabel('|S| dB','FontSize',12);axis([0 182 -50 2]);box on
% 画图 —— 极坐标系figure('NumberTitle', 'off', 'Name', 'S Parameter (dB) - Polar');S_pol = SS / max(SS);polarplot(theta,S_all,'k','LineWidth',1.5);thetalim([0 180]);rmin = min(S_all);rmax = max(S_all);rlim([-50 rmax]);
上述测试的N=10的10个阵列,侧射阵(θ=0),副瓣电平SLL=26dB,结果如下:
经过比较,结果较为标准。
更改一下theta0的值,改为120读,即偏离法相30度:
来源:cnblogs
matlab浮点数求绝对值_MATLAB仿真阵列天线切比雪夫综合法(附代码)相关推荐
- matlab浮点数求绝对值_matlab绝对值函数怎么表示 matlab绝对值怎么表示?
matlab绝对值怎么表示? matlab中绝对值是用abs()函数来取的,对于正数来说,返回值为自己:对负数来说,返回值为相反数. 下面是一些例子. 比如求整数的绝对值: a=abs(10) a?= ...
- matlab浮点数求绝对值_在matlab编程中遇到绝对值符号怎么办
fun=@(x,y)y.*log10(y)-x创建关于x,y的匿名函数 log10(x)由于Matlab中自然对数的符号是log(与我们数学上的不同),以10为底的对数符号是log10,同理有log2 ...
- c语言实现台球位置计算,台球碰撞 (C语言代码)三角函数、浮点数求绝对值函数...
解题思路:复合运动,将运动分解为水平方向和垂直方向两个分运动,分别求出横坐标和纵坐标. 注意事项:fabs()是用于浮点数求绝对值:三角函数的输入值为弧度制数:π=acos(-1)--反三角函数求pi ...
- 相控阵天线(六):直线阵列天线特殊综合方法(变形泰勒综合法、贝利斯综合法、伍德沃德抽样法)
目录 简介 变形泰勒综合法 贝利斯综合法 伍德沃德-劳森抽样法 配相抵消法 简介 阵列天线的综合问题是其分析的逆问题,即是在预先给定辐射特性(如方向图形状.副瓣电平等)的情况下,综合出阵列激励幅度和相 ...
- matlab中求绝对值的函数图像,ex的图像(绝对值的函数图像口诀)
由电场E-X图像可以得到1.纵轴表示电场强度的大小.方向2.横轴表示距离3.根据图像可以确定电场强度随距离的变化情况,斜率k=0是匀强电场,k≠0是变化的电场4.. 你这是什么意思? 是e*x还是e^ ...
- matlab 离散点求导_Matlab的离散点曲线导数曲率数值模拟方法
模拟方法: 复制内容到剪贴板 代码 : clc; clear all; close all; x0 = linspace(0, 1); y0 = sin(x0).*cos(x0); h = abs(d ...
- hfss matlab联合仿真_一文搞定matlab 与 STK 联合调试仿真环境配置(文末附软件下载链接)...
最近在做导师给的课题:卫星星座的快速优化设计. 需要用到matlab 和 STK 来进行联合调试仿真,但是这第一步的环境配置就让我头疼了几天.在好几次重装,失败和查找资料之后,我终于成功实现了matl ...
- matlab 平滑曲线连接_平滑轨迹插值方法之多项式插值(附代码)
前言 今天我们来聊聊轨迹插值,在机器人的运动规划和控制领域,参考轨迹的生成是一个历史悠久的问题,已经发展出了一系列的方法.今天我们就来聊一聊轨迹插值领域中最常见的轨迹插值方法:多项式插值. 说明:本文 ...
- Matlab 数字滤波器设计大报告(数字信号处理课程设计)附代码
文章中滤波器设计与信号产生.处理的全部代码,几乎每行都有注释,在这里下载:https://download.csdn.net/download/weixin_42845306/18877388 摘要 ...
最新文章
- MySQL Online DDL的改进与应用
- 【Django】模型层说明
- 一次公司内部的Tech Talk中涉及到的关于语言的发展问题
- assign, retain, copy, weak, strong
- 【优化SQL Server循环更新、插入耗时长的问题】
- 计算文本相似度_Python文本相似性计算
- 文件描述符和fcntl及阻塞非阻塞
- Linux进程全解2——进程环境(环境变量、进程运行的虚拟地址空间)
- 信息学奥赛一本通 2034:【例5.1】反序输出
- KindEditor在eclipse里的配置方法
- 用 radial-gradient 实现波浪效果
- VirtualBox 安装 Linux 时遇到的问题及解决办法
- Linux下源码安装ElasticResearch
- Windows通用文件关联错误恢复
- 【SLAM小车】硬件_TB6612FNG设计记录
- 西安交大计算机专业考研复试,2018西安交通大学计算机考研复试分数线_计算机考研分数线...
- OpenHarmony更新编译问题及解决办法
- 2020-01-04
- 按Enter键调用登录按钮
- HYS与七夕节大作战 vijos