如果您使用了这种方法，您的代码将更具可读性(因此也更容易调试)
cosd (x)
和
信德(x)
而不是
因为(π* x / 180)
和
罪(π* x / 180)。
用fsolve求解非线性方程
3次查看(最近30天)
显示旧的评论
函数F = root2d(φ)
F = 0 (2, 1);
F(1) =(π* sin((π*(φ(1)+ theta_n)) / 180) * cos((π*φ(2))/ 180)* (cos((π* theta_n) / 180) +棕褐色((π*我)/ 180)* tan((π* theta_i) / 180))) / (180 * sin((π*φ(1))/ 180)* (cos((π*(φ(1)+ theta_n)) / 180) * cos((φ(2)*π)/ 180)+罪((φ(2)*π)/ 180)* tan ((theta_i *π)/ 180))^ 2)——(π* cos((π*φ(1))/ 180)* (cos((π* theta_n) / 180) +谭((π*我)/ 180)* tan((π* theta_i) / 180))) /(180 *罪((φ(1)*π)/ 180)^ 2 * (cos((π*(φ(1)+ theta_n)) / 180) * cos((π*φ(2))/ 180)+罪((π*φ(2))/ 180)* tan((π* theta_i) / 180)));
F (2) = ((cos((π* theta_n) / 180) +棕褐色((π*我)/ 180)* tan((π* theta_i) / 180)) *((π* cos((π*(φ(1)+ theta_n)) / 180) *罪((π*φ(2))/ 180))/ 180(π* cos((π*φ(2))/ 180)* tan((π* theta_i) / 180)) / 180)) /(罪((π*φ(1))/ 180)* (cos((π*(φ(1)+ theta_n)) / 180) * cos((φ(2)*π)/ 180)+罪((φ(2)*π)/ 180)* tan ((theta_i *π)/ 180))^ 2);
结束
beta_a = 10;
Alpha_n = 10;
η= 18;
我= 18;
theta_i = 10;
= 10;
有趣= @root2d;
phi0 = (0,0);
φ= fsolve(有趣,phi0);
phi_n =φ(1);
phi_i =φ(2);
我得到错误:
未被认可的函数或变量'theta_n'。
错误在root2d(第4行)
F(1) =(π* sin((π*(φ(1)+ theta_n)) / 180) * cos((π*φ(2))/ 180)* (cos((π* theta_n) / 180) +
谭((π*我)/ 180)* tan((π* theta_i) / 180))) /(180 *罪((π*φ(1))/ 180)* (cos((π*(φ(1)+
theta_n)) / 180) * cos((φ(2)*π)/ 180)+罪((φ(2)*π)/ 180)* tan ((theta_i *π)/ 180))^ 2)
(π* cos((π*φ(1))/ 180)* (cos((π* theta_n) / 180) +
谭((π*我)/ 180)* tan((π* theta_i) / 180))) /(180 *罪((φ(1)*π)/ 180)^ 2 * (cos((π*(φ(1)+
theta_n)) / 180) * cos((π*φ(2))/ 180)+罪((π*φ(2))/ 180)* tan((π* theta_i) / 180)));
错误在fsolve(第258行)
熔化炉=函数宏指令(funfcn {3}, x,变长度输入宗量{:});
错误在最近(第12行)
φ= fsolve(有趣,phi0);
引起的由:
失败在初始目标函数评估中。FSOLVE无法继续。
接受的答案
马特·J
2021年2月2日
编辑:马特·J
2021年2月2日
函数solve_it ()
beta_a = 10;
Alpha_n = 10;
η= 18;
我= 18;
theta_i = 10;
= 10;
有趣= @root2d;
phi0 = (10,10);
φ= fsolve(有趣,phi0);
phi_nφ= (1),
phi_iφ= (2),
函数F = root2d(φ)
F = 0 (2, 1);
F(1) =(π* sin((π*(φ(1)+ theta_n)) / 180) * cos((π*φ(2))/ 180)* (cos((π* theta_n) / 180) +棕褐色((π*我)/ 180)* tan((π* theta_i) / 180))) / (180 * sin((π*φ(1))/ 180)* (cos((π*(φ(1)+ theta_n)) / 180) * cos((φ(2)*π)/ 180)+罪((φ(2)*π)/ 180)* tan ((theta_i *π)/ 180))^ 2)——(π* cos((π*φ(1))/ 180)* (cos((π* theta_n) / 180) +谭((π*我)/ 180)* tan((π* theta_i) / 180))) /(180 *罪((φ(1)*π)/ 180)^ 2 * (cos((π*(φ(1)+ theta_n)) / 180) * cos((π*φ(2))/ 180)+罪((π*φ(2))/ 180)* tan((π* theta_i) / 180)));
F (2) = ((cos((π* theta_n) / 180) +棕褐色((π*我)/ 180)* tan((π* theta_i) / 180)) *((π* cos((π*(φ(1)+ theta_n)) / 180) *罪((π*φ(2))/ 180))/ 180(π* cos((π*φ(2))/ 180)* tan((π* theta_i) / 180)) / 180)) /(罪((π*φ(1))/ 180)* (cos((π*(φ(1)+ theta_n)) / 180) * cos((φ(2)*π)/ 180)+罪((φ(2)*π)/ 180)* tan ((theta_i *π)/ 180))^ 2);
结束
结束
更多的答案(0)
你也可以从以下列表中选择一个网站: