姜启源《数学模型》——血管分支

案例选自《数学模型》血管分支,水平有限,如有错误敬请谅解。

syms r1 r theta alpha k q H L b
fun=@(r1,r,theta)(k*q^2/r^4+b*r^alpha)*(L-H/tan(theta))+(0.25*k*q^2/(r1)^4+b*(r1)^alpha)*2*H/sin(theta)
res1=diff(fun,r1)
res2=diff(fun,r)
res3=diff(fun,theta)
a=1:0.01:2;
y=power(2,(a-4)./(a+4));
miny=min(y);
maxy=max(y);
rr=power(4,1.0./(a+4));
minrr=min(rr);
maxrr=max(rr);
theta1=acos(miny)*180/pi;
theta2=acos(maxy)*180/pi;
fprintf('theta的范围大概是%f~%f',theta2,theta1)
fprintf('r/r1的范围大概是%f~%f',minrr,maxrr)


运行结果如下: (实时编辑器内运行结果)

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值