| 查看: 2556 | 回复: 15 | ||
[求助]
matlan求解含有三角函数的隐函数问题？ 已有2人参与
|
||
|
已知条件： a=2.856; m^2=1+(a+1)^2/(a*tan(x1))^2; tan(x2)=tan(x1)/(1+a*(1+tan(x1))*tan(x1)); tan(x2)=(m^2*sin(x3+x2)*sin(x3+x2)-1)*cot(x3+x2)/(m^2*((a+1)/2-sin(x3+x2)*sin(x3+x2))-1); 求解：当x1的取值为x1=0:pi/20:pi/2时，角度x3的大小？ 并且画出图像，横坐标是x1；纵坐标是x3。 |
回复此楼» 本帖已获得的红花（最新10朵）
ju5200» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- fortran中如何调用matlab中求解特征值函数？
已经有9人回复
- 求助：Matlab中如何绘制三维隐函数图像
已经有20人回复
- MATLAB-pdepe函数使用问题
已经有5人回复
- Matlab画隐函数，ezplot
已经有3人回复
- 求助matlab的发fgoalattain函数
已经有4人回复
- 求matlab画一个函数的图
已经有6人回复
- 小弟新学matlab，想借用 fmincon求个函数的最小值，求助错误该法，感谢！
已经有4人回复
- matlab的函数调用问题
已经有3人回复
- 问个关于matlab里find函数的问题！
已经有4人回复
- matlab 求解隐式常微分方程时用ode15i 函数出现的问题求助
已经有3人回复
- 求助：能不能帮忙写一个用matlab求一个函数的傅里叶变换的程序
已经有16人回复
- matlab做函数图象的问题
已经有6人回复
- 【求助】求解隐函数
已经有6人回复
- 【求助】关于matlab中的函数代换
已经有10人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
|
x2和x3都在（0，pi/2)内: 循环常数 x1 x2 x3 m 0.0314159265358979 0.0287555833592285 0.0351971731972987 -42.9737897661545 0.109955742875643 0.0815900212903 0.111275584618701 -12.2702325504199 0.188495559215388 0.115188464315113 0.177382315237603 7.14797754409201 0.267035375555132 0.136273601407194 0.238672389800242 5.0355731493838 0.345575191894877 0.148996450783062 0.298087032110649 -3.88119660345004 0.424115008234622 0.155948935701275 0.357421640180446 -3.15300936364694 0.502654824574367 0.158777851292948 0.417836258978744 2.65168312922698 0.581194640914112 0.15854561247721 0.480127514682949 2.2857461469141 0.659734457253857 0.15594302634527 0.544877738661448 -2.0073997605745 0.738274273593601 0.151416611093451 0.612533446816156 -1.78930617394594 0.816814089933346 0.14524587471559 0.683439393079181 1.61477055675776 0.895353906273091 0.137590506742645 0.757841229168019 -1.47309255064489 0.973893722612836 0.128518804248998 0.835861942536153 1.35716829219462 1.05243353895258 0.118023779705002 0.91745182742723 1.26216627285968 1.13097335529233 0.10603055952273 1.00230778118668 -1.18475367370879 1.20951317163207 0.0923969359011175 1.08975826774206 -1.12262146467709 1.28805298797182 0.0769077582425526 1.17863687334423 -1.07417949665181 1.36659280431156 0.0592629423303635 1.26731483330163 -1.03835278050158 1.4451326206513 0.0390580710377663 1.35463513951086 1.0144406021796 1.52367243699105 0.015755537356128 1.44515570849173 -1.00202494300619 1.5707963267949 4.29923650476096E-31 1.5707963179239 -1 |
» 本帖已获得的红花（最新10朵）
7楼2014-01-07 17:55:21
ju5200
木虫 (正式写手)
- 应助: 78 (初中生)
- 金币: 1778.6
- 散金: 163
- 红花: 12
- 帖子: 511
- 在线: 65.2小时
- 虫号: 1734187
- 注册: 2012-04-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
感谢参与，应助指数 +1
|
楼上正解 clc clear format long syms x2 x3 a=2.856; x1=pi/20:pi/20:pi/2; m=sqrt(1+(a+1)^2./(a*tan(x1)).^2); %y=solve('tan(x2)-tan(x1).*tan(x1)=0',x2) %此处用点乘不对 x2=subs(solve('tan(x2)=tan(x1)/(1+a*(1+tan(x1))*tan(x1))',x2)); Y=subs(solve('tan(x2)=(m^2*sin(x3+x2)*sin(x3+x2)-1)*cot(x3+x2)/(m^2*((a+1)/2-sin(x3+x2)*sin(x3+x2))-1)',x3)); Y' ANSS=vpa(Y',5) 结果： ans = 2*conj(atan(z)) + 2*pi*conj(k) - 1865493290623271/18014398509481984 2*conj(atan(z)) + 2*pi*conj(k) - 1303509379291781/9007199254740992 2*conj(atan(z)) + 2*pi*conj(k) - 5694820166643655/36028797018963968 2*conj(atan(z)) + 2*pi*conj(k) - 5665006901631221/36028797018963968 2*conj(atan(z)) + 2*pi*conj(k) - 666077805630529/4503599627370496 2*conj(atan(z)) + 2*pi*conj(k) - 2383599045855489/18014398509481984 2*conj(atan(z)) + 2*pi*conj(k) - 7999487286693907/72057594037927936 2*conj(atan(z)) + 2*pi*conj(k) - 6005511245412801/72057594037927936 2*conj(atan(z)) + 2*pi*conj(k) - 6842333245735417/144115188075855872 2*conj(atan(z)) + 2*pi*conj(k) - 3478821848754559/162259276829213363391578010288128 ANSS = 2.0*conj(atan(z)) + 6.2832*conj(k) - 0.10356 2.0*conj(atan(z)) + 6.2832*conj(k) - 0.14472 2.0*conj(atan(z)) + 6.2832*conj(k) - 0.15806 2.0*conj(atan(z)) + 6.2832*conj(k) - 0.15724 2.0*conj(atan(z)) + 6.2832*conj(k) - 0.1479 2.0*conj(atan(z)) + 6.2832*conj(k) - 0.13232 2.0*conj(atan(z)) + 6.2832*conj(k) - 0.11102 2.0*conj(atan(z)) + 6.2832*conj(k) - 0.083343 2.0*conj(atan(z)) + 6.2832*conj(k) - 0.047478 2.0*conj(atan(z)) + 6.2832*conj(k) - 2.144*10^(-17) |
4楼2014-01-07 12:32:35
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2014-01-07 09:58:10
3楼2014-01-07 10:05:03
ju5200
木虫 (正式写手)
- 应助: 78 (初中生)
- 金币: 1778.6
- 散金: 163
- 红花: 12
- 帖子: 511
- 在线: 65.2小时
- 虫号: 1734187
- 注册: 2012-04-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2014-01-07 12:34:16
6楼2014-01-07 17:49:06
8楼2014-01-07 19:40:02
ju5200
木虫 (正式写手)
- 应助: 78 (初中生)
- 金币: 1778.6
- 散金: 163
- 红花: 12
- 帖子: 511
- 在线: 65.2小时
- 虫号: 1734187
- 注册: 2012-04-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2014-01-07 19:43:18
ju5200
木虫 (正式写手)
- 应助: 78 (初中生)
- 金币: 1778.6
- 散金: 163
- 红花: 12
- 帖子: 511
- 在线: 65.2小时
- 虫号: 1734187
- 注册: 2012-04-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2014-01-07 19:43:58
