设为首页收藏本站联系我们
切换到宽版

声振论坛

 找回密码
 我要加入

QQ登录

只需一步，快速开始

声振论坛»论坛 › 算法编程 › Matlab讨论区 › 用matlab求解数值积分
返回列表 发新帖
查看: 2080|回复: 11

[编程技巧] 用matlab求解数值积分

[复制链接]
zhihaikou
发表于 2008-9-22 12:25 | 显示全部楼层 |阅读模式

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要 登录 才可以下载或查看，没有账号？我要加入

x
用matlab求解数值积分

函数：
function du=eqns2(x,u)

du=zeros(4,1);
du(1)=u(2); du(2)=u(3);
du(3)=u(4);
du(4)=-3/[u(2)*u(4)/u(1)-u(3)/u(1)^4+u(2)^2/u(1)^5-
... 1/((u(1)+1)*u(1)^3)];

命令窗口：

y0=[1 0 0 0]; [x,num_y2]=ode45('eqns2',[0,10],y0);



我想求num_y2(:,1)（也就是u(x)）在区间[0,10]的积分怎么办？ 怎么写函数？

用quad计算积分，需要给出函数名， I=quad('num_y2(:,1)',1,10) 这样写就出错了

本帖被以下淘专辑推荐:

回复
分享到：

使用道具 举报

sigma665
发表于 2008-9-22 12:39 | 显示全部楼层

回复 楼主 zhihaikou 的帖子

num_y2(:,1)是什么
回复 支持 反对

使用道具 举报

messenger
发表于 2008-9-22 12:42 | 显示全部楼层
用trapz即可
回复 支持 反对

使用道具 举报

zhihaikou
 楼主| 发表于 2008-9-22 12:43 | 显示全部楼层
num_y2(:,1)是什么
——————————————
num_y2(:,1)是求解的函数

num_y2(:,2)是求解函数的二阶导数

依此类推
回复 支持 反对

使用道具 举报

zhihaikou
 楼主| 发表于 2008-9-22 12:44 | 显示全部楼层

回复 板凳 messenger 的帖子

直接
trapz('num_y2(:,1)',1,10)

？？？
回复 支持 反对

使用道具 举报

zhihaikou
 楼主| 发表于 2008-9-22 12:45 | 显示全部楼层
>> trapz('num_y2(:,1)',1,10)
??? Error using ==> trapz
LENGTH(X) must equal the length of the DIM'th dimension of Y.
回复 支持 反对

使用道具 举报

sigma665
发表于 2008-9-22 13:43 | 显示全部楼层

回复 地板 zhihaikou 的帖子

我问的是里面的东西
数据？
把他贴出来
回复 支持 反对

使用道具 举报

messenger
发表于 2008-9-22 14:02 | 显示全部楼层
自己看看帮助trapz怎么用，很简单

最讨厌告诉命令以后，再问命令怎么用了


原帖由 zhihaikou 于 2008-9-22 12:45 发表

                               
登录/注册后可看大图

>> trapz('num_y2(:,1)',1,10)
??? Error using ==> trapz
LENGTH(X) must equal the length of the DIM'th dimension of Y.
回复 支持 反对

使用道具 举报

zhihaikou
 楼主| 发表于 2008-9-22 14:04 | 显示全部楼层
x =

            0
  3.4276e-026
  6.8552e-026
  1.0283e-025
  1.3710e-025
  3.0848e-025
  4.7986e-025
  6.5124e-025
  8.2262e-025
  1.6795e-024
  2.5364e-024
  3.3933e-024
  4.2502e-024
  8.5347e-024
  1.2819e-023
  1.7104e-023
  2.1388e-023
  4.2811e-023
  6.4233e-023
  8.5656e-023
  1.0708e-022
  2.1419e-022
  3.2130e-022
  4.2842e-022
  5.3553e-022
  1.0711e-021
  1.6067e-021
  2.1422e-021
  2.6778e-021
  5.3556e-021
  8.0334e-021
  1.0711e-020
  1.3389e-020
  2.6778e-020
  4.0167e-020
  5.3556e-020
  6.6945e-020
  1.3389e-019
  2.0084e-019
  2.6778e-019
  3.3473e-019
  6.6945e-019
  1.0042e-018
  1.3389e-018
  1.6736e-018
  3.3473e-018
  5.0209e-018
  6.6945e-018
  8.3682e-018
  1.6736e-017
  2.5105e-017
  3.3473e-017
  4.1841e-017
  8.3682e-017
  1.2552e-016
  1.6736e-016
  2.0920e-016
  4.1841e-016
  6.2761e-016
  8.3682e-016
  1.0460e-015
  2.0920e-015
  3.1381e-015
  4.1841e-015
  5.2301e-015
  1.0460e-014
  1.5690e-014
  2.0920e-014
  2.6151e-014
  5.2301e-014
  7.8452e-014
  1.0460e-013
  1.3075e-013
  2.6151e-013
  3.9226e-013
  5.2301e-013
  6.5376e-013
  1.3075e-012
  1.9613e-012
  2.6151e-012
  3.2688e-012
  6.5376e-012
  9.8064e-012
  1.3075e-011
  1.6344e-011
  3.2688e-011
  4.9032e-011
  6.5376e-011
  8.1720e-011
  1.6344e-010
  2.4516e-010
  3.2688e-010
  4.0860e-010
  6.5860e-010
  9.0860e-010
  1.1586e-009
  1.4086e-009
  1.6586e-009
  1.9086e-009
  2.1586e-009
  2.4086e-009
  2.6586e-009
  2.9086e-009
  3.1586e-009
  3.4086e-009
  3.6586e-009
  3.9086e-009
  4.1586e-009
  4.4086e-009
  4.6586e-009
  4.9086e-009
  5.1586e-009
  5.4086e-009
  5.6586e-009
  5.9086e-009
  6.1586e-009
  6.4086e-009
  6.6586e-009
  6.9086e-009
  7.1586e-009
  7.4086e-009
  7.6586e-009
  7.9086e-009
  8.1586e-009
  8.4086e-009
  8.6586e-009
  8.9086e-009
  9.1586e-009
  9.4086e-009
  9.5565e-009
  9.7043e-009
  9.8522e-009
  1.0000e-008


num_y2 =

  1.7182e-009            0            0            0
  1.7182e-009  9.8369e-057  8.6097e-031  5.0238e-005
  1.7182e-009  7.8695e-056  3.4439e-030  1.0048e-004
  1.7182e-009  2.6560e-055  7.7488e-030  1.5071e-004
  1.7182e-009  6.2956e-055  1.3776e-029  2.0095e-004
  1.7182e-009  7.1711e-054  6.9739e-029  4.5214e-004
  1.7182e-009  2.6993e-053  1.6875e-028  7.0333e-004
  1.7182e-009  6.7471e-053  3.1081e-028  9.5452e-004
  1.7182e-009  1.3599e-052  4.9592e-028  1.2057e-003
  1.7182e-009  1.1573e-051  2.0672e-027  2.4616e-003
  1.7182e-009  3.9862e-051  4.7147e-027  3.7176e-003
  1.7182e-009  9.5448e-051  8.4384e-027  4.9735e-003
  1.7182e-009  1.8755e-050  1.3238e-026  6.2295e-003
  1.7182e-009  1.5186e-049  5.3381e-026  1.2509e-002
  1.7182e-009  5.1461e-049  1.2043e-025  1.8789e-002
  1.7182e-009  1.2223e-048  2.1438e-025  2.5069e-002
  1.7182e-009  2.3901e-048  3.3524e-025  3.1348e-002
  1.7182e-009  1.9167e-047  1.3431e-024  6.2747e-002
  1.7182e-009  6.4739e-047  3.0236e-024  9.4146e-002
  1.7182e-009  1.5352e-046  5.3768e-024  1.2554e-001
  1.7182e-009  2.9991e-046  8.4026e-024  1.5694e-001
  1.7182e-009  2.4004e-045  3.3621e-023  3.1394e-001
  1.7182e-009  8.1028e-045  7.5655e-023  4.7093e-001
  1.7182e-009  1.9208e-044  1.3451e-022  6.2792e-001
  1.7182e-009  3.7518e-044  2.1017e-022  7.8491e-001
  1.7182e-009  3.0017e-043  8.4074e-022  1.5699e+000
  1.7182e-009  1.0131e-042  1.8917e-021  2.3548e+000
  1.7182e-009  2.4015e-042  3.3631e-021  3.1398e+000
  1.7182e-009  4.6904e-042  5.2548e-021  3.9248e+000
  1.7182e-009  3.7524e-041  2.1020e-020  7.8496e+000
  1.7182e-009  1.2664e-040  4.7294e-020  1.1774e+001
  1.7182e-009  3.0020e-040  8.4079e-020  1.5699e+001
  1.7182e-009  5.8632e-040  1.3137e-019  1.9624e+001
  1.7182e-009  4.6906e-039  5.2550e-019  3.9248e+001
  1.7182e-009  1.5831e-038  1.1824e-018  5.8872e+001
  1.7182e-009  3.7525e-038  2.1020e-018  7.8496e+001
  1.7182e-009  7.3291e-038  3.2844e-018  9.8121e+001
  1.7182e-009  5.8633e-037  1.3137e-017  1.9624e+002
  1.7182e-009  1.9789e-036  2.9559e-017  2.9436e+002
  1.7182e-009  4.6906e-036  5.2550e-017  3.9248e+002
  1.7182e-009  9.1613e-036  8.2109e-017  4.9060e+002
  1.7182e-009  7.3291e-035  3.2844e-016  9.8121e+002
  1.7182e-009  2.4736e-034  7.3898e-016  1.4718e+003
  1.7182e-009  5.8633e-034  1.3137e-015  1.9624e+003
  1.7182e-009  1.1452e-033  2.0527e-015  2.4530e+003
  1.7182e-009  9.1614e-033  8.2109e-015  4.9060e+003
  1.7182e-009  3.0920e-032  1.8475e-014  7.3590e+003
  1.7182e-009  7.3291e-032  3.2844e-014  9.8121e+003
  1.7182e-009  1.4315e-031  5.1318e-014  1.2265e+004
  1.7182e-009  1.1452e-030  2.0527e-013  2.4530e+004
  1.7182e-009  3.8649e-030  4.6186e-013  3.6795e+004
  1.7182e-009  9.1614e-030  8.2109e-013  4.9060e+004
  1.7182e-009  1.7893e-029  1.2830e-012  6.1325e+004
  1.7182e-009  1.4315e-028  5.1318e-012  1.2265e+005
  1.7182e-009  4.8312e-028  1.1547e-011  1.8398e+005
  1.7182e-009  1.1452e-027  2.0527e-011  2.4530e+005
  1.7182e-009  2.2367e-027  3.2074e-011  3.0663e+005
  1.7182e-009  1.7893e-026  1.2830e-010  6.1325e+005
  1.7182e-009  6.0390e-026  2.8866e-010  9.1988e+005
  1.7182e-009  1.4315e-025  5.1318e-010  1.2265e+006
  1.7182e-009  2.7958e-025  8.0185e-010  1.5331e+006
  1.7182e-009  2.2367e-024  3.2074e-009  3.0663e+006
  1.7182e-009  7.5487e-024  7.2166e-009  4.5994e+006
  1.7182e-009  1.7893e-023  1.2830e-008  6.1325e+006
  1.7182e-009  3.4948e-023  2.0046e-008  7.6657e+006
  1.7182e-009  2.7958e-022  8.0185e-008  1.5331e+007
  1.7182e-009  9.4359e-022  1.8042e-007  2.2997e+007
  1.7182e-009  2.2367e-021  3.2074e-007  3.0663e+007
  1.7182e-009  4.3685e-021  5.0115e-007  3.8328e+007
  1.7182e-009  3.4948e-020  2.0046e-006  7.6657e+007
  1.7182e-009  1.1795e-019  4.5104e-006  1.1499e+008
  1.7182e-009  2.7958e-019  8.0185e-006  1.5331e+008
  1.7182e-009  5.4606e-019  1.2529e-005  1.9164e+008
  1.7182e-009  4.3685e-018  5.0115e-005  3.8328e+008
  1.7182e-009  1.4744e-017  1.1276e-004  5.7493e+008
  1.7182e-009  3.4948e-017  2.0046e-004  7.6657e+008
  1.7182e-009  6.8257e-017  3.1322e-004  9.5821e+008
  1.7182e-009  5.4606e-016  1.2529e-003  1.9164e+009
  1.7182e-009  1.8429e-015  2.8190e-003  2.8746e+009
  1.7182e-009  4.3685e-015  5.0115e-003  3.8328e+009
  1.7182e-009  8.5322e-015  7.8305e-003  4.7910e+009
  1.7182e-009  6.8257e-014  3.1322e-002  9.5821e+009
  1.7182e-009  2.3037e-013  7.0475e-002  1.4373e+010
  1.7182e-009  5.4606e-013  1.2529e-001  1.9164e+010
  1.7182e-009  1.0665e-012  1.9576e-001  2.3955e+010
  1.7182e-009  8.5324e-012  7.8306e-001  4.7912e+010
  1.7182e-009  2.8797e-011  1.7619e+000  7.1871e+010
  1.7182e-009  6.8260e-011  3.1324e+000  9.5833e+010
  1.7182e-009  1.3332e-010  4.8946e+000  1.1980e+011
  1.7182e-009  1.0671e-009  1.9584e+001  2.3975e+011
  1.7182e-009  3.6020e-009  4.4087e+001  3.5999e+011
  1.7182e-009  8.5407e-009  7.8433e+001  4.8067e+011
  1.7182e-009  1.6690e-008  1.2266e+002  6.0194e+011
  1.7182e-009  7.0159e-008  3.1998e+002  9.7815e+011
  1.7182e-009  1.8476e-007  6.1265e+002  1.3655e+012
  1.7182e-009  3.8467e-007  1.0040e+003  1.7686e+012
  1.7182e-009  6.9524e-007  1.4986e+003  2.1918e+012
  1.7182e-009  1.1431e-006  2.1021e+003  2.6402e+012
  1.7182e-009  1.7561e-006  2.8213e+003  3.1189e+012
  1.7182e-009  2.5640e-006  3.6645e+003  3.6332e+012
  1.7182e-009  3.5993e-006  4.6414e+003  4.1892e+012
  1.7183e-009  4.8968e-006  5.7631e+003  4.7931e+012
  1.7183e-009  6.4941e-006  7.0425e+003  5.4519e+012
  1.7183e-009  8.4324e-006  8.4942e+003  6.1732e+012
  1.7183e-009  1.0757e-005  1.0135e+004  6.9651e+012
  1.7183e-009  1.3517e-005  1.1983e+004  7.8368e+012
  1.7183e-009  1.6768e-005  1.4061e+004  8.7982e+012
  1.7183e-009  2.0569e-005  1.6391e+004  9.8603e+012
  1.7183e-009  2.4986e-005  1.9000e+004  1.1035e+013
  1.7183e-009  3.0095e-005  2.1919e+004  1.2337e+013
  1.7183e-009  3.5975e-005  2.5181e+004  1.3779e+013
  1.7183e-009  4.2717e-005  2.8822e+004  1.5379e+013
  1.7183e-009  5.0421e-005  3.2885e+004  1.7156e+013
  1.7183e-009  5.9198e-005  3.7416e+004  1.9128e+013
  1.7183e-009  6.9172e-005  4.2467e+004  2.1320e+013
  1.7184e-009  8.0480e-005  4.8096e+004  2.3755e+013
  1.7184e-009  9.3274e-005  5.4367e+004  2.6461e+013
  1.7184e-009  1.0772e-004  6.1352e+004  2.9471e+013
  1.7184e-009  1.2402e-004  6.9131e+004  3.2817e+013
  1.7185e-009  1.4236e-004  7.7791e+004  3.6537e+013
  1.7185e-009  1.6299e-004  8.7434e+004  4.0675e+013
  1.7185e-009  1.8617e-004  9.8168e+004  4.5278e+013
  1.7186e-009  2.1218e-004  1.1012e+005  5.0397e+013
  1.7187e-009  2.4134e-004  1.2341e+005  5.6090e+013
  1.7187e-009  2.7401e-004  1.3821e+005  6.2423e+013
  1.7188e-009  3.1059e-004  1.5469e+005  6.9466e+013
  1.7189e-009  3.5151e-004  1.7301e+005  7.7300e+013
  1.7190e-009  3.9727e-004  1.9341e+005  8.6013e+013
  1.7191e-009  4.4840e-004  2.1610e+005  9.5702e+013
  1.7191e-009  4.8142e-004  2.3071e+005  1.0194e+014
  1.7192e-009  5.1667e-004  2.4626e+005  1.0857e+014
  1.7193e-009  5.5429e-004  2.6283e+005  1.1564e+014
  1.7194e-009  5.9444e-004  2.8048e+005  1.2317e+014
回复 支持 反对

使用道具 举报

zhihaikou
 楼主| 发表于 2008-9-22 14:08 | 显示全部楼层
原帖由 messenger 于 2008-9-22 14:02 发表

                               
登录/注册后可看大图

自己看看帮助trapz怎么用，很简单

最讨厌告诉命令以后，再问命令怎么用了
那个命令需要给出函数

我这个是求出来的数值解

最后得到这个数值解在该区间上的积分
回复 支持 反对

使用道具 举报

messenger
发表于 2008-9-22 14:16 | 显示全部楼层
trapz就是求数值解积分的

你在哪看到trapz命令需要给出函数的？？？:@o


原帖由 zhihaikou 于 2008-9-22 14:08 发表

                               
登录/注册后可看大图

那个命令需要给出函数

我这个是求出来的数值解

最后得到这个数值解在该区间上的积分

评分

1

查看全部评分

回复 支持 反对

使用道具 举报

zhihaikou
 楼主| 发表于 2008-9-22 14:18 | 显示全部楼层
原帖由 messenger 于 2008-9-22 14:16 发表

                               
登录/注册后可看大图

trapz就是求数值解积分的

你在哪看到trapz命令需要给出函数的？？？:@o
谢谢老兄了 明白了 搞定了
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-4 13:26 , Processed in 0.061659 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表