应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlan求解含有三角函数的隐函数问题？
162/2返回列表上一页12
查看: 2557  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ju5200

木虫 (正式写手)
送红花一朵
自己琢磨了一会 这样应该可以了 （包括图）
clc
clear
format long
syms x2 x3 T
a=2.856;
X3=[];M=[];x2=[];
x1=pi/20:pi/20:pi/2
for i=1:length(x1)
    x00=x1(i);
    m=sqrt(1+(a+1)^2/(a*tan(x00))^2);
    M(i)=m;
    x20=atan(tan(x00)/(1+a*(1+tan(x00))*tan(x00)));
    %atan(x)范围-pi/2~pi/2，满足要求
    x2(i)=x20;
    [yy,fval,exitflag]=fmincon(@(x3)tan(x20)-(m^2*sin(x20+x3)*sin(x20+x3)-1)*cot(x20+x3)/(m^2*((a+1)/2-sin(x20+x3)*sin(x20+x3))-1),0,[],[],[],[],0,pi/2);
    X3(i)=yy;
end
M
x1
x2
X3
plot(x1,X3);
xlabel('x1')
ylabel('x3')
赞一下 回复此楼
11楼2014-01-07 20:51:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
引用回帖:
10楼: Originally posted by ju5200 at 2014-01-07 19:43:58
可否告知具体实现方法...

用1stOpt，代码如下：
CODE:
Constant a=2.856;
Parameter x2=[0,pi/2], x3=[0,pi/2];
LoopConstant x1=[pi/100:pi/40:pi/2];
PlotLoopData x1[x],x3;
Function 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);
赞一下 回复此楼
12楼2014-01-07 21:01:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ju5200

木虫 (正式写手)
至于为何答案与版主有出入 还请自己斟酌
赞一下 回复此楼
13楼2014-01-07 21:01:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ju5200

木虫 (正式写手)
引用回帖:
12楼: Originally posted by dingd at 2014-01-07 21:01:04
用1stOpt，代码如下：

Constant a=2.856;
Parameter x2=, x3=;
LoopConstant x1=;
PlotLoopData x1,x3;
Function m^2=1+(a+1)^2/(a*tan(x1))^2;
         tan(x2)=tan(x1)/(1+a*(1+tan(x1))*tan(x1));
   ...

那我这段代码有什么问题吗 结果稍有出入
赞一下 回复此楼
14楼2014-01-07 21:02:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ju5200

木虫 (正式写手)
M =
  Columns 1 through 3
   8.582903223762685   4.273938448777693   2.832213806724004
  Columns 4 through 6
   2.110286703931021   1.680142306733627   1.400796860655276
  Columns 7 through 9
   1.213774556514105   1.091991927151632   1.022608453435410
  Column 10
   1.000000000000000
x1 =
  Columns 1 through 3
   0.157079632679490   0.314159265358979   0.471238898038469
  Columns 4 through 6
   0.628318530717959   0.785398163397448   0.942477796076938
  Columns 7 through 9
   1.099557428756428   1.256637061435917   1.413716694115407
  Column 10
   1.570796326794897
x2 =
  Columns 1 through 3
   0.103555680176685   0.144718612570458   0.158063011752686
  Columns 4 through 6
   0.157235527421285   0.147898983200563   0.132316327109166
  Columns 7 through 9
   0.111015187136047   0.083343210741282   0.047478224447335
  Column 10
   0.000000000000000
X3 =
  Columns 1 through 3
   0.868368679046074   0.850111857513349   0.873011826888741
  Columns 4 through 6
   0.921909356113148   0.991041512984862   1.077308302925381
  Columns 7 through 9
   1.181128432161135   1.306571217046848   1.570796326794897
  Column 10
   1.570796326794897
>>
赞一下 回复此楼
15楼2014-01-07 21:02:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ju5200

木虫 (正式写手)
引用回帖:
14楼: Originally posted by ju5200 at 2014-01-07 21:02:16
那我这段代码有什么问题吗 结果稍有出入...

matlan求解含有三角函数的隐函数问题？
赞一下 回复此楼
16楼2014-01-07 21:04:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 815292578 的主题更新
162/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕326求调剂 +4 墨煜姒莘 2026-03-15 4/200 2026-03-15 11:02 by dyw
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[基金申请] 现在如何回避去年的某一个专家，不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 295复试调剂 +5 简木ChuFront 2026-03-09 5/250 2026-03-14 01:29 by JourneyLucky
[考研] 一志愿郑大070303，338分，求调剂 +4 dadawaf 2026-03-10 5/250 2026-03-14 01:20 by lsw010101
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 0805，333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 282分材料专业求调剂院校 +18 枫桥ZL 2026-03-09 25/1250 2026-03-13 10:47 by 白夜悠长
[考研] 一志愿华中师范071000，325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
[硕博家园] 木虫好像不热闹了，是不是？ +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
信息提示
关闭
请填处理意见
关闭