I get an error in the loop

2017 3 29
1 个回答
回答已采纳
更新时间：2017 3 30
5 次查看（30 天）

0 个投票

K=[27 27320.5963 4223312.4563
32 25496.3840 4221760.5030
34 28874.9170 4221756.7649
39 27284.2609 4220235.3943
133 26812.2676 4221811.7328]
olc=[0.00000
70.76351
105.09584
0.00000
42.40162
89.83755
0.00000
50.28006
98.61182
0
36.51504
106.45528
0
80.90273
160.68555
276.73136]
dn=[27;27;27;32;32;32;34;34;34;39;39;39;133;133;133;133]
bn=[34;133;32;27;133;39;39;133;27;32;133;34;27;34;39;32]
dog_olc=[dn bn olc]
R=6373394
for i = 1:length(dog_olc)
dn=dog_olc(i,1); bn=dog_olc(i,2); olc=dog_olc(i,3);
xi=K(find(K(:,1)==dn,2)); yi=K(find(K(:,1)==dn,3));
xk=K(find(K(:,1)==bn,2)); yk=K(find(K(:,1)==bn,3));
dt=((-200/pi)/(6*R^2))*((xk-xi)*(2*yi+yk));
Dog_olc(i)=dog_olc(3,i)+dt
end
There are results visually. I do not succeed in getting these results.
http://i.hizliresim.com/7qyVGa.jpg

请先登录，再回答此问题。

 采纳的回答

KSSV
KSSV 2017-3-30

1 个投票

clc; clear all ;
K=[27 27320.5963 4223312.4563
32 25496.3840 4221760.5030
34 28874.9170 4221756.7649
39 27284.2609 4220235.3943
133 26812.2676 4221811.7328] ;
olc=[0.00000
70.76351
105.09584
0.00000
42.40162
89.83755
0.00000
50.28006
98.61182
0
36.51504
106.45528
0
80.90273
160.68555
276.73136] ;
dn=[27;27;27;32;32;32;34;34;34;39;39;39;133;133;133;133] ;
bn=[34;133;32;27;133;39;39;133;27;32;133;34;27;34;39;32] ;
dog_olc=[dn bn olc] ;
N = size(dog_olc,1) ;
R=6373394 ;
Dog_olc = zeros(N,1) ;
for i = 1:N
dn=dog_olc(i,1); bn=dog_olc(i,2); olc=dog_olc(i,3);
xi=K(find(K(:,1)==dn,2)); yi=K(find(K(:,1)==dn,3));
xk=K(find(K(:,1)==bn,2)); yk=K(find(K(:,1)==bn,3));
dt=((-200/pi)/(6*R^2))*((xk-xi)*(2*yi+yk));
Dog_olc(i)=dog_olc(i,3)+dt ;
end
Dog_olc

4 个评论
显示 2更早的评论 隐藏 2更早的评论

Muhendisleksi
Muhendisleksi 2017-3-30
dt=[0.339375914
0.319298746
0.324864452
-0.317469398
-0.010411849
0.311830257
0.337917926
-0.01214174
-0.345692049
-0.318952756
-0.335085174
-0.331596985
-0.317306091
0.011845578
0.333141557
0.010587941]
Dog_olc=[0.000034
70.763542
105.095872
-0.000032
42.401619
89.837581
0.000034
50.280059
98.611785
-0.000032
36.515006
106.455247
-0.000032
80.902731
160.685583
276.731361]
The result should be like this. I think there is a slight mistake.
KSSV
KSSV 2017-3-30
The result matches with fourth column of your attached image..
Muhendisleksi
Muhendisleksi 2017-3-30
编辑：Muhendisleksi 2017-3-30
dt(1,1)= ((-200/pi*10000)/(6*R^2))*((4221756.7650-4223312.4563)*(2*27320.5963+28874.917))
Dog_olc(1,1)=0.339375914157676
This should be the process.
Muhendisleksi
Muhendisleksi 2017-3-30
I saw my mistake! thank you very, very much.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by