| 查看: 3456 | 回復(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) 希望可以提供程序代碼,或者解!非常感謝! |
版主 (知名作家)
|
這個(gè)需要嘗試,可以用 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 |

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

銅蟲 (小有名氣)
送紅花一朵 |
我是去百度了solve用法,但是運(yùn)行的時(shí)候一直顯示錯(cuò)誤, 我最先給a,b賦值進(jìn)行計(jì)算,但是一直顯示運(yùn)算符錯(cuò)誤 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)橹粫?huì)一點(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)這個(gè)初始值是自己根據(jù)實(shí)際情況設(shè)置,還是隨意設(shè)置的呢?這次我運(yùn)行出來了,但是感覺結(jié)果不是我想象中的那樣。 發(fā)自小木蟲Android客戶端 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
|
百度網(wǎng)盤 |
360云盤 |
千易網(wǎng)盤 |
華為網(wǎng)盤
在新窗口頁面中打開自己喜歡的網(wǎng)盤網(wǎng)站,將文件上傳后,然后將下載鏈接復(fù)制到帖子內(nèi)容中就可以了。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304分材料專碩求調(diào)劑 +8 | qiuzhigril 2026-03-03 | 10/500 |
|
|---|---|---|---|---|
|
[論文投稿]
100+3
|
Stray2021 2026-03-03 | 3/150 |
|
|
[考研] 一志愿西交大材料學(xué)碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
[考研] 0857求調(diào)劑 +3 | CDDD1013 2026-03-04 | 3/150 |
|
|
[考研]
材料325求調(diào)劑
30+5
|
mariusuki 2026-03-02 | 10/500 |
|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +5 | 1一11 2026-03-03 | 5/250 |
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
[考研] 085700資環(huán)求調(diào)劑,初始279,六級(jí)已過,英語能力強(qiáng) +3 | 085700資環(huán)調(diào)劑 2026-03-03 | 4/200 |
|
|
[考研]
|
glwshine 2026-03-02 | 5/250 |
|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級(jí)已過 求調(diào)劑 +3 | ccc! 2026-03-03 | 3/150 |
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
[考研] 0856材料求調(diào)劑 +12 | hyf hyf hyf 2026-02-28 | 13/650 |
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
[考研] 材料與化工328求調(diào)劑 +3 | 。,。,。,。i 2026-03-02 | 3/150 |
|
|
[考研] 264求調(diào)劑 +4 | 巴拉巴拉根556 2026-02-28 | 4/200 |
|
|
[考研] 299求調(diào)劑 +3 | Y墨明棋妙Y 2026-02-28 | 5/250 |
|
|
[考研] 311求調(diào)劑 +6 | 亭亭亭01 2026-03-01 | 6/300 |
|
|
[考研] 307求調(diào)劑 +4 | 73372112 2026-02-28 | 6/300 |
|
|
[考研] 304求調(diào)劑 +3 | 52hz~~ 2026-02-28 | 5/250 |
|