| 查看: 383 | 回复: 0 | |||
[求助]
循环解一个方程将值存入两个数组
|
|
想把这个方程解出来存入两个数组中，但是方程怎么都报错..... 求大神帮助啊...... 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
» 猜你喜欢
读博的好处
已经有3人回复
-
什么时候出结果，有咨询渠道？？？
已经有3人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有26人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有7人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有8人回复
-
分享一下我之前已中青C的计划书的filecode
已经有4人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有3人回复
-
为什么网上很多人说本周 12号出结果
已经有6人回复
-
听说今天filecode变了
已经有49人回复
