| 3 | 1/1 | 返回列表 |
| 查看: 430 | 回復(fù): 2 | ||
寧靜心空銀蟲 (著名寫手)
|
[求助]
matlab 積分計(jì)算無法得到結(jié)果 已有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 |
|
問題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的代碼不做大改動(dòng),則需把上面兩行代碼改為: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函數(shù)相當(dāng)于定義了f1(a,m,n,s,x,xr) =...這樣一個(gè)函數(shù),這樣后面就可以直接使用了,如f1(0.45,1,1,1,0.01,xr),注意是6個(gè)變量。 或者可以不用inline函數(shù),將所有的f1(s,a,n,m,x)的形式均去掉,而使用subs賦值函數(shù) 問題2:y3,y4的積分是沒有解析解的,要用數(shù)值解quad去求 |

銀蟲 (著名寫手)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
簡(jiǎn)木ChuFront 2026-03-19 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工 +7 | 安全上岸! 2026-03-16 | 7/350 |
|
|
[考研] 材料080500調(diào)劑求收留 +8 | 一顆meteor 2026-03-13 | 8/400 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 14/700 |
|
|
[考研] 304求調(diào)劑 +5 | 曼殊2266 2026-03-18 | 5/250 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|