24小时热门版块排行榜    

查看: 359  |  回复: 0

wonderlong

铁虫 (初入文坛)

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

想把这个方程解出来存入两个数组中,但是方程怎么都报错..... 求大神帮助啊......
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=1length(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

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wonderlong 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料学硕297已过四六级求调剂推荐 +9 adaie 2026-03-19 9/450 2026-03-20 21:12 by 无际的草原
[考研] 材料080500调剂求收留 +8 一颗meteor 2026-03-13 8/400 2026-03-20 20:00 by zhukairuo
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +3 solanXXX 2026-03-20 3/150 2026-03-20 16:41 by fxue1114
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
信息提示
请填处理意见