本帖產(chǎn)生 1 個(gè) 程序強(qiáng)帖，點(diǎn)擊這里進(jìn)行查看

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

zhmdream

應(yīng)助: 0 (幼兒園)
金幣: 2535.9
紅花: 1
帖子: 776
在線: 292.2小時(shí)
蟲(chóng)號(hào): 591487
注冊(cè): 2008-09-02
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

[求助] Matlab兩個(gè)小問(wèn)題

您好，請(qǐng)教兩個(gè)問(wèn)題：
（1）我用Matlab解一個(gè)一元4次方程，得到4個(gè)解。怎么才能讓其輸出我想要的實(shí)數(shù)解，如第一個(gè)解（或者解的條件滿足：y>400)。

y =
  1.0e+002 *
4.3541
2.9607
  -3.6574 + 5.2191i
  -3.6574 - 5.2191i

(2)第2個(gè)問(wèn)題是在第1個(gè)問(wèn)題基礎(chǔ)上的。我的目的是為了畫(huà)出x與y的關(guān)系圖，所以我必須得對(duì)連續(xù)的x進(jìn)行求解，然后求出y，作圖。以一個(gè)簡(jiǎn)單例子吧：

x=1:1:10; %x取1到10
syms y;
a=y+1;
b=y-1;

f=a*b+a+b*x+y+6*x+2

%解方程f=0
equation = [char(f) ' = ' num2str(0,9)];
A=solve(equation, 'y')

double(A)

plot(x,y)

但是運(yùn)行出錯(cuò)。有沒(méi)有其他方法可以解決？謝謝！

回復(fù)此樓

人生如夢(mèng)，歲月無(wú)情~~~

1樓 2011-10-11 20:24:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sudo

木蟲(chóng) (正式寫手)

程序強(qiáng)帖: 16
應(yīng)助: 6 (幼兒園)
金幣: 1297.6
散金: 1486
紅花: 20
帖子: 588
在線: 641小時(shí)
蟲(chóng)號(hào): 1211394
注冊(cè): 2011-02-24
性別: GG
專業(yè): 文藝美學(xué)

★
jjdg(金幣+1): 好建議 2011-10-15 01:56:10

要善用help命令和google啊...

上面這個(gè)程序里面，需要注意的手法只有一個(gè)：

如何把一個(gè)符號(hào)表達(dá)式，轉(zhuǎn)化為函數(shù)句柄（因?yàn)镸ATLAB里面很多解方程/優(yōu)化函數(shù)里面都有一個(gè)“函數(shù)句柄”的參數(shù)）

贊一下(1人)

回復(fù)此樓

9樓2011-10-14 13:07:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 11 個(gè)回答

sudo

木蟲(chóng) (正式寫手)

程序強(qiáng)帖: 16
應(yīng)助: 6 (幼兒園)
金幣: 1297.6
散金: 1486
紅花: 20
帖子: 588
在線: 641小時(shí)
蟲(chóng)號(hào): 1211394
注冊(cè): 2011-02-24
性別: GG
專業(yè): 文藝美學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★
xzhdty(金幣+2): 歡迎常來(lái)程序語(yǔ)言看看 2011-10-12 15:31:47
微塵、夢(mèng)想(金幣+2, 程序強(qiáng)帖+1): sudo總是很熱心的，祝福哦！ 2011-10-12 20:41:46
zhmdream(金幣+40): 謝謝啊，我試試 2011-10-12 20:45:04
zhmdream(金幣+10): 謝謝了，呵呵 2011-10-14 22:14:37

你不用每次都放那么多分上來(lái)的...

(1)使用isreal函數(shù)和循環(huán)

CODE:

for i=1:length(y)

    if isreal(y(i)) && y(i)>400 %判斷y(i)是否是實(shí)數(shù)，并且是否大于400

        fprintf(1, '%.9f\n', y(i));

    end

end

(2)

CODE:

syms x y; %注意用到的符號(hào)都要聲明一下

a=y+1;

b=y-1;

f=a*b+a+b*x+y+6*x+2;

%解方程f=0

equation = [char(f) ' = 0']; 

A=solve(equation, 'y');

%上面解得的結(jié)果是

%A =

%  - x/2 - (x^2 - 16*x - 4)^(1/2)/2 - 1

%   (x^2 - 16*x - 4)^(1/2)/2 - x/2 - 1

%對(duì)于符號(hào)式子，使用subs函數(shù)代入數(shù)值

%注意上面的方程因?yàn)橛虚_(kāi)方，所以不一定能解出實(shí)數(shù)解，這里直接取一段有實(shí)數(shù)解的

x = 20:80;

y1 = zeros(1, length(x)); %習(xí)慣上循環(huán)前先給定存儲(chǔ)空間

y2 = zeros(1, length(x));

for i=1:length(x)

    y1(i) = subs(A(1), x(i));

    y2(i) = subs(A(2), x(i));

end

plot(x, y1, 'b', x, y2, 'r');

贊一下(2人)

回復(fù)此樓

2樓2011-10-12 10:12:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sudo

木蟲(chóng) (正式寫手)

程序強(qiáng)帖: 16
應(yīng)助: 6 (幼兒園)
金幣: 1297.6
散金: 1486
紅花: 20
帖子: 588
在線: 641小時(shí)
蟲(chóng)號(hào): 1211394
注冊(cè): 2011-02-24
性別: GG
專業(yè): 文藝美學(xué)

★ ★
微塵、夢(mèng)想(金幣+2): 謝謝參與應(yīng)助~ 2011-10-12 20:42:05

其實(shí)第二個(gè)問(wèn)題里面，算出A之后，直接

CODE:

x = 20:80;

Y= subs(A, x);

plot(x, Y(1, :), 'b', x, Y(2, :), 'r');

就可以的...

不過(guò)這么寫的話覺(jué)得你可能不容易理解

贊一下(1人)

回復(fù)此樓

3樓2011-10-12 10:16:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhmdream

木蟲(chóng) (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2535.9
紅花: 1
帖子: 776
在線: 292.2小時(shí)
蟲(chóng)號(hào): 591487
注冊(cè): 2008-09-02
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

2樓: Originally posted by sudo at 2011-10-12 10:12:19:

你不用每次都放那么多分上來(lái)的...

(1)使用isreal函數(shù)和循環(huán)
[code]
for i=1:length(y)
if isreal(y(i)) && y(i)>400 %判斷y(i)是否是實(shí)數(shù)，并且是否大于400
fprintf(1, '%.9f ...

呵呵，謝謝幫忙。問(wèn)題解決差不多了，但是還有些小問(wèn)題：

問(wèn)題（2）的方程當(dāng)時(shí)是為了方便，隨便設(shè)的。我要解的是一元四次方程，解的結(jié)果有兩個(gè)實(shí)解和兩個(gè)虛解。舍去虛解，并只有實(shí)解>400，才是我所要求的解，并用于作圖。下面我的設(shè)計(jì)與思路：

%%%之前程序省略
%求解Tg
A=solve(equation, 'Tg');
Tp=450:1:453;    %變量Tp取值，并將Tp代入，求出Tg
for i=1:length(Tp)
Tg1(i) = subs(A(1), Tp(i));
Tg2(i) = subs(A(2), Tp(i));
Tg3(i) = subs(A(3), Tp(i));
Tg4(i) = subs(A(4), Tp(i));
end
%得到兩實(shí)解與兩虛解

%現(xiàn)在對(duì)這兩實(shí)解與兩虛解，進(jìn)行比較。大于400的實(shí)數(shù)解才是我最終想要的
for i=1:4
   j=1:length(Tg)
if isreal(Tg(i)(j)) && Tg(i)(j)>400 %判斷Tg(i)是否是實(shí)數(shù)，并且是否大于400
      fprintf(1, '%.9f\n', Tg(i)(j));
end
end

plot(Tp, Tg) %對(duì)Tp與Tg作圖

我運(yùn)行時(shí)，到了第2部分就出錯(cuò)。謝謝！

贊一下

回復(fù)此樓

人生如夢(mèng)，歲月無(wú)情~~~

4樓2011-10-13 20:23:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 22 350 本科985求調(diào)劑，求老登收留 +4	李軼男003 2026-03-20	4/200	2026-03-26 16:05 by 哇啦啦啦xtj
[考研] 一志愿河工大 081700 276求調(diào)劑 +4	地球繞著太陽(yáng)轉(zhuǎn) 2026-03-23	4/200	2026-03-26 14:27 by zzll406
[考研] 材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦 +13	adaie 2026-03-19	13/650	2026-03-26 14:05 by 小小麥片
[考研] 總分293求調(diào)劑 +6	加一一九 2026-03-25	8/400	2026-03-26 13:30 by yujianx
[考研] 309求調(diào)劑 +4	gajsj 2026-03-25	5/250	2026-03-26 00:27 by Dyhoer
[考研] 311求調(diào)劑 +4	勇敢的小吳 2026-03-20	4/200	2026-03-25 18:12 by xcjcqu
[考研] 招08考數(shù)學(xué) +8	laoshidan 2026-03-20	17/850	2026-03-25 17:52 by 一個(gè)紅太陽(yáng)
[考研] 311求調(diào)劑 +3	冬十三 2026-03-24	3/150	2026-03-24 21:31 by peike
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +7	天天好運(yùn)來(lái)上岸?/a> 2026-03-24	7/350	2026-03-24 20:26 by peike
[考研] 化工專碩求調(diào)劑 +3	question挽風(fēng) 2026-03-24	3/150	2026-03-24 18:48 by jhhcooi
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +5	leaflight 2026-03-24	5/250	2026-03-24 16:25 by laoshidan
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +6	zxts12 2026-03-21	9/450	2026-03-23 18:37 by macy2011
[考研] 276求調(diào)劑 +3	YNRYG 2026-03-21	4/200	2026-03-23 08:31 by 醉在風(fēng)里
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +4	酥酥魚(yú).. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +4	幸運(yùn)的醬醬 2026-03-22	5/250	2026-03-22 20:15 by edmund7
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5	vv迷 2026-03-21	8/400	2026-03-22 14:20 by ColorlessPI
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 329求調(diào)劑 +9	想上學(xué)吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] 295復(fù)試調(diào)劑 +8	簡(jiǎn)木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo