| 3 | 1/1 | 返回列表 |
| 查看: 414 | 回復(fù): 2 | ||
DZ木林森新蟲 (初入文壇)
|
[求助]
Matlab 解方程問題求教!。
|
|
代碼1,求解sym類型參數(shù)方程中未知數(shù)。 以下是MATLAB程序代碼,我運(yùn)行之后提示“FZERO cannot continue because user supplied function_handle ==> @(x)(sum((y'./beta')./(1/alpha_t*ones(size(gamma))-1./gamma))-K) failed with the error below.”,作為新手有點(diǎn)看不懂,不會改了,求大神幫忙,先謝謝了。 需要輸入的內(nèi)容如下: the number of the variables:8 K:12.2 y:[0.04 0.44 6.16 24.98 37.23 25.43 5.69 0.05] alpha_t:0.527 d:[60.7 28.3 13.2 6.1 2.9 1.3 0.6 0.3] 程序代碼 clc;clear;close all n=input('the number of the variables:'); K=input('K:'); y=input('y:'); alpha_t=input('alpha_t:'); d=input('d:'); if (length(y)~=n||length(d)~=n) disp('please input the same number as n'); break; end syms beta; gamma=sym(zeros(n,1)); for i=1:n sum1=0; sum2=0; if i==1 for j=i+1:n aij=sqrt(1-(1-d(j)/d(i))^1.02); sum2=sum2+aij/beta*y(j); end gamma(1)=1/(y(1)/beta+sum2+sum1); end if i==n for j=1:i-1 bij=1-(1-d(i)/d(j))^1.5; sum1=sum1+(1-bij*(1-1/beta))*y(j); gamma(i)=1/(y(i)/beta+sum1+sum2); end end if i<n&&1<i for j=1:i-1 bij=1-(1-d(i)/d(j))^1.5; sum1=sum1+(1-bij*(1-1/beta))*y(j); end for j=i+1:n aij=sqrt(1-(1-d(j)/d(i))^1.02); sum2=sum2+aij/beta*y(j); end gamma(i)=1/(y(i)/beta+sum1+sum2); end end fun=@(x)(sum((y'./beta')./(1/alpha_t*ones(size(gamma))-1./gamma))-K); x1=0; x2=1; x=fzero(fun,0.6) abs(fun(x)) beta1=x; disp(['beta:']) beta1; |
鐵蟲 (初入文壇)

新蟲 (初入文壇)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 085600材料與化工 +8 | 安全上岸! 2026-03-16 | 8/400 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 13/650 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +4 | 一瓶苯甲酸 2026-03-14 | 4/200 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +10 | chaojifeixia 2026-03-19 | 11/550 |
|
|
[考研] 265求調(diào)劑 +8 | 梁梁校校 2026-03-17 | 8/400 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
[基金申請] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|