还是关于matlab的问题!

来源:百度知道 编辑:UC知道 时间:2024/06/15 07:12:43
i=-50;
for ii=1:13;
i=i+30;
if i>90&i<270
x(ii)=pi+atan((650+92*sin(i*pi/180))/(92*cos(i*pi/180)));
else
x(ii)=atan((650+92*sin(i*pi/180))/(92*cos(i*pi/180)));
end;
y=asin((1119-1125*sin(x))/281);
s3(ii)=92*cos(i*pi/180)*(1./cos(x(ii)));
end;
上述可以求出x,y,s3的值,想问一下接着怎样将这些值传递到以下矩阵中并求出a,b,c,d的值
cos(x) -s3*sin(x) 0 0 a sin(i*pi/180)
[ sin(x) s3*cos(x) 0 0 ] [ b ]= cos(i*pi/180)
0 -s3*sin(x) sin(y) -1 c 0
0 s3*cos(x) cos(y) 0 d 0
矩阵不好写,每纵行对齐可以看出来的,不好意思
对就是求你说的X的

你写的真是。。。
A=[cos(x) -s3.*sin(x) 0 0
sin(x) s3.*cos(x) 0 0
0 -s3.*sin(x) sin(y) -1
0 s3.*cos(x) cos(y) 0 ];
b=[sin(i*pi/180);cos(i*pi/180);0;0];
Ax=b,求x,对不对?

i=-50;
for ii=1:13;
i=i+30;
if i>90&i<270
x(ii)=pi+atan((650+92*sin(i*pi/180))/(92*cos(i*pi/180)));
else
x(ii)=atan((650+92*sin(i*pi/180))/(92*cos(i*pi/180)));
end;
y=asin((1119-1125*sin(x))/281);
s3(ii)=92*cos(i*pi/180)*(1./cos(x(ii)));
A=[cos(x(ii)) -s3(ii).*sin(x(ii)) 0 0
sin(x(ii)) s3(ii).*cos(x(ii)) 0 0
0 -s3(ii).*sin(x(ii)) sin(y(ii)) -1
0 s3(ii).*cos(x(ii)) cos(y(ii)) 0 ];
b=[sin(i*pi/180);cos(i*pi/180);0;0];
res=inv(A)*b
end;
你试试看

在想

.不懂!