不兼容的数组大小2 x1列向量

1视图(30天)

显示旧的评论

蒂莫西·邓宁 2022年12月2日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/1869452-incompatible-array-sizes-for-2x1-column-vectors

回答: 蒂莫西·邓宁2022年12月7日

答:接受蒂莫西·邓宁

我有一个数组,z, 2行和迭代扩展的列第一个分配: z (: 1) = (0, 2500);

我通过当前(n个数)值数组的z = ivpsolver (z (:, n)),然后再次与z (:, n + 1) = rungekutta (z (:, n));

在这个函数,我得到一个错误行说:znext = z + (1/6) * (Az + 2 * +商务+ 2 * Cz + Dz);阿兹,Bz Cz和Dz 1 x2列向量(z应该)。错误状态“数组已经为这个操作不兼容的大小。”

我缺少一些明显的或错误可能别的地方吗?

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

蒂莫西·邓宁 2022年12月2日

下面是函数,决定就把整件事情(不知道这是坏etiquete,抱歉)。第一个函数应该被称为:

[z n t, q] = ivpSolver (1 [0, 2500] (0,0), 0.25, 90)

函数[z n t, q] = ivpSolver (n, t, z, q, dt,θ)

% ivpSolver解决一个初值问题(IVP)和阴谋的结果

% (T, Z) = ivpSolver (T0, Z0, DT, TE)计算IVP解决方案使用步骤

%大小DT,开始时刻T0和初始状态Z0结束时间

%倾向。解决方案是输出的时间T向量和矩阵的状态

% Z向量。

%设置初始条件

(1)= 25000;

%继续走,直到超过结束时间

而 z (1, n) > = 0 | | z (2 n) > 0

%增量向量由一个时间步的时间

t (n + 1) = t (n) + dt;

% RK4方法申请一次一步

[z: n + 1)、问(:,n + 1)] = stepRungeKutta (z (:, n)、问(:,n), dt, M (n),θ);

如果 0.5米(n) >

M (n + 1) = M (n) - (dM * dt);

其他的

M (n + 1) = M (n);

结束

n = n + 1;

结束

函数 [znext, qnext] = stepRungeKutta (z, q, dt, M,θ)

%计算导数从当前状态

(dz, dq) = stateDeriv2 (z, q, M,θ);

%计算单个近似1日,2日,3日和4日的顺序

%的平均

阿兹= dt * (dz);

Bz = dt * (dz + (Az / 2));

Cz = dt * (dz + (Bz / 2));

Dz = dt * (Dz + Cz);

Aq = dt * (dq);

Bq = dt * (dq + (Aq / 2));

Cq = dt * (dq + (Bq / 2));

Dq = dt * (Dq + Cq);

%结合未来价值平均近似

znext = z + (1/6) * (Az + 2 * +商务+ 2 * Cz + Dz);

qnext = q + (1/6) * (Aq + 2 * Bq + 2 * Cq + Dq);

函数 (dz, dq) = stateDeriv2 (z, q, m,θ)

%计算状态微分系统质-弹

% DZ = stateDeriv (T, Z)计算导数DZ = [V;的)

%状态向量Z = [X;V V),其中X是位移,速度,

%,是加速度。

Cp = 1.2;

Cc = 0.1;

Ac =π;

美联社= 6 * 4;

G = 6.674 e-11;

M = 5.972 e24;

R = 6378年e3;

ρ= atmosEarth (z);

θ=θ*π/ 180;

如果 z (2) > 0 | | z (: 1) > 0

dz1 z = (2);

时= (- m * G / (R + z(1:)) ^ 2)——(罪(θ)* (Cc * Ac + Cp *美联社)*ρ* (z(2:) ^ 2)) /(2 *米);

dq1 q = (2:);

dq2 = (cos(θ)* (Cc * Ac + Cp *美联社)*ρ* (q(2:) ^ 2)) /(2 *米);

其他的

dz1 z = (2);

时= (- m * G / (R + z(1:)) ^ 2)——(罪(θ)* Cc * Ac *ρ* (z(2:) ^ 2)) /(2 *米);

dq1 q = (2:);

dq2 = (cos(θ)* Cc * Ac *ρ* (q(2:) ^ 2)) /(2 *米);

结束

dz = [dz1;时);

dq = [dq1;dq2];

Askic V 2022年12月7日

首先你的第一个函数调用以下方法:

[z n t, q] = ivpSolver (1 [0, 2500] (0,0), 0.25, 90)

但是函数的定义是这样的:

函数[z n t, q] = ivpSolver (n, t, z, q, dt,θ)

功能预计6输入参数调用时,你打个电话5参数。

当函数开始执行时,它将有以下变量在它的范围:

因此,θ的定义。后来,在ivpSolver函数有一行代码:

                                [z: n + 1)、问(:,n + 1)] = stepRungeKutta (z (:, n)、问(:,n), dt, M (n),θ);
                               

你在哪里调用另一个函数提供θ作为输入参数,这将产生第一个错误。

一件事,如果你有一个列,是很不寻常的 ”,因为在所有行第1列声明”。资料片列元素和它的大小,所以我通常的预期

z (: 1) = [0;2500]%,而不是z (: 1) = [0, 2500)

托马斯。 2022年12月7日

只是想添加,您可以使用函数 isrow 和 iscolumn 验证什么MATLAB分别定义为行和列向量。

登录置评。

在回答这个问题。

接受的答案

蒂莫西·邓宁 2022年12月7日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/1869452-incompatible-array-sizes-for-2x1-column-vectors answer_1122612

谢谢,你都是对的,我经过错了其他地方的变量,但我也有另一个问题是导致数组,z, 1 x3代替1 x2

0评论
显示1年长的评论藏1年长的评论

登录置评。

类别

MATLAB 编程

找到更多的在编程在帮助中心和文件交换

s manbetx 845

MATLAB

释放

R2022a

社区寻宝

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

开始狩猎!