| 查看: 339 | 回复: 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人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有7人回复
-
压汞仪和BET测气凝胶孔隙率
已经有4人回复
-
博士申请都是内定的吗？
已经有14人回复
-
谈谈两天一夜的“延安行”
已经有13人回复
-
氨基封端PDMS和HDI反应快速固化
已经有11人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有11人回复
-
论文投稿求助
已经有4人回复
-
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛
已经有3人回复
-
投稿精细化工
已经有6人回复