24小时热门版块排行榜    

查看: 1293  |  回复: 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的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

去掉ii=ii+1;
在最后的end前加入
if x(ii)>20
break;
end
showmethemoney
4楼2012-12-21 16:13:44
已阅   回复此楼   关注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的回帖

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工调剂 324分 +7 llllkkkhh 2026-03-18 7/350 2026-03-18 17:18 by 216的故事
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +7 步川酷紫123 2026-03-13 7/350 2026-03-18 17:12 by 尽舜尧1
[考研] 281求调剂(0805) +5 烟汐忆海 2026-03-16 13/650 2026-03-18 14:30 by stone_128
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 085600材料与化工求调剂 +5 绪幸与子 2026-03-17 5/250 2026-03-17 16:40 by laoshidan
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 化学调剂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
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
信息提示
请填处理意见