Unable to perform assignment because the indices on the left side are not compatible with the size of the right side.

9 次查看（过去 30 天）
AeroEng
AeroEng 2020-3-29
编辑： Adam Danz 2020-3-31
When I go to run this code, it says that my A(i) in the second for loop cannot match up with the size of the right side. Can someone help me debug what is happening here and why they arent compatible.
Thank you!
R0 = [8598.311720315167; 5076.086378182427; 46.47171438826854];
y = 2016; % Year
m = 2; % Month of February
d = 20; % Day
h = 12; % Hour
min = 28; % Minute
s = 25; % Second
tmax = 24*3600; % Seconds
t = s:10:tmax; % Seconds
r_t = R0;
for i = 1:length(t)
d2000(i) = (367*y) - floor((7*(y + floor((m+9)/12)))/4) + floor((275*m)/9)...
+ ((h + (min/60) + ((s+t(i))/3600))/24) + d - 730531.5;
theta(i) = 280.46061837 + 360.98564736628*d2000(i);
end
for i = 1:length(t)
A(i) = [cos(theta(i)) sin(theta(i)) 0;...
-sin(theta(i)) cos(theta(i)) 0;...
0 0 1];
rECEF(i) = A(i)*r_t;
end
xx = rECEF(1,:);
yy = rECEF(2,:);
zz = rECEF(3,:);
a = 6378.137*10^3; % Meters --> Ellipsoidal Equatorial Radius
e = sqrt(0.00669437999); % Eccentricity of Ellipsoid
b = a*sqrt(1-e^2);
w = sqrt(xx^2 + yy^2);
l = e^2/2;
m = (w/a)^2;
n = ((1-e^2)*zz/b)^2;
i = -(2*l^2 + m + n)/2;
k = l^2*(l^2 - m - n);
q = (m + n - 4*l^2)^3/216 + m*n*l^2;
D = sqrt((2*q - m*n*l^2)*m*n*l^2);
beta = i/3 - (q+D)^(1/3) - (q-D)^(1/3);
t = sqrt(sqrt(beta^2-k) - (beta+i)/2) - sign(m-n)*sqrt((beta-i)/2);
w1 = w/(t+l);
z1 = (1-e^2)*zz/(t-l);
lat = atan(z1/((1-e^2)*w1));
long = 2*atan((w-xx)/yy);
alt = sign(t-1+l)*sqrt((w-w1)^2 + (zz-z1)^2);
figure()
plot(long,lat,'r*')
xlabel('Longitudinal')
ylabel('Latitudinal')
grid on
grid minor

请先登录，再回答此问题。

回答（2 个）

MaryD
MaryD 2020-3-29
You're trying to asign 3x3 matrix to A(i) which is vector of size i.
try A(1:3,1:3,i)=[cos(theta(i)) sin(theta(i)) 0;...
-sin(theta(i)) cos(theta(i)) 0;...
0 0 1];
Adam Danz
Adam Danz 2020-3-29
编辑：Adam Danz 2020-3-31
The transformation matrix produces a 3x3 matrix and you're attempting to store the 3x3 matrix into a single 1x1 space in A(i).
It doesn't look like you need to store each iteration of A so you could probably just store the output in "A" instead of "A(i)". However, when you multiply "A" with "rECEF" it will result in a 3x1 vector and you'll have the same error since you're storing that 3x1 vector in a 1x1 space "rECEF(i)".
I recommend the following
rECEF = nan(3, length(t)); %pre-allocate your loop variable
for i = 1:length(t)
A = [cos(theta(i)) sin(theta(i)) 0;...
-sin(theta(i)) cos(theta(i)) 0;...
0 0 1];
rECEF(:,i) = A*r_t; %Store the 3x1 values by column
end
xx = rECEF(1,:);
yy = rECEF(2,:);
zz = rECEF(3,:);

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

产品


版本

R2019a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by