24小时热门版块排行榜    

查看: 516  |  回复: 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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

宁静心空

银虫 (著名写手)

月只蓝: 提醒:为了表示对他人回复的尊重,请及时给予应助者金币 2017-02-05 11:01:21
引用回帖:
2楼: Originally posted by wdlgezi at 2015-01-07 22:04:13
问题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 ...

已按照所提问题进行修改,问题已解决。Matlab新手上路,感谢指教,学习了
3楼2015-01-08 11:02:26
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

wdlgezi

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 长年满屏的广告,版主太不责任了。基金也等的急 +4 gltch 2026-08-12 4/200 2026-08-12 16:08 by eulota
[基金申请] filecode +13 documentary 2026-08-10 15/750 2026-08-12 15:56 by 云上清扬
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
请填处理意见