24小时热门版块排行榜

返回列表

wonderlong

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 12
帖子: 5
在线: 43分钟
虫号: 1506306
注册: 2011-11-23
专业: 计算机软件

[求助] 循环解一个方程将值存入两个数组

想把这个方程解出来存入两个数组中，但是方程怎么都报错..... 求大神帮助啊......
load '037.txt';
load 'std.txt'
figure(1);
x1=X037(:,1);
y1=-X037(:,2);

subplot(2,2,1);
R=plot(x1,y1,'*');
title('实测');
% hold on;

x2=std(:,1);
y2=std(:,2);
subplot(2,2,2);
S=plot(x2,y2,'r*');
title('标准');

real_x=x1(find((x1>20)&(x1<50)));  % 查找符合要求的数据坐标
real_xh=find((x1>20)&(x1<50)); %返回范围内横坐标对应序号
real_y=y1(real_xh(1):real_xh(length(real_xh))); %符合要求的纵坐标

for j=1

length(real_xh)-1)
real_dis(j)=sqrt((real_x(j)-real_x(j+1))*(real_x(j)-real_x(j+1))+(real_y(j)-real_y(j+1))*(real_y(j)-real_y(j+1)));
end;

real_dis_max_xh=find(real_dis ==max(real_dis)); %距离最大值起始点位置

k=1;
for t=1:4             %real_D() 从MXH往后4个点距相邻最大距离两点的起点的距离
real_D(k)=sqrt((real_x(real_dis_max_xh)-real_x(real_dis_max_xh+t))*(real_x(real_dis_max_xh)-real_x(real_dis_max_xh+t))+(real_y(real_dis_max_xh)-real_y(real_dis_max_xh+t))*(real_y(real_dis_max_xh)-real_y(real_dis_max_xh+t)));
k=k+1;
end;

% %
% %用距离real_D[]在标准上确定点选点范围为(a,b),此线段在标准上为一直线
a=55.9;b=59.5;

std_xh=find(x2==a);
std_start_x=x2(std_xh);std_start_y=y2(std_xh); %设为标准轨上的起点

%用标准上那段直线和real_D（）来找点
%直线方程：if (x >= 55.83) and (x < 59.58) then
%       result := ((-35.34 + 12.83) / (59.58 - 55.83)) * (x - 55.83) - 12.83;

%双层循环，外层循环变起点，内层循环变距离
X=zeros(36,4);
Y=zeros(36,4);

for i=1:36
for j=1:4
      X=solve(y(i,j)+12.83==(-35.34+12.83)/(59.58-55.83)*(x(i,j)-55.83),real_D(j)*real_D(j)==(y(i,j)-y2(std_xh+i-1))*(y(i,j)-y2(std_xh+i-1))+(x(i,j)-x2(std_xh+i-1))*(x(i,j)-x2(std_xh+i-1)),x(i,j));
      Y=solve(y(i,j)+12.83==(-35.34+12.83)/(59.58-55.83)*(x(i,j)-55.83),real_D(j)*real_D(j)==(y(i,j)-y2(std_xh+i-1))*(y(i,j)-y2(std_xh+i-1))+(x(i,j)-x2(std_xh+i-1))*(x(i,j)-x2(std_xh+i-1)),y(i,j));
end;
end;

问题在于最后的这个解方程这始终有问题，求大神帮助。。
问题程序从  for r=1:36开始
想要两个数组X,Y分别存x,y的值

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : std.txt

2015-07-23 10:53:59, 19.35 K

» 猜你喜欢

1楼 2015-07-23 10:54:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 wonderlong 的主题更新

返回列表