24小时热门版块排行榜    

查看: 2842  |  回复: 15

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Qj98d49dHFb9 2026-08-25 4/200 2026-08-26 03:08 by cNXvBfCpiZOM
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +17 木水思豆 2026-08-25 20/1000 2026-08-26 01:26 by chengyan1220
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见