[求救] 叠代循环写不出来

楼主: ccam (C.C.C)   2017-12-03 16:09:03
如题,叠代的算法卡住了QQ
0~90度之内,每一度对应到一个 r 值
已知 r(1)=0 ,用循环叠代方式从前一个 r 算出下一个 r
假定 r(n+1)=sqrt(A+B)
其中 A 是由 r(n) 算出
B 则是依据 r(n) 的大小来决定其中一个变量
B 的算式中有一个 C 值
C 是 1X50 的数列(单纯的正自然数,无规则)
然后现在有一串跟 r(n) 比大小的数列 R 是 1X51 (由0开始的等差数列)
如果 r(n) 的大小落在 R 的第一数值~第二数值的区间内,B 的计算就使用 C(1)
再求出 r(n+1)
如果 r(n) 的大小落在 R 的第二数值~第三数值的区间内,B 的计算就使用 C(2)
再求出 r(n+1)
以次类推到最后一个区间,求出 r(91),总共91个 r 值
注:R 的一个区间内可能同时对应到 好几个角度
这个问题已经卡住两礼拜了,还请各位指点迷津
目前自己只写出一个很怪的东西
for i=0:1:90;
for k=1:50;
n = n+1;
if r(n) >= R(k) && r(n) <= R(k+1)
r(n+1) = sqrt(((sind(i)^2)/(C(k)))+r(n)^2);
else
r(n) > R(k+1);
break
end
end
end
作者: LiamIssac (Madchester)   2017-12-03 16:32:00
有跑过吗? 问题是?你的n要换个位置 而且没有初始值然后else里面有没有运算 试着在每个运算后面加pause让程式暂停一下 看看跟你用手算的也不一样然后把if的条件 也独立出来跑看看 确定程式会不会进入那个condition去算你要的东西
作者: tml (流刑人形)   2017-12-11 06:41:00
for i=0:90r(i+1)=sqrt(sind(i)^2/C(find(r(i)>=R,1))+r(i)^2);end

Links booklink

Contact Us: admin [ a t ] ucptt.com