 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 1603 | 回復: 9 | |||
chuzhenhui金蟲 (小有名氣)
|
[交流]
【求助】如何解一個含定積分的方程組
|
|
我想解一個方程組，表達式如下： x1+x2-0.091=0; x1-a/(1-k1*a+c1*k1*a)=0; x2-f(t)-0.0778=0; t=quad(@(a)((k1^2*(c1-1).*a.^2+1)./((1-k1.*a).^2.*(c1*k1.*a-k1.*a+1).^2)),0.1404,a); 其中x1,x2,t,a是四個未知數(shù)，其他是常數(shù)。f(t)是關于未知數(shù)t的表達式。 這個方程組是我想解的方程組的簡化，只是用來表示一個類似的方程組表達式，最終我是想得到t的值。 我用solve函數(shù)解： [x1,x2,t,a]=solve('x1+x2-0.091=0','x1-a/(1-k1*a+c1*k1*a)=0','x2-f(t)-0.0778=0','t=quad(@(a)((k1^2*(c1-1).*a.^2+1)./((1-k1.*a).^2.*(c1*k1.*a-k1.*a+1).^2)),0.1404,a)') 但是解不出來，matlab報錯說我其中的定積分的表達式 t=quad(@(a)((k1^2*(c1-1).*a.^2+1)./((1-k1.*a).^2.*(c1*k1.*a-k1.*a+1).^2)),0.1404,a); is not a valid expression or equation. 意思應該是說這個方程的表達式有問題。 我想請教各位的是這種含有定積分的方程組怎么解的？是化成其他形式，用int 積分好像也不行，還是用別的函數(shù)解？能否給我舉個簡單的例子，各位大俠拜托了，謝謝！��！ |
回復此樓» 猜你喜歡
070300化學求調劑
已經(jīng)有3人回復
-
311求調劑
已經(jīng)有3人回復
-
生物學調劑
已經(jīng)有5人回復
-
269專碩求調劑
已經(jīng)有6人回復
-
08工科 320總分 求調劑
已經(jīng)有8人回復
-
尋找調劑
已經(jīng)有3人回復
-
291 求調劑
已經(jīng)有3人回復
-
一志愿北京化工大學070300 學碩336求調劑
已經(jīng)有8人回復
-
一志愿東華大學化學070300，求調劑
已經(jīng)有6人回復
-
318求調劑
已經(jīng)有7人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 解線性代數(shù)方程組
已經(jīng)有10人回復
- 請問如何求解二元一階微分方程組
已經(jīng)有24人回復
- 高金求助matlab解微分方程組
已經(jīng)有12人回復
- 求matlab求解一個方程組的問題
已經(jīng)有6人回復
- 都不知道怎么命名標題：迭代法解方程組？
已經(jīng)有16人回復
- 就解方程組，未知數(shù)不收斂，要怎么辦？誠求！
已經(jīng)有12人回復
- 求一定積分的解析解或數(shù)值解(數(shù)值解則要求給出matlab程序)
已經(jīng)有5人回復
- 【求助】如何解這樣的方程組【有效期至2011年4月30日】
已經(jīng)有9人回復
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題
已經(jīng)有21人回復
- 【求助】求助解決這樣一個非線性方程組
已經(jīng)有8人回復
- 【求助】matlab怎么求解偏微分方程組啊，先謝謝了
已經(jīng)有13人回復
- 一個方程組的解法
已經(jīng)有8人回復
2樓2010-10-27 21:42:43
chuzhenhui
金蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 600.7
- 散金: 35
- 帖子: 129
- 在線: 10.2小時
- 蟲號: 763256
- 注冊: 2009-05-05
- 性別: GG
- 專業(yè): 食品&包裝
3樓2010-10-28 13:24:51
4樓2010-10-28 14:33:10
lijinfeng042
木蟲 (小有名氣)
Matlab
- 仿真EPI: 2
- 應助: 1 (幼兒園)
- 金幣: 2156.1
- 散金: 115
- 帖子: 291
- 在線: 31.5小時
- 蟲號: 1019062
- 注冊: 2010-05-15
- 性別: GG
- 專業(yè): 通信理論與系統(tǒng)
5樓2010-10-28 21:06:56
chuzhenhui
金蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 600.7
- 散金: 35
- 帖子: 129
- 在線: 10.2小時
- 蟲號: 763256
- 注冊: 2009-05-05
- 性別: GG
- 專業(yè): 食品&包裝
|
我用fsolve試了，還是不行，我把我的原方程貼出來，項數(shù)挺多，但方程結構很簡單，你幫我看下吧： function H=myfun(x) syms w1 c1 k1 d p pa s hum w sumy x %藍色字體是常數(shù)，sumy是一個求和項，包含未知數(shù)x(4),但在列方程的時候不知道怎么表達出來，所以就先分開列了。 %4個方程 H=[x(1)+x(2)-0.091; x(3)-quad(((w1*c1*k1*w*d)/(p*pa*s*3600))*((k1^2*(c1-1).*x(4)^2+1)./((hum-x(4)).*(1-k1.*x(4)).^2.*(c1*k1.*x(4)-k1.*x(4)+1).^2)),0.1404,x(4)); x(2)-0.0778-0.0518*sumy; x(1)-(w1*c1*k1.*x(4))./((1-k1.*x(4)).*(1-k1.*x(4)+c1*k1.*x(4)))]; 其中的的quad(((w1*c1*k1*w*d)/(p*pa*s*3600))*((k1^2*(c1-1).*x(4)^2+1)./((hum-x(4)).*(1-k1.*x(4)).^2.*(c1*k1.*x(4)-k1.*x(4)+1).^2)),0.1404,x(4));就是我的定積分部分的方程，我用quad表示了。 %常數(shù)取值 w1=0.0378;c1=9.992;k1=0.9247;p=2.1e-14;pa=2810.94;w=3.7;hum=0.9;s=64;d=30e-4; sumy=0; for n=0:6 sumy=sumy+8/((2*n+1)^2*3.141596^2)*exp(-(2*n+1)^2*3.141596^2*0.5*10e-4*36*x(4)); end x0=[0.001;0.001;0.001;0.001]; [x] = fsolve(@myfun,x0); 這樣解不出來，問題出在哪？麻煩看看吧，謝謝了！ [ Last edited by chuzhenhui on 2010-10-29 at 13:20 ] |
6樓2010-10-29 13:00:34
lijinfeng042
木蟲 (小有名氣)
Matlab
- 仿真EPI: 2
- 應助: 1 (幼兒園)
- 金幣: 2156.1
- 散金: 115
- 帖子: 291
- 在線: 31.5小時
- 蟲號: 1019062
- 注冊: 2010-05-15
- 性別: GG
- 專業(yè): 通信理論與系統(tǒng)
7樓2010-10-29 22:49:49
chuzhenhui
金蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 600.7
- 散金: 35
- 帖子: 129
- 在線: 10.2小時
- 蟲號: 763256
- 注冊: 2009-05-05
- 性別: GG
- 專業(yè): 食品&包裝
8樓2010-10-30 10:17:39
lijinfeng042
木蟲 (小有名氣)
Matlab
- 仿真EPI: 2
- 應助: 1 (幼兒園)
- 金幣: 2156.1
- 散金: 115
- 帖子: 291
- 在線: 31.5小時
- 蟲號: 1019062
- 注冊: 2010-05-15
- 性別: GG
- 專業(yè): 通信理論與系統(tǒng)
9樓2010-10-30 19:31:42
chuzhenhui
金蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 600.7
- 散金: 35
- 帖子: 129
- 在線: 10.2小時
- 蟲號: 763256
- 注冊: 2009-05-05
- 性別: GG
- 專業(yè): 食品&包裝
10樓2010-11-01 09:44:26
相關版塊跳轉
我要訂閱樓主 chuzhenhui 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學071000，求調劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 297求調劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 求調劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 材料學碩333求調劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
|
[考研] 277材料科學與工程080500求調劑 +6 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境總分308求調劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
|
[考研] 307求調劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿中國石油大學（華東） 本科齊魯工業(yè)大學 +3 | 石能偉 2026-03-17 | 3/150 |
|
|
|
[考研] 南京大學化學376求調劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿南昌大學，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 294求調劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
|
[考研] 308求調劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 304求調劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研] 生物學調劑招人�。�！ +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 考研求調劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 308求調劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 302求調劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
|
[考研] 0856求調劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|