24小时热门版块排行榜    

查看: 1291  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +9 yangfz 2026-03-17 9/450 2026-03-18 12:38 by 尽舜尧1
[考研] 0703化学调剂 +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 341求调剂 +5 捣蛋猪猪 2026-03-11 7/350 2026-03-17 19:09 by 捣蛋猪猪
[考研] 求调剂,总分315,考的生物医药,一志愿湖南师范大学。调剂到任何专业都可以 +4 小丁想进步 2026-03-11 5/250 2026-03-17 16:05 by 外星文明
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
信息提示
请填处理意见