| 查看: 3459 | 回復(fù): 10 | ||
鐘聲入翠微銅蟲 (小有名氣)
|
[求助]
用matlab求解方程式 已有1人參與
|
|
|
假設(shè)a和b已知,a,b,x,y>0,求解x和y a=[1+(0.0172*x)^2]^(-y/2) b=y*arctan(0.0172*x) 希望可以提供程序代碼,或者解!非常感謝! |
版主 (知名作家)
|
這個需要嘗試,可以用 fsolve 函數(shù)隨機(jī)賦值求解。 結(jié)果不理想就要考慮方程是否正確,或者a,b取值是否合理。 function example_0708 clear; a=0.709133065; b=0.08555; x0 = [1,1]; x = fsolve(@root2d,x0); fprintf('\tx=%.8f\n',x(1)); fprintf('\ty=%.8f\n',x(2)); function F = root2d(x) F(1) = (1+(0.0172*x(1))^2)^(-x(2)/2)-a; F(2) = x(2)*atan(0.0172*x(1))-b; end end |

版主 (知名作家)
|
這個盡量不要做伸手黨,求解符號方程組,看看solve函數(shù)用法 發(fā)自小木蟲Android客戶端 |

銅蟲 (小有名氣)
送紅花一朵 |
我是去百度了solve用法,但是運(yùn)行的時(shí)候一直顯示錯誤, 我最先給a,b賦值進(jìn)行計(jì)算,但是一直顯示運(yùn)算符錯誤 syms x y; >>[x ,y]=solve(' 0.709133=[1+(0.0172*x)^2]^(-y/2)','0.08555=y*arctan(0.0172*x)'); 然后使用了未賦值 clc,clear syms a b y x sola=solve(a==(1+(0.0172*x)^2)^(-y/2),x) %待求解的變量是x sol=solve(b==y*atan(0.0172*x),y) %待求解的變量是y 還是運(yùn)算不出想要的東西 因?yàn)橹粫稽c(diǎn)點(diǎn)matlab,就一直解決不了問題。 |
版主 (知名作家)

銅蟲 (小有名氣)
銅蟲 (小有名氣)
送紅花一朵 |
我自己又代入了一下,得到這樣的結(jié)果,是什么意思呢? Warning: The solutions are parameterized by the symbols: z, z1. To include parameters and conditions in the solution, specify the 'ReturnConditions' option. > In solve>warnIfParams (line 510) In solve (line 367) Warning: The solutions are valid under the following conditions: 1/((1849*z^2)/6250000 + 1)^(z1/2) - 6387302814580195/9007199254740992 == 0 & z1*atan((43*z)/2500) - 1711/20000 == 0 & 0 < z & 0 < z1. To include parameters and conditions in the solution, specify the 'ReturnConditions' option. > In solve>warnIfParams (line 517) In solve (line 367) X = z Y = z1 |
版主 (知名作家)

銅蟲 (小有名氣)
送紅花一朵 |
謝謝。30000,0.05)這個初始值是自己根據(jù)實(shí)際情況設(shè)置,還是隨意設(shè)置的呢?這次我運(yùn)行出來了,但是感覺結(jié)果不是我想象中的那樣。 發(fā)自小木蟲Android客戶端 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級已過 求調(diào)劑 +4 | ccc! 2026-03-03 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 292求調(diào)劑 +9 | yhk_819 2026-02-28 | 9/450 |
|
|
[考研] 材料328求調(diào)劑 +7 | 一個蘿卜02 2026-03-03 | 7/350 |
|
|
[考研] 武理材料工程302調(diào)劑 +3 | Doleres 2026-03-03 | 6/300 |
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
[考研] 材料專碩290求調(diào)劑 +3 | 杰尼龜aaa 2026-03-04 | 3/150 |
|
|
[考研] 一志愿西交大材料學(xué)碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
[考研] 322,求調(diào)劑 +3 | 菜菜愛玩 2026-03-04 | 3/150 |
|
|
[考研] 理學(xué),工學(xué),農(nóng)學(xué)調(diào)劑,少走彎路,這里歡迎您! +8 | likeihood 2026-03-02 | 11/550 |
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
[考研] 290求調(diào)劑 +9 | ErMiao1020 2026-03-02 | 9/450 |
|
|
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5 | Qwertyuop 2026-03-01 | 5/250 |
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
[考研]
材料270求調(diào)劑
6+6
|
Eiiiio 2026-03-01 | 11/550 |
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +4 | 青科11 2026-03-02 | 4/200 |
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
[考研] 調(diào)劑材料學(xué)碩 +4 | 詞凝Y 2026-03-02 | 4/200 |
|
|
[考研] 0856調(diào)劑 +10 | 劉夢微 2026-02-28 | 10/500 |
|
|
[考研] 281求調(diào)劑 +5 | 2026計(jì)算機(jī)_誠心 2026-03-01 | 8/400 |
|
|
[考研] 化工299分求調(diào)劑 一志愿985落榜 +5 | 嘻嘻(*^ω^*) 2026-03-01 | 5/250 |
|