24小时热门版块排行榜    

查看: 362  |  回复: 0

xiayou89

银虫 (正式写手)

[求助] 解方程老是出错 呜呜呜

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)
就是方程逐个迭代  求大神知道额  错了一上午了
回复此楼
我们都要好好的
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiayou89 的主题更新
信息提示
请填处理意见