| 查看: 322 | 回复: 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循环。。。 |
回复此楼» 猜你喜欢
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
假如你的研究生提出不合理要求
已经有5人回复
-
所感
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有28人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68123.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.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