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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 求解一个计算结果里有符号的数值定积分
391/1返回列表
查看: 1913  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭