| 查看: 386 | 回复: 2 | |||
[求助]
matlab 积分计算无法得到结果 已有1人参与
|
|
function FXX clear;clc; syms ks s a n m x y xr=3000; f1(s,a,n,m,x)=s*(1-log(1+x/xr)/log(1+1e6/xr))/(log(exp(1)+(x/a)^n))^m; % f1(s,a,n,m,x)=s+a+n+m+x; f2(s,a,n,m,x)=diff(f1,x); y1=(f1(0.45,1,1,1,y)-f1(0.45,1,1,1,0.01))*f2(0.45,1,1,1,y)/y^2 y2=(f1(0.45,1,1,1,y)-0.45)*f2(0.45,1,1,1,y)/y^2 y3=int(y1,y,0.01,1e6) y4=int(y2,y,1,1e6) return |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 为什么matlab算出的数不正确，在线等
已经有5人回复
- matlab “Inner matrix dimensions must agree.”一直不知道哪里出错了 求助
已经有6人回复
- matlab求函数积分出错
已经有12人回复
- 对于符号积分matlab积分积不出来怎么解决+有代码
已经有14人回复
- matlab循环计算转为矩阵 30金币求助
已经有7人回复
- 请教大家一个反卷积的问题~~帮忙推导或者matlab编程计算~~
已经有5人回复
- matlab计算面积求助！！难
已经有8人回复
- asp.net调用matlab的制作网站，在vs2010上正常运行，但在服务器IIS无法显示计算的结果
已经有3人回复
- Matlab如何实现这样的计算：赋值——计算——判断结果——重复
已经有10人回复
- matlab 积分求助
已经有8人回复
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 关于Matlab算符号积分
已经有6人回复
- 求教，matlab并行运算加速
已经有12人回复
- matlab求积分，出现奇点怎么办？
已经有8人回复
- matlab求数值积分？
已经有10人回复
- 【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？
已经有3人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】请教Matlab高手根据积分画图问题
已经有4人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
宁静心空(月只蓝代发): 金币+10, 感谢应助！ 2017-02-05 11:00:58
感谢参与，应助指数 +1
宁静心空(月只蓝代发): 金币+10, 感谢应助！ 2017-02-05 11:00:58
|
问题1：代码“f1(s,a,n,m,x)=...;f2(s,a,n,m,x)=...;”中f1(s,a,n,m,x)和f2(s,a,n,m,x)这样的写法有误，若想后面y1,y2的代码不做大改动，则需把上面两行代码改为：f1=inline('s*(1-log(1+x/xr)/log(1+1e6/xr))/(log(exp(1)+(x/a)^n))^m'); f2=inline('diff((s*(1-log(1+x/xr)/log(1+1e6/xr))/(log(exp(1)+(x/a)^n))^m),x)');其中，inline函数相当于定义了f1(a,m,n,s,x,xr) =...这样一个函数，这样后面就可以直接使用了，如f1(0.45,1,1,1,0.01，xr)，注意是6个变量。 或者可以不用inline函数，将所有的f1(s,a,n,m,x)的形式均去掉，而使用subs赋值函数 问题2:y3,y4的积分是没有解析解的，要用数值解quad去求 |
2楼2015-01-07 22:04:13
3楼2015-01-08 11:02:26