| 查看: 2640 | 回复: 0 | ||
[求助]
如何迭代计算编写下面程序(matlab), 求大神帮助.....
|
|
已知条件: Dj=4.14; d=1; k=0.235; t=1/(2*k); angle=46.5*pi/180; 下面思路是：将theta进行n等分，分别迭代计算： Dj=4.14; d=1; k=0.235; t=1/(2*k); angle=46.5*pi/180; L=d*tan(angle) a=angle; b=pi/2; n=10; h=(b-a)/n; i=1:1:n; theta(i)=angle+i*h; %根据下面两个式子分别求得x(0)...x(i-1)的值 tan(theta(i-1))*sqrt((1-(2*z-1)/z^2)*(1-1/((2*z-1)^t)^2))=1/((2*z-1)^t)+sqrt((2*z-1)/z^2) cos(x(i-1))/sin(psaiI-x(j))=sqrt(z^2/(2*z-1)) 其中迭代公式为; 当n=1时： L1=(tan(theta(1))*d-L)/(1+tan(x(0))*tan(theta(1)); %根据上面的式子可以预先求出x(0)、x(1)、x(2)、...、x(i-1)的值 h1=L1*tan(x(0)); H1=h1; 当n=2时; L2=(tan(2)*(d-L1*tan(x(0)))-(L+L1))/(1+tan(x(1))*tan(theta(2))); h2=L2*tan(x(2)); H2=h1+h2; ..... 当n=i时： L(i)=tan(i)*(d-(L1*tan(x(0))+L2*tan(x(1)+...+L(i-1)*tan(x(i-2))))-(L+L1+L2+...+L(i-1))/(1+tan(x(i-1))*tan(theta(i))) h(i)=L(i)*tan(x(i-1)); H(i)=h1+h2+...+h(i); 最终目的：求得随着theta的变化，H(i)值的变化，画出曲线。 |
回复此楼» 猜你喜欢
散金币祈福
已经有187人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有135人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有115人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有45人回复