| 查看: 401 | 回复: 2 | ||||
tianpengusst铁虫 (小有名气)
|
[求助]
一个简单的程序 不知道错哪了 求大神点拨！！！ 已有1人参与
|
|
i=1;cd=0;xc=0.25; a=[0 0.07126 0.02046 0.1128 0.01073 -0.05203]; for s=1:1:100 while i<=181 fai(i)=a(1).*(1-cos(cd))+a(2).*sin(cd)+a(3).*(1-cos(cd)).^2+a(4).*(sin(cd)).^2+... a(5).*(1-cos(cd)).^3+a(6).*(sin(cd)).^3; i=i+1; cd=cd+pi/90; end cd=0;ii=1; while ii<=181 r=xc*exp(fai(ii)); x(ii)=(r+(xc^2)/r)*cos(cd)+2*xc; y(ii)=(r-(xc^2)/r)*sin(cd); cd=cd+pi/90; ii=ii+1; end c=[x' y']; a(2)=a(2)+0.1; end fai的值为什么在迭代过程中不变？应该变化的啊？就一个for循环。。。 |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
两块石头
已经有8人回复
-
初秋的晨风
已经有6人回复
-
学科评审组评审是指会评吗？
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有14人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有7人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68173.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
jjdg: 金币+2, 感谢参与 2017-05-02 01:37:08
感谢参与，应助指数 +1
jjdg: 金币+2, 感谢参与 2017-05-02 01:37:08
|
结果： fai = Columns 1 through 15 0 0.0026 0.0055 0.0086 0.0120 0.0155 0.0192 0.0231 0.0272 0.0313 0.0356 0.0399 0.0443 0.0488 0.0532 Columns 16 through 30 0.0577 0.0622 0.0667 0.0711 0.0755 0.0799 0.0841 0.0883 0.0925 0.0965 0.1005 0.1044 0.1082 0.1119 0.1155 Columns 31 through 45 0.1190 0.1224 0.1257 0.1290 0.1322 0.1353 0.1384 0.1413 0.1442 0.1471 0.1499 0.1527 0.1554 0.1580 0.1606 Columns 46 through 60 0.1632 0.1658 0.1682 0.1707 0.1731 0.1754 0.1777 0.1799 0.1821 0.1842 0.1862 0.1881 0.1899 0.1917 0.1933 Columns 61 through 75 0.1948 0.1961 0.1974 0.1985 0.1994 0.2002 0.2008 0.2012 0.2015 0.2016 0.2015 0.2012 0.2008 0.2001 0.1993 Columns 76 through 90 0.1983 0.1971 0.1958 0.1942 0.1926 0.1908 0.1888 0.1868 0.1846 0.1823 0.1800 0.1776 0.1752 0.1727 0.1702 Columns 91 through 105 0.1677 0.1652 0.1628 0.1604 0.1580 0.1558 0.1536 0.1516 0.1497 0.1479 0.1462 0.1447 0.1433 0.1420 0.1410 Columns 106 through 120 0.1400 0.1393 0.1386 0.1381 0.1378 0.1375 0.1374 0.1374 0.1375 0.1376 0.1378 0.1380 0.1383 0.1385 0.1387 Columns 121 through 135 0.1389 0.1391 0.1391 0.1391 0.1389 0.1386 0.1381 0.1375 0.1367 0.1357 0.1344 0.1330 0.1313 0.1294 0.1272 Columns 136 through 150 0.1248 0.1221 0.1192 0.1160 0.1126 0.1089 0.1051 0.1010 0.0968 0.0923 0.0877 0.0830 0.0781 0.0732 0.0682 Columns 151 through 165 0.0631 0.0581 0.0530 0.0480 0.0430 0.0381 0.0333 0.0287 0.0242 0.0199 0.0159 0.0120 0.0085 0.0052 0.0022 Columns 166 through 180 -0.0005 -0.0029 -0.0050 -0.0067 -0.0080 -0.0090 -0.0097 -0.0100 -0.0099 -0.0095 -0.0087 -0.0076 -0.0062 -0.0044 -0.0023 Column 181 -0.0000 为什么会说fai没有变化。 此外，代码中 for s=1:1:100 这一句有什么意义，加这一句会把中间两个while循环重复做上100遍，而100遍的结果都一样。 |
2楼2017-04-28 16:07:48
tianpengusst
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 131.8
- 散金: 81
- 帖子: 109
- 在线: 28.8小时
- 虫号: 4187746
- 注册: 2015-10-31
- 专业: 流体力学
3楼2017-05-01 16:41:49