MATLAB的答案

避免根开关使用的根源

3视图(30天)
简•菲利普
简•菲利普 2021年9月13日
编辑: 简•菲利普2021年9月13日
你好,我有以下问题。
这是一个非常简单。
首先制定象征性地一些矩阵 。
然后我想学习 。特别是,我很感兴趣 只有 虚轴的零交叉 而做一些清扫的变量矩阵,即
gm2
gm3
到目前为止那么简单。
得到的分子和denumerator了行列式在每个迭代中,然后使用函数根找到自己的根。
但问题是,在代码的最后一部分
为页= 1:尺寸(real_parts_of_poles, 1)
polecross_index = zci (real_parts_of_poles (pp、:));
如果~ isempty (polecross_index)
ω= imag_parts_of_poles (pp、polecross_index);
毫米=(ω,gm2 gm1vec (polecross_index) polecross_index]
结束
结束
当我发现这些穿越事件,我不能继续像这样因为很明显根和他们的订单混合在几乎每个调用的函数。
所以顺利进化的像下面的次要情节
有这些我不能检测所需的穿越事件的轨迹。
问: 有另一种极如何遵循混凝土轨道/零?
R = 1;
L = 2;
C = 2;
gm1 = 0;
gm2 = 0;
gm3 = 0;
信谊年代
gm1vec = 0:0.1:2.5;
gm2vec = 0:0.5:2;
Z1 =南(长度(gm2vec), (gm1vec))
real_parts_of_poles = [];
imag_parts_of_poles = [];
为m = 1:长度(gm2vec)
为n = 1:长度(gm1vec)
gm2 = gm2vec (m);
gm3 = gm1vec (n);
Y = [1 / R + 2 / (s * L), 2 / (s * L), 0, 0, 0, 0, 0, 0;
2 / (s * L), 2 / (s * L) + s * C + 1 / (s * L), 1 / (s * L), 0, 0, 0, gm3, 0;
0,1 / (s * L), s * C + 2 / (s * L), 1 / (s * L), 0, gm2, 0, 0;
0,0,1 / (s * L), s * C + 2 / (s * L), 1 / (s * L) + gm1, 0, 0, 0;
0,0,0,1 / (s * L), s * C + 2 / (s * L), 1 / (s * L), 0, 0;
0,0,0,0,1 / (s * L), s * C + 2 / (s * L), 1 / (s * L), 0;
0,0,0,0,0,1 / (s * L), s * C + 2 / (s * L), 1 / (s * L);
0,0,0,0,0,0,2 / (s * L), 1 / R + 2 / (s * L);
];
dY =检波器(发票(Y));
[N、D] = numden (dY);
%合成(N D s)
sys =特遣部队(sym2poly (N), sym2poly (D));
P =杆(系统)
% % %表面情节13.9.2021
RR =图像放大(根(sym2poly (D)));
rr =实际(根(sym2poly (D)));
real_parts_of_poles =圆([real_parts_of_poles rr], 2);
imag_parts_of_poles =圆([imag_parts_of_poles RR], 2);
OO =圆(abs (RR (RR > = 0)), 2);
如果isempty (OO)
Z1 (m, n) =南;
其他的% & &(长度(OO) < 2)
Z1 (m, n) = min (OO (1));
% counter_first_crossing_first_surface = 1;
%如果~ counter_first_crossing_second_surface & &(长度(OO) = = 2)
% Z2 (m, n) = OO (2);
% counter_first_crossing_second_surface = 1
%结束
%
%如果~ counter_first_crossing_third_surface & &(长度(OO) = = 3)
% Z3 (m, n) = OO (3);
% counter_first_crossing_third_surface = 1
%结束
结束
%
% %
图(1)
次要情节(2,1,1)
情节(真正的根(sym2poly (D))),图像放大(根(sym2poly (D))),“x”,“颜色”、ps {m})
持有在
情节(真正的根(sym2poly (N))),图像放大(根(sym2poly (N))),“o”,“颜色”、ps {m})
网格在
包含(“\再保险(s)美元”,“翻译”,“乳胶”)
ylabel (“\ Im (s)美元”,“翻译”,“乳胶”)
%暂停(0.002)
甘氨胆酸hax =;
hax。TickLabelInterpreter =“乳胶”;
甘氨胆酸ax =;%获取当前轴
斧子。XAxisLocation =“起源”;%将它们设置为零
斧子。YAxisLocation =“起源”;%将它们设置为零
斧子。盒=“关闭”;%关闭周围的框
轴平等的
%
%
%
图(1)
次要情节(2,1,2)
情节(gm1vec (n)、真实(根(sym2poly (D))),“x”,“颜色”、ps {m})
持有在
情节(gm1vec (n)、真实(根(sym2poly (n))),“o”,“颜色”、ps {m})
ylabel (“美元\再保险\ mathrm{侦破}\,\ mathbf {Z} (s) _ {\ Omega_0 / \ Omega_ \ infty} $”,“翻译”,“乳胶”)
包含(“美元g_ {\ mathrm {m} 3}的美元,“翻译”,“乳胶”)
网格在
ylim ([-0.2, 0.1])
结束
% %
% real_parts_of_poles;
% imag_parts_of_poles;
%图
%的阴谋(gm1vec real_parts_of_poles (6,:))
为页= 1:尺寸(real_parts_of_poles, 1)
polecross_index = zci (real_parts_of_poles (pp、:));
如果~ isempty (polecross_index)
ω= imag_parts_of_poles (pp、polecross_index);
毫米=(ω,gm2 gm1vec (polecross_index) polecross_index]
结束
结束
结束

在回答这个问题。

答案(0)

类别

找到更多的在频域分析在帮助中心和文件交换

标签

下载188bet金宝搏


释放

R2019a

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

翻译的