我如何修复错误的数组大小不兼容?

3视图(30天)

显示旧的评论

安德烈Khosrow 2023年3月17日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/1930880-how-do-i-fix-the-error-of-arrays-having-incompatible-sizes

编辑: 安德烈Khosrow2023年4月18日

我想计算四阶龙格库塔方法来获得速度和我需要做它通过KC间隔。我试着用地板上但不确定,甚至更改我的错误。我应该使用一个for循环KC间隔吗?

                          嗯= 1;%振幅的振荡流m / s
                         
                          D = 0.1;%气缸的直径
                         
                          m = 50;%的公斤
                         
                          ρ= 1024;%密度公斤/立方米的水
                         
                          K = 200;弹簧刚度总额% N / m
                         
                          c = 100;% d6amping系数
                         
                          Ca = 1;%附加质量系数
                         
                          Cd = 1.8;%阻力系数
                         
                          L = 1;%缸的长度
                         
                          美联社= D * L;%投影面积
                         
                          P =π;
                         
                          md =ρ* P * (D * D) / (4 * L);%附加质量
                         
                          %对第4部分KC间隔
                         
                          KC = 2:2:20;
                         
                          T = (KC。* 0.1) / 1;%时期由于KC
                         
                          Om = (2 * P)。/ T;%ω
                         
                          dt = T / 40;%时间步
                         
                          无损检测= (T / dt) * 10;%的总计算步骤
                         
                          X(1:地板(无损检测+ 1)= 0;
                         
                          V(1:地板(无损检测+ 1)= 0;
                         
                          时间(1:地板(无损检测+ 1)=(0:地板(无损检测))。* dt;
                         
                          为n = 1:无损检测
                         
                          %计算kx1和kv1
                         
                          ta =时间(n);
                         
                          aw =嗯* Om。* (cos (Om。* ta) - ((2/3) * cos (2 * Om。* ta)));%加速度的液体
                         
                          大众=嗯* sin (Om。* ta) -((嗯/ 3)* sin (2 * Om。* ta));%的速度流畅
                         
                          弗吉尼亚州= V (n);
                         
                          Xa = X (n);
                         
                          t1 = Ca * md *啊;
                         
                          t2ρ= 0.5 * * Cd *美联社* abs (Vw-Va) * (Vw-Va);
                         
                          t3 = - k * Xa-c * Va;
                         
                          kx1 = V (n);
                         
                          kv1 = (t1 + t2 + t3) / (m + Ca * md);
                         
                          %计算kx2和kv2
                         
                          助教(n) + dt. / 2 =时间;
                         
                          aw =嗯* Om。* (cos (Om。* ta) - ((2/3) * cos (2 * Om。* ta)));%加速度的液体
                         
                          大众=嗯* sin (Om。* ta) -((嗯/ 3)* sin (2 * Om。* ta));%的速度流畅
                         
                          弗吉尼亚州= V (n) + 0.5 * dt。* kv1;
                         
                          Xa = X (n) + 0.5 * dt。* kx1;
                         
                          t1 = Ca * md *啊;
                         
                          t2ρ= 0.5 * * Cd *美联社* abs (Vw-Va) * (Vw-Va);
                         
                          t3 = - k * Xa-c * Va;
                         
                          kx2 = V (n) + 0.5 * dt。* kv1;
                         
                          kv2 = (t1 + t2 + t3) / (m + Ca * md);
                         
                          %计算kx3和kv3
                         
                          助教(n) + dt. / 2 =时间;
                         
                          aw =嗯* Om * (cos (Om * ta) - ((2/3) * cos (2 * Om * ta)));%加速度的液体
                         
                          大众=嗯* sin (Om * ta) -((嗯/ 3)* sin (2 * Om * ta));%的速度流畅
                         
                          弗吉尼亚州= V (n) + 0.5 * dt * kv2;
                         
                          Xa = X (n) + 0.5 * dt * kx2;
                         
                          t1 = Ca * md *啊;
                         
                          t2ρ= 0.5 * * Cd *美联社* abs (Vw-Va) * (Vw-Va);
                         
                          t3 = - k * Xa-c * Va;
                         
                          kx3 = V (n) + 0.5 * dt * kv2;
                         
                          kv3 = (t1 + t2 + t3) / (m + Ca * md);
                         
                          %计算kx4和kv4
                         
                          助教(n) + dt / 2 =时间;
                         
                          aw =嗯* Om * (cos (Om * ta) - ((2/3) * cos (2 * Om * ta)));%加速度的液体
                         
                          大众=嗯* sin (Om * ta) -((嗯/ 3)* sin (2 * Om * ta));%的速度流畅
                         
                          弗吉尼亚州= V (n) + 0.5 * dt * kv3;
                         
                          Xa = X (n) + 0.5 * dt * kx3;
                         
                          t1 = Ca * md *啊;
                         
                          t2ρ= 0.5 * * Cd *美联社* abs (Vw-Va) * (Vw-Va);
                         
                          t3 = - k * Xa-c * Va;
                         
                          kx4 = V (n) + dt * kv3;
                         
                          kv4 = (t1 + t2 + t3) / (m + Ca * md);
                         
                          %下一步值
                         
                          (n + 1) = X (n) + (dt / 6) * (kx1 + 2 * kx2 + 2 * kx3 + kx4);
                         
                          V (n + 1) = (n) + (dt / 6) * (kv1 + 2 * kv2 + 2 * kv3 + kv4);
                         
                          结束

4评论
显示3年长的评论隐藏3年长的评论

安德烈Khosrow 2023年3月21日

编辑:安德烈Khosrow 2023年3月21日

我明白了,所以有无论如何计算它的for循环,消除了数组imcompatibility吗?还是容易分别计算它?[我。e KC = 2与1的元素,这将使dt固定获得V然后手动改变KC = 4 V和获取等。

登录置评。

在回答这个问题。

答案(1)

沃尔特·罗伯森 2023年3月21日

1
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/1930880-how-do-i-fix-the-error-of-arrays-having-incompatible-sizes answer_1197295

是的,你需要遍历KC值。

你的无损检测是一个矢量值计算KC。这意味着

X(1:地板(无损检测+ 1)= 0;

不是你所希望的方式工作。无损检测是一个矢量不要创建一个行无损检测X为每个不同的值的行是不同的长度不将会发生什么。相反,当你有一个:运营商有一种基本价值作为运营商的参数之一,然后第一个元素是用剩下的忽视,如果你写了

X(1:地板(无损检测(1)+ 1)= 0;

6个评论
显示5年长的评论隐藏5年长的评论

安德烈Khosrow 2023年3月23日

我按你的建议做一个了;

                                   
                                   %的代码简单离岸能源开采的初值问题
                                  
                                   %的设备。
                                  
                                   嗯= 1;%振幅的振荡流m / s
                                  
                                   D = 0.1;%气缸的直径
                                  
                                   m = 50;%的公斤
                                  
                                   ρ= 1024;%密度公斤/立方米的水
                                  
                                   K = 200;弹簧刚度总额% N / m
                                  
                                   c = 100;% d6amping系数
                                  
                                   Ca = 1;%附加质量系数
                                  
                                   Cd = 1.8;%阻力系数
                                  
                                   L = 1;%缸的长度
                                  
                                   美联社= D * L;%投影面积
                                  
                                   P =π;
                                  
                                   md =ρ* P * (D * D) / (4 * L);%附加质量
                                  
                                   %对第4部分KC间隔
                                  
                                   T = 0 (10);
                                  
                                   Om = 0 (10);
                                  
                                   dt = 0 (10);
                                  
                                   无损检测= 0 (10);
                                  
                                   KC = 2:2:20;
                                  
                                   为k = 1:10
                                  
                                   T (k) = (KC (k) * 0.1) / 1;
                                  
                                   Om (k) = (2 * P) / T (k);
                                  
                                   dt (k) = T (k) / 40;
                                  
                                   无损检测(k) = (T (k) / dt (k)) * 10;
                                  
                                   结束
                                  
                                   X(1:无损检测(k) + 1) = 0;
                                  
                                   V(1:无损检测(k) + 1) = 0;
                                  
                                   时间(1:无损检测(k) + 1) =(0:无损检测(k)) * dt (k);
                                  
                                   为n = 1:无损检测(k)
                                  
                                   %计算kx1和kv1
                                  
                                   ta =时间(n);
                                  
                                   aw =嗯* Om (k) * (cos (Om (k) * ta) - ((2/3) * cos (2 * Om (k) * ta)));%加速度的液体
                                  
                                   大众=嗯* sin (Om (k) * ta) -((嗯/ 3)* sin (2 * Om (k) * ta));%的速度流畅
                                  
                                   弗吉尼亚州= V (n);
                                  
                                   Xa = X (n);
                                  
                                   t1 = Ca * md *啊;
                                  
                                   t2ρ= 0.5 * * Cd *美联社* abs (Vw-Va) * (Vw-Va);
                                  
                                   t3 = - k * Xa-c * Va;
                                  
                                   kx1 = V (n);
                                  
                                   kv1 = (t1 + t2 + t3) / (m + Ca * md);
                                  
                                   %计算kx2和kv2
                                  
                                   ta =时间(n) + dt (k) / 2;
                                  
                                   aw =嗯* Om (k) * (cos (Om (k) * ta) - ((2/3) * cos (2 * Om (k) * ta)));%加速度的液体
                                  
                                   大众=嗯* sin (Om (k) * ta) -((嗯/ 3)* sin (2 * Om (k) * ta));%的速度流畅
                                  
                                   弗吉尼亚州= V (n) + 0.5 * dt (k) * kv1;
                                  
                                   Xa = X (n) + 0.5 * dt (k) * kx1;
                                  
                                   t1 = Ca * md *啊;
                                  
                                   t2ρ= 0.5 * * Cd *美联社* abs (Vw-Va) * (Vw-Va);
                                  
                                   t3 = - k * Xa-c * Va;
                                  
                                   kx2 = V (n) + 0.5 * dt (k) * kv1;
                                  
                                   kv2 = (t1 + t2 + t3) / (m + Ca * md);
                                  
                                   %计算kx3和kv3
                                  
                                   ta =时间(n) + dt (k) / 2;
                                  
                                   aw =嗯* Om (k) * (cos (Om (k) * ta) - ((2/3) * cos (2 * Om (k) * ta)));%加速度的液体
                                  
                                   大众=嗯* sin (Om (k) * ta) -((嗯/ 3)* sin (2 * Om (k) * ta));%的速度流畅
                                  
                                   弗吉尼亚州= V (n) + 0.5 * dt (k) * kv2;
                                  
                                   Xa = X (n) + 0.5 * dt (k) * kx2;
                                  
                                   t1 = Ca * md *啊;
                                  
                                   t2ρ= 0.5 * * Cd *美联社* abs (Vw-Va) * (Vw-Va);
                                  
                                   t3 = - k * Xa-c * Va;
                                  
                                   kx3 = V (n) + 0.5 * dt (k) * kv2;
                                  
                                   kv3 = (t1 + t2 + t3) / (m + Ca * md);
                                  
                                   %计算kx4和kv4
                                  
                                   ta =时间(n) + dt (k) / 2;
                                  
                                   aw =嗯* Om (k) * (cos (Om (k) * ta) - ((2/3) * cos (2 * Om (k) * ta)));%加速度的液体
                                  
                                   大众=嗯* sin (Om (k) * ta) -((嗯/ 3)* sin (2 * Om (k) * ta));%的速度流畅
                                  
                                   弗吉尼亚州= V (n) + 0.5 * dt (k) * kv3;
                                  
                                   Xa = X (n) + 0.5 * dt (k) * kx3;
                                  
                                   t1 = Ca * md *啊;
                                  
                                   t2ρ= 0.5 * * Cd *美联社* abs (Vw-Va) * (Vw-Va);
                                  
                                   t3 = - k * Xa-c * Va;
                                  
                                   kx4 = V (n) + dt (k) * kv3;
                                  
                                   kv4 = (t1 + t2 + t3) / (m + Ca * md);
                                  
                                   %下一步值
                                  
                                   (n + 1) = X (n) + (dt (k) / 6) * (kx1 + 2 * kx2 + 2 * kx3 + kx4);
                                  
                                   V (n + 1) = (n) + (dt (k) / 6) * (kv1 + 2 * kv2 + 2 * kv3 + kv4);
                                  
                                   结束

我可以看到循环工作但我还是得到1 X的值(n + 1) & V (n + 1)而不是10。

安德烈Khosrow 2023年4月17日

抱歉迟回复,我得到了我想要(调整原始代码),直到评论下一步的计算值。我得到错误“无法执行任务,因为左右两边有不同数量的元素”。我设法得到kx1十值,kx2 . .kx4但我只需要下一步值计算为每个值kx1, kx2等等。我尝试你的建议但是没有帮助与原代码。

                                   %的代码简单离岸能源开采的初值问题
                                  
                                   %的设备。
                                  
                                   嗯= 1;%振幅的振荡流m / s
                                  
                                   D = 0.1;%气缸的直径
                                  
                                   m = 50;%的公斤
                                  
                                   ρ= 1024;%密度公斤/立方米的水
                                  
                                   K = 200;弹簧刚度总额% N / m
                                  
                                   c = 100;% d6amping系数
                                  
                                   Ca = 1;%附加质量系数
                                  
                                   Cd = 1.8;%阻力系数
                                  
                                   L = 1;%缸的长度
                                  
                                   美联社= D * L;%投影面积
                                  
                                   P =π;
                                  
                                   md =ρ* P * (D * D) / (4 * L);%附加质量
                                  
                                   %对第4部分KC间隔
                                  
                                   T = 0 (10);
                                  
                                   Om = 0 (10);
                                  
                                   dt = 0 (10);
                                  
                                   无损检测= 0 (10);
                                  
                                   KC = 2:2:20;
                                  
                                   为k = 1:10
                                  
                                   T (k) = (KC (k) * 0.1) / 1;
                                  
                                   Om (k) = (2 * P) / T (k);
                                  
                                   dt (k) = T (k) / 40;
                                  
                                   无损检测(k) = (T (k) / dt (k)) * 10;
                                  
                                   结束
                                  
                                   X(1:无损检测(k) + 1) = 0;
                                  
                                   V(1:无损检测(k) + 1) = 0;
                                  
                                   (1:无损检测(k) + 1) = 0;
                                  
                                   时间(1:无损检测(k) + 1) =(0:无损检测(k)) * dt (k);
                                  
                                   为n = 1:无损检测(k)
                                  
                                   %计算kx1和kv1
                                  
                                   ta =时间(n);
                                  
                                   aw =嗯* Om。* (cos (Om。* ta) - ((2/3) * cos (2 * Om。* ta)));%加速度的液体
                                  
                                   大众=嗯* sin (Om。* ta) -((嗯/ 3)* sin (2 * Om。* ta));%的速度流畅
                                  
                                   弗吉尼亚州= V (n);
                                  
                                   Xa = X (n);
                                  
                                   t1 = Ca * md *啊;
                                  
                                   Cd t2 = 0.5 *ρ* *。* abs (Vw-Va)。* (Vw-Va);
                                  
                                   t3 = - k * Xa-c * Va;
                                  
                                   kx1 = V (n);
                                  
                                   kv1 = (t1 + t2 + t3) / (m + Ca * md);
                                  
                                   %计算kx2和kv2
                                  
                                   ta =时间(n) + dt (k) / 2;
                                  
                                   aw =嗯* Om。* (cos (Om。* ta) - ((2/3) * cos (2 * Om。* ta)));%加速度的液体
                                  
                                   大众=嗯* sin (Om。* ta) -((嗯/ 3)* sin (2 * Om。* ta));%的速度流畅
                                  
                                   弗吉尼亚州= V (n) + 0.5 * dt (k) * kv1;
                                  
                                   Xa = X (n) + 0.5 * dt (k) * kx1;
                                  
                                   t1 = Ca * md *啊;
                                  
                                   Cd t2 = 0.5 *ρ* *。* abs (Vw-Va)。* (Vw-Va);
                                  
                                   t3 = - k * Xa-c * Va;
                                  
                                   kx2 = V (n) + 0.5 * dt (k) * kv1;
                                  
                                   kv2 = (t1 + t2 + t3) / (m + Ca * md);
                                  
                                   %计算kx3和kv3
                                  
                                   ta =时间(n) + dt (k) / 2;
                                  
                                   aw =嗯* Om。* (cos (Om。* ta) - ((2/3) * cos (2 * Om。* ta)));%加速度的液体
                                  
                                   大众=嗯* sin (Om。* ta) -((嗯/ 3)* sin (2 * Om。* ta));%的速度流畅
                                  
                                   弗吉尼亚州= V (n) + 0.5 * dt (k) * kv2;
                                  
                                   Xa = X (n) + 0.5 * dt (k) * kx2;
                                  
                                   t1 = Ca * md *啊;
                                  
                                   Cd t2 = 0.5 *ρ* *。* abs (Vw-Va)。* (Vw-Va);
                                  
                                   t3 = - k * Xa-c * Va;
                                  
                                   kx3 = V (n) + 0.5 * dt (k) * kv2;
                                  
                                   kv3 = (t1 + t2 + t3) / (m + Ca * md);
                                  
                                   %计算kx4和kv4
                                  
                                   ta =时间(n) + dt (k) / 2;
                                  
                                   aw =嗯* Om。* (cos (Om。* ta) - ((2/3) * cos (2 * Om。* ta)));%加速度的液体
                                  
                                   大众=嗯* sin (Om。* ta) -((嗯/ 3)* sin (2 * Om。* ta));%的速度流畅
                                  
                                   弗吉尼亚州= V (n) + 0.5 * dt (k) * kv3;
                                  
                                   Xa = X (n) + 0.5 * dt (k) * kx3;
                                  
                                   t1 = Ca * md *啊;
                                  
                                   Cd t2 = 0.5 *ρ* *。* abs (Vw-Va)。* (Vw-Va);
                                  
                                   t3 = - k * Xa-c * Va;
                                  
                                   kx4 = V (n) + dt (k) * kv3;
                                  
                                   kv4 = (t1 + t2 + t3) / (m + Ca * md);
                                  
                                   %下一步值尝试修复错误
                                  
                                   (n + 1) = X (n) + (dt (k) / 6) * (kx1 + 2 * kx2 + 2 * kx3 + kx4);
                                  
                                   V (n + 1) = (n) + (dt (k) / 6) * (kv1 + 2 * kv2 + 2 * kv3 + kv4);
                                  
                                   电源(n + 1) = c * V (n + 1) ^ 2;
                                  
                                   结束
                                  
                                      无法执行任务,因为左右有不同数量的元素。

                                   图(01);
                                  
                                   情节(时间(1:无损检测(k)), X(1:无损检测(k)),“- r”);
                                  
                                   包含(“时间(s)”);% honrizontal轴的标签是“x”
                                  
                                   ylabel (“X (m)”);%纵轴的标签是“T”
                                  
                                   传奇(“RK方法”)
                                  
                                   图(02);
                                  
                                   情节(时间(1:无损检测(k)), V(1:无损检测(k)),“- r”);
                                  
                                   包含(“时间(s)”);% honrizontal轴的标签是“x”
                                  
                                   ylabel (“V (m)”);%纵轴的标签是“T”
                                  
                                   传奇(“RK方法”)
                                  
                                   图(03)
                                  
                                   情节(时间(1:无损检测(k))、电力(1:无损检测(k)),“- r”);
                                  
                                   包含(“时间(s)”);% honrizontal轴的标签是“x”
                                  
                                   ylabel (“权力”);
                                  
                                   传奇(“权力”)

沃尔特·罗伯森 2023年4月17日

Om = 0 (10);

Om 10是一个向量的长度。

                                   aw =嗯* Om。* (cos (Om。* ta) - ((2/3) * cos (2 * Om。* ta)));%加速度的液体
                                  
                                   大众=嗯* sin (Om。* ta) -((嗯/ 3)* sin (2 * Om。* ta));%的速度流畅
                                  
                                   弗吉尼亚州= V (n);

aw定义的向量Om,所以啊是一个矢量。同样,大众。但在这一点上,我们被分配一个标量。

                                   t1 = Ca * md *啊;
                                  
                                   Cd t2 = 0.5 *ρ* *。* abs (Vw-Va)。* (Vw-Va);

aw和大众是向量t1和t2是向量。

kv1 = (t1 + t2 + t3) / (m + Ca * md);

这让kv1向量。

                                   aw =嗯* Om。* (cos (Om。* ta) - ((2/3) * cos (2 * Om。* ta)));%加速度的液体
                                  
                                   大众=嗯* sin (Om。* ta) -((嗯/ 3)* sin (2 * Om。* ta));%的速度流畅

由于未知原因,aw和大众都重新计算。但出于同样的原因,在他们成为向量。

弗吉尼亚州= V (n) + 0.5 * dt (k) * kv1;

弗吉尼亚州也重新计算,但这一次它被定义的向量kv1,现在我们变成了一个向量。

                                   t1 = Ca * md *啊;
                                  
                                   Cd t2 = 0.5 *ρ* *。* abs (Vw-Va)。* (Vw-Va);

t1和t2重新计算,但因为同样的原因之前他们成为向量。之前的差异是一个标量现在Va是一个向量的t1比第一次可能有所不同。

总之,你继续和重新计算变量多次,但t1和t2继续向量

kv3 = (t1 + t2 + t3) / (m + Ca * md);

是一个向量,和几个kx *变量向量。

                                   (n + 1) = X (n) + (dt (k) / 6) * (kx1 + 2 * kx2 + 2 * kx3 + kx4);
                                  
                                   V (n + 1) = (n) + (dt (k) / 6) * (kv1 + 2 * kv2 + 2 * kv3 + kv4);

但随着kx *和kv *变量向量,右边的两个计算向量,但你分配一个标量位置在每种情况下。

安德烈Khosrow 2023年4月18日

编辑:安德烈Khosrow 2023年4月18日

哦我明白了。我将如何解决这些问题呢?所以我就第一个值了,跑过。我将做同样的为其他9值。

                                   (n + 1) = X (n) + (dt (k) / 6) * (kx1 + 2 * kx2 (1:1) + 2 * kx3 (1:1) + kx4 (1:1));
                                  
                                   V (n + 1) = (n) + (dt (k) / 6) * (kv1 (1:1) + 2 * kv2 (1:1) + 2 * kv3 (1:1) + kv4 (1:1));
                                  
                                   电源(n + 1) = c * V (n + 1) ^ 2;