| 查看: 428 | 回复: 0 | ||
[求助]
解方程老是出错 呜呜呜
|
|
clear rt=0.05; rr=0.10; alpha=0.08722; gama=0.08722; theta=0.010; d0=0.30; n=1.4; c=300000; Q=0.00005; tp=5*10^(-9); k=0.5; beta=0.0026; %z=0.5:0.001:30; t=0:0.5:60; z=c.*t.^2.*0.5/(n.*(1+cos(theta))); r=z.*cos(theta); d=d0-r.*tan(theta); r1=(rt+z.*tan(alpha))*cos(theta); r2=rr+r.*tan(gama); w1=2.*acos((d.^2+r1.^2-r2.^2)./(2.*r1.*d)); w2=2.*acos((d.^2+r2.^2-r1.^2)./(2.*r2.*d)); y=(r1.^2.*(w1-sin(w1))+r2.^2.*(w2-sin(w2)))./(2.*pi*r1.^2); B=n./c.*pi.*rr.^2.*beta.*Q./tp.*((1+cos(theta))./(cos(theta)).^2); p=B.*int('y.*e^(-k.*c.*t./n).*(1/t.^2)',-tp,t); plot(t,p) 就是方程逐个迭代 求大神知道额 错了一上午了 |
回复此楼» 猜你喜欢
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有134人回复
-
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有107人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有62人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
