不工作

15的观点(30天)
Donghun李
Donghun李 2020年5月12日
编辑: Mehmed萨德 2020年5月12日
一个= 0.06;
k_l = 26400;%的线性刚度
m = 483;%的质量
f = @ (t, x, Om, l, k_s, d) [x (2);
(2 * k_s * (x (1) - (A * sin (Om * t)))) *
(√(l-d) ^ 2 + (x (1) - (A * sin (Om * t))) ^ 2) - l) /
(m * (√(l-d) ^ 2 + (x (1) - (A * sin (Om * t))) ^ 2))));
% %
Om_array = linspace(0, 20日,21日);%在rad / s - 1
l_array = linspace (0.2、1、21);
[om_array, L_array] = meshgrid (om_array L_array);
d = linspace (-0.005, -0.03, 10);
% d = 1;
Response_amp = 0([大小(Om_array),元素个数(d)));
T = 150;
x0 = (0,0);
k = 1:元素个数(d)
i = 1:元素个数(Om_array)
j = 1:元素个数(l_array)
Om = om_array (i, j);
l = L_array (i, j);
k_s = - (k_l * (l-d)) / (4 * d);%弹簧刚度
(t, x) =数值(@ (t, x) f (t, x, Om, l, k_s, d (k)), (100 t), x0);
Response_amp (i, j, k) = (max (x (: 1)) - min (x (: 1))) / 2;
结束
结束
结束
% %的阴谋
图(1);
ax =轴();
视图(3);
持有(ax);
视图(33 [30]);
网格
i = 1:尺寸(Response_amp, 3)
网格(om_array /(2 *π)L_array, Response_amp(:,:我));
结束
持有
d = 1;
Response_amp = 0([大小(Om_array),元素个数(d)));
T = 150;
x0 = (0,0);
k = 1:元素个数(d)
i = 1:元素个数(Om_array)
j = 1:元素个数(l_array)
Om = om_array (i, j);
l = L_array (i, j);
k_s = - (k_l * (l-d)) / (4 * d);%弹簧刚度
(t, x) =数值(@ (t, x) f (t, x, Om, l, k_s, d (k)), (100 t), x0);
Response_amp (i, j, k) = (max (x (: 1)) - min (x (: 1))) / 2;
结束
结束
结束
% %的阴谋
图(1)%;
ax =轴();
视图(3);
持有(ax);
视图(33 [30]);
网格
i = 1:尺寸(Response_amp, 3)
网格(om_array /(2 *π)L_array, Response_amp(:,:我));
结束
包含(的频率(赫兹))
ylabel (弹簧(m)的长度)
zlabel (“响应振幅(m)”)
集(gca),“字形大小”15)
持有
嗨,这个代码代表不同的图表取决于d的价值。首先,d将linspace (-0.005, -0.03, 10)。通过坚持,我增加了单值d是1。然而,这段代码给我分开2图而不是合并。我怎么解决这个问题?
感谢你的阅读。

接受的答案

Mehmed萨德
Mehmed萨德 2020年5月12日
编辑:Mehmed萨德 2020年5月12日
ax =轴();
创建新轴电流图。所以叫它一次
(ax)
将轴或推迟。当你调用它第一次添加nextplot,当你叫它再次取代nextpot
例如
图中,
ax =轴;
情节(兰德(10))
(ax)
现在检查如果更换或添加nextplot属性
ax.NextPlot
ans =
“添加”
现在又叫它
(ax)
现在检查Nextplot财产
ax.NextPlot
ans =
“替换”
所以你需要从你的代码注释掉这两条线(情节第二次,第一次是不错的)

更多的答案(0)

社区寻宝

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

开始狩猎!