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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 求解一个计算结果里有符号的数值定积分
391/1返回列表
查看: 1832  |  回复: 38
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

雪径踏青

铁杆木虫 (著名写手)

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

如题，具体方程式如下：
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不可积，用数值法也没有成功
请各位赐教，不胜感激，谢谢！
回复此楼

» 猜你喜欢

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

查看全部散金贴
1楼 2017-11-30 19:42:18
已阅   回复此楼   关注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客户端
zhaoyan19783楼
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): 谢谢参与
敏敏的爸13楼
2017-11-30 22:22   回复  
雪径踏青(金币+1): 谢谢参与
发自小木虫Android客户端
纹哉il14楼
2017-11-30 22:34   回复  
雪径踏青(金币+1): 谢谢参与
敏敏的妈15楼
2017-11-30 22:38   回复  
雪径踏青(金币+1): 谢谢参与
发自小木虫Android客户端
dongmings16楼
2017-11-30 22:59   回复  
雪径踏青(金币+1): 谢谢参与
发自小木虫IOS客户端
大海的尽头17楼
2017-11-30 23:14   回复  
雪径踏青(金币+1): 谢谢参与
bkkkmze8919楼
2017-11-30 23:26   回复  
雪径踏青(金币+1): 谢谢参与
Emily爱生活21楼
2017-11-30 23:32   回复  
雪径踏青(金币+1): 谢谢参与
淡然云霞22楼
2017-12-01 00:19   回复  
雪径踏青(金币+1): 谢谢参与
心静茶凉23楼
2017-12-01 06:35   回复  
雪径踏青(金币+1): 谢谢参与
hcmeux75224楼
2017-12-01 07:35   回复  
雪径踏青(金币+1): 谢谢参与
丢了的欣欣25楼
2017-12-01 08:03   回复  
雪径踏青(金币+1): 谢谢参与
zyburga83129楼
2017-12-01 08:38   回复  
雪径踏青(金币+1): 谢谢参与
蝶舞水岸30楼
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): 谢谢参与
矿泉水和鱼35楼
2017-12-01 08:49   回复  
雪径踏青(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 雪径踏青 的主题更新
391/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +45 医学老男孩 2026-08-13 99/4950 2026-08-18 12:56 by dotnetchen9056
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 KZyU1THERRVo 2026-08-18 3/150 2026-08-18 12:52 by vHWVV9D5Deb5
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 5/250 2026-08-18 11:25 by 無關想念
[基金申请] 时间戳变了，能看出什么问题？ +5 基诺咪客 2026-08-17 6/300 2026-08-18 11:06 by hunter无悔
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
关闭
请填处理意见
关闭