24小时热门版块排行榜    

查看: 1610  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gaojinggigi

新虫 (正式写手)

[求助] 求助matlab写循环

各位大侠,我现在有以下程序,不知道怎么写循环,能不能帮我看一下?
a=0.8;
n=1.45;
f=100;
x0=a;y0=0;
gamma1=atan((x0+a)/(y0+f))
alpha1=gamma1/(n-1)
beta1=alpha1*n
l1=2*a*sin(alpha1)/cos(beta1)

x1=x0+(2*a+l1*cos(gamma1))
y1=y0+(l1*sin(gamma1))

gamma2=atan((x1+a)/(y1+f))
alpha2=gamma2/(n-1)
beta2=alpha2*n
l2=2*a*sin(alpha2)/cos(beta2)

x2=x1+(2*a+l2*cos(gamma2))
y2=y1+(l2*sin(gamma2))

gamma3=atan((x2+a)/(y2+f))
alpha3=gamma3/(n-1)
beta3=alpha3*n
l3=2*a*sin(alpha3)/cos(beta3)

x3=x2+(2*a+l3*cos(gamma3))
y3=y2+(l3*sin(gamma3))

在x<20,y<0.5时要一直算出xn,yn来,最后得出一组(x,y)的数据,请问怎么写循环语句啊?这么多变量我已经完全晕掉了,求赐教~
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

加油加油~
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaojinggigi

新虫 (正式写手)

引用回帖:
2楼: Originally posted by csgt0 at 2012-12-21 09:25:31
a=0.8;
n=1.45;
f=100;
i=1;
x(1)=a;
y(1)=0;

while x(end)<20&&y(end)<0.5
     gamma1=atan((x(end)+a)/(y(end)+f));
     alpha1=gamma1/(n-1);
     beta1=alpha1*n;
      ...

大侠,我后来发现每计算一次出来y>0.4,需要改变x(i)和y(i)和a的初值,x(i)和a的初值是下面程序的值,y(i)=0,再进行之前的循环,再求出y>0.4,改变x(i)和y(i)和a,一直计算到x>20结束。可是我加上if条件之后不知道怎么让它返回继续循环,还有x>20停止循环的命令怎么加?
我不会加循环,每次计算完了找到y>0.4的值再用公式求出a和x再代入,发现做了十几步还没做完~愁死我了~

a0=0.8;
a=-0.8;
n=1.411;
f=100;
x(1)=a;y(1)=0;

for ii=2:1:100
    gamma(ii)=atan((x(ii-1)-a)/(y(ii-1)+f));
    alpha(ii)=atan(sin(gamma(ii)/(n-cos(gamma(ii)))));
    beta(ii)=asin(sin(alpha(ii))*n);
    l(ii)=(a0-a)*sin(alpha(ii))/cos(beta(ii));

    x(ii)=x(ii-1)+(a0-a+l(ii)*sin(gamma(ii)))
    y(ii)=y(ii-1)+(l(ii)*cos(gamma(ii)))
  if y(ii)>0.4
    x(ii)=(0.4-y(ii-1))/tan(alpha(ii));
    a=x(ii)-(f+0.4)*tan(gamma(ii));
    y(ii)=0;
    ii=ii+1;
   end
  
end
加油加油~
3楼2012-12-21 15:59:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
gaojinggigi: 金币+5, ★★★很有帮助 2012-12-21 19:03:33
CODE:
a=0.8;
n=1.45;
f=100;
i=1;
x(1)=a;
y(1)=0;

while x(end)<20&&y(end)<0.5
     gamma1=atan((x(end)+a)/(y(end)+f));
     alpha1=gamma1/(n-1);
     beta1=alpha1*n;
     l1=2*a*sin(alpha1)/cos(beta1);     
     i=i+1;
     x(i)=x(end)+(2*a+l1*cos(gamma1));
     y(i)=y(end)+(l1*sin(gamma1));
end

showmethemoney
2楼2012-12-21 09:25:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

去掉ii=ii+1;
在最后的end前加入
if x(ii)>20
break;
end
showmethemoney
4楼2012-12-21 16:13:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaojinggigi

新虫 (正式写手)

引用回帖:
4楼: Originally posted by csgt0 at 2012-12-21 16:13:44
去掉ii=ii+1;
在最后的end前加入
if x(ii)>20
break;
end

这样算出来的y在第一次大于0.4后,后面都是0了,不对啊
加油加油~
5楼2012-12-21 16:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息,中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[教师之家] 两姐妹,一个辛辛苦苦30岁工作当大学教师,另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 7/350 2026-08-18 16:42 by LNP@mRNA
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见