24小时热门版块排行榜    

查看: 428  |  回复: 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的回帖

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的回帖

宁静心空

银虫 (著名写手)

月只蓝: 提醒:为了表示对他人回复的尊重,请及时给予应助者金币 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的回帖
相关版块跳转 我要订阅楼主 宁静心空 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +3 西南交通专材355 2026-03-19 3/150 2026-03-20 09:34 by 每天只摆一小会
[考研] 304求调剂 +5 曼殊2266 2026-03-18 5/250 2026-03-20 09:00 by ZHANG0tao
[考研] 294求调剂材料与化工专硕 +14 陌の森林 2026-03-18 14/700 2026-03-19 22:38 by 学员8dgXkO
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by !本暗一次!
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +5 吃吃吃才有意义 2026-03-19 5/250 2026-03-19 16:18 by 30660438
[考研] 324分 085600材料化工求调剂 +3 llllkkkhh 2026-03-18 3/150 2026-03-19 14:22 by houyaoxu
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 一志愿天大材料与化工(085600)总分338 +5 蔡大美女 2026-03-13 5/250 2026-03-19 10:44 by 是小刘呀~
[考研] 0703化学调剂,求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
信息提示
请填处理意见