24小时热门版块排行榜    

查看: 1984  |  回复: 38

雪径踏青

铁杆木虫 (著名写手)


[交流] 求解一个计算结果里有符号的数值定积分

如题,具体方程式如下:
syms Yt0 a
Z= -((56384229794665*((13*Yt0*exp((3*a)/20))/150000 - 1)^2*(1/20 - (13*Yt0*exp((3*a)/20))/3000000)^(27/100))/476053504 + (56384229794665*((13*Yt0*exp((3*a)/20))/150000 + 1)^2*((13*Yt0*exp((3*a)/20))/3000000 + 1/20)^(27/100))/476053504 - (146598997466129*Yt0*exp((3*a)/20)*((13*Yt0*exp((3*a)/20))/150000 - 1)*(1/20 - (13*Yt0*exp((3*a)/20))/3000000)^(27/100))/7990149120000 - (146598997466129*Yt0*exp((3*a)/20)*((13*Yt0*exp((3*a)/20))/150000 + 1)*((13*Yt0*exp((3*a)/20))/3000000 + 1/20)^(27/100))/7990149120000 + 1591911009283099/140737488355328)/(40000000*((13*Yt0*exp((3*a)/20))/3000000 - 1)^3)
A=int(Z,a,0,3.14/3)

其实我是想把Z方程中的a当作被积函数的未知数,将Yt0当作积分结果中的符号变量,也就是说我想把Z方程以a为自变量积分,积分上下限分别为0与3.14/6,最终的结果里面含有Yt0这个符号。
然而我用int不可积,用数值法也没有成功
请各位赐教,不胜感激,谢谢!
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

jackplay

金虫 (小有名气)


★
雪径踏青(金币+1): 谢谢参与
应该先简化下Z的表达式吧
12楼2017-11-30 22:12:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小橘子呆呆

银虫 (正式写手)


★
雪径踏青(金币+1): 谢谢参与
::::
18楼2017-11-30 23:21:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)


★
雪径踏青(金币+1): 谢谢参与
数值法需要确定参数具体的数值,可用匿名函数进行参数传递 z=quadl(@Yt0-((5638…^3),a,0,3.14/3)如果需要具体表达式可用maple计算,结果很复杂写出来有意义吗?

发自小木虫IOS客户端
20楼2017-11-30 23:30:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪径踏青

铁杆木虫 (著名写手)


引用回帖:
20楼: Originally posted by 独孤神宇 at 2017-11-30 23:30:48
数值法需要确定参数具体的数值,可用匿名函数进行参数传递 z=quadl(@Yt0-((5638…^3),a,0,3.14/3)如果需要具体表达式可用maple计算,结果很复杂写出来有意义吗?
...

你好,很感谢你的回复!有意义,我是想最终绘出Z-Yt0的曲线图,然后看这条曲线的规律。这个式子是之前其他的方程联立用MATLAB算出的结果,结果确实看上去很不规则,然而原始的算式比较正常了。
你刚刚说的这个方法,可以赐教得详细点吗?
26楼2017-12-01 08:07:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪径踏青

铁杆木虫 (著名写手)


引用回帖:
12楼: Originally posted by jackplay at 2017-11-30 22:12:20
应该先简化下Z的表达式吧

你好,谢谢你的回复。我觉得表达式没必要简化吧?毕竟最终结果是要算出来的
27楼2017-12-01 08:07:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪径踏青

铁杆木虫 (著名写手)


引用回帖:
20楼: Originally posted by 独孤神宇 at 2017-11-30 23:30:48
数值法需要确定参数具体的数值,可用匿名函数进行参数传递 z=quadl(@Yt0-((5638…^3),a,0,3.14/3)如果需要具体表达式可用maple计算,结果很复杂写出来有意义吗?
...

刚刚说错了,应该是得出A-Yt0的曲线图,谢谢
28楼2017-12-01 08:08:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)


引用回帖:
26楼: Originally posted by 雪径踏青 at 2017-12-01 08:07:05
你好,很感谢你的回复!有意义,我是想最终绘出Z-Yt0的曲线图,然后看这条曲线的规律。这个式子是之前其他的方程联立用MATLAB算出的结果,结果确实看上去很不规则,然而原始的算式比较正常了。
你刚刚说的这个方法 ...

Yto范围呢?

发自小木虫IOS客户端
36楼2017-12-02 00:02:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪径踏青

铁杆木虫 (著名写手)


引用回帖:
36楼: Originally posted by 独孤神宇 at 2017-12-02 00:02:30
Yto范围呢?
...

0,360

发自小木虫Android客户端
37楼2017-12-04 00:25:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
可用 quadl+arrayfun 求解,自己可以找找例子学习一下,以后还会用得到的

发自小木虫IOS客户端
38楼2017-12-04 19:21:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪径踏青

铁杆木虫 (著名写手)


引用回帖:
38楼: Originally posted by 独孤神宇 at 2017-12-04 19:21:27
可用 quadl+arrayfun 求解,自己可以找找例子学习一下,以后还会用得到的

谢谢你朋友,我明天试试!!!

发自小木虫Android客户端
39楼2017-12-04 21:52:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
nicenicey2楼
2017-11-30 19:47   回复  
雪径踏青(金币+1): 谢谢参与
发自小木虫Android客户端
2017-11-30 20:00   回复  
雪径踏青(金币+1): 谢谢参与
发自小木虫Android客户端
tzynew4楼
2017-11-30 20:05   回复  
雪径踏青(金币+1): 谢谢参与
毓肯315楼
2017-11-30 20:11   回复  
雪径踏青(金币+1): 谢谢参与
knwxwm6906楼
2017-11-30 20:44   回复  
雪径踏青(金币+1): 谢谢参与
顶
rvwmea3137楼
2017-11-30 20:57   回复  
雪径踏青(金币+1): 谢谢参与
afyrux388楼
2017-11-30 21:14   回复  
雪径踏青(金币+1): 谢谢参与
谆蚀mr9楼
2017-11-30 21:28   回复  
雪径踏青(金币+1): 谢谢参与
syhorchid10楼
2017-11-30 21:30   回复  
雪径踏青(金币+1): 谢谢参与
pkftp55111楼
2017-11-30 21:57   回复  
雪径踏青(金币+1): 谢谢参与
2017-11-30 22:22   回复  
雪径踏青(金币+1): 谢谢参与
发自小木虫Android客户端
纹哉il14楼
2017-11-30 22:34   回复  
雪径踏青(金币+1): 谢谢参与
2017-11-30 22:38   回复  
雪径踏青(金币+1): 谢谢参与
发自小木虫Android客户端
dongmings16楼
2017-11-30 22:59   回复  
雪径踏青(金币+1): 谢谢参与
发自小木虫IOS客户端
2017-11-30 23:14   回复  
雪径踏青(金币+1): 谢谢参与
bkkkmze8919楼
2017-11-30 23:26   回复  
雪径踏青(金币+1): 谢谢参与
2017-11-30 23:32   回复  
雪径踏青(金币+1): 谢谢参与
2017-12-01 00:19   回复  
雪径踏青(金币+1): 谢谢参与
2017-12-01 06:35   回复  
雪径踏青(金币+1): 谢谢参与
hcmeux75224楼
2017-12-01 07:35   回复  
雪径踏青(金币+1): 谢谢参与
2017-12-01 08:03   回复  
雪径踏青(金币+1): 谢谢参与
zyburga83129楼
2017-12-01 08:38   回复  
雪径踏青(金币+1): 谢谢参与
2017-12-01 08:41   回复  
雪径踏青(金币+1): 谢谢参与
BACKTOHO31楼
2017-12-01 08:44   回复  
雪径踏青(金币+1): 谢谢参与
nono200932楼
2017-12-01 08:47   回复  
雪径踏青(金币+1): 谢谢参与
`
nono200933楼
2017-12-01 08:47   回复  
`
sxccrbz8534楼
2017-12-01 08:47   回复  
雪径踏青(金币+1): 谢谢参与
2017-12-01 08:49   回复  
雪径踏青(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 雪径踏青 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:27 by xzcv0Tehd8lQ
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:27 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:01 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 19:50 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 16:58 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:45 by BCdvBAEpt3zm
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 14:50 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 14:50 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 10:43 by BCdvBAEpt3zm
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
信息提示
请填处理意见