24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

王风0304

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 55.8
散金: 20
帖子: 24
在线: 1.7小时
虫号: 21952660
注册: 2020-04-22
专业: 机械结构强度学

[求助] 小白请教大佬，我用Matlab求解下图中的k怎么求解不出？已有1人参与

>> syms x,k
t=(k*51*0.9/(sqrt(72.5^2-x^2)-72.5+51+20))
28585.35=1327767.26-int(pi*(sqrt((72.5-t)^2-x^2)-((72.5-t)-51-(20-t)))^2,-50,50)
vpa(k)

发自小木虫Android客户端

回复此楼

» 猜你喜欢

球磨粉体时遇到了大的问题，请指教！已经有9人回复
江汉大学解明教授课题组招博士研究生/博士后已经有3人回复

1楼 2020-04-29 22:11:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

应助: 387 (硕士)
金币: 17803.5
红花: 53
帖子: 2879
在线: 606.3小时
虫号: 1517335
注册: 2011-11-30
性别: GG
专业: 常微分方程与动力系统

【答案】应助回帖

for k = 0.90241673:1e-10:0.90241674
syms x
t=(k*51*0.9/(sqrt(72.5^2-x^2)-72.5+51+20));
eq=pi*(sqrt((72.5-t)^2-x^2)-((72.5-t)-51-(20-t)))^2;
func=@(x)eval(eq);
S=1327767.26-quad(func,-50,50);%double();
if abs( S-28585.35)<=1e-6
k
end
end