版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

本帖產(chǎn)生 1 個程序強帖，點擊這里進行查看

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

zhmdream

木蟲 (正式寫手)

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

[求助] Matlab兩個小問題

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

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

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

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)

但是運行出錯。有沒有其他方法可以解決？謝謝！

回復(fù)此樓

人生如夢，歲月無情~~~

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

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

sudo

木蟲 (正式寫手)

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

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

其實第二個問題里面，算出A之后，直接

CODE:

x = 20:80;

Y= subs(A, x);

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

就可以的...

不過這么寫的話覺得你可能不容易理解

贊一下(1人)

回復(fù)此樓

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

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

查看全部 11 個回答

sudo

木蟲 (正式寫手)

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

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

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

你不用每次都放那么多分上來的...

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

CODE:

for i=1:length(y)

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

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

    end

end

(2)

CODE:

syms x y; %注意用到的符號都要聲明一下

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

%對于符號式子，使用subs函數(shù)代入數(shù)值

%注意上面的方程因為有開方，所以不一定能解出實數(shù)解，這里直接取一段有實數(shù)解的

x = 20:80;

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

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的回帖

zhmdream

木蟲 (正式寫手)

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

引用回帖:

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

你不用每次都放那么多分上來的...

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

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

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

%%%之前程序省略
%求解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
%得到兩實解與兩虛解

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

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

我運行時，到了第2部分就出錯。謝謝！

贊一下

回復(fù)此樓

人生如夢，歲月無情~~~

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

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

sudo

木蟲 (正式寫手)

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

還是貼出你的全部程序吧

贊一下

回復(fù)此樓

5樓2011-10-13 20:28:33

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

查看全部 11 個回答

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

[考研] 297求調(diào)劑 +11	田洪有 2026-03-26	11/550	2026-03-29 13:14 by peike
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +11	c297914 2026-03-28	11/550	2026-03-29 13:11 by peike
[考研] 壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7	迷人的哈哈 2026-03-28	7/350	2026-03-28 16:47 by 催化大白
[考研] 283求調(diào)劑 +3	A child 2026-03-28	3/150	2026-03-28 15:41 by ms629
[考研] 286求調(diào)劑 +12	PolarBear11 2026-03-26	12/600	2026-03-28 12:14 by zllcz
[考研] 【求調(diào)劑】085601材料工程專碩 \| 總分272 \| +6	腳滑的守法公民 2026-03-27	6/300	2026-03-28 11:02 by gjlllb
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6	烤麥芽 2026-03-27	8/400	2026-03-28 10:40 by 烤麥芽
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7	hanamiko 2026-03-27	7/350	2026-03-28 09:56 by 李上岸0921
[有機交流] 高溫高壓反應(yīng)求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 化學(xué)調(diào)劑 +4	愛吃番茄的旭 2026-03-24	5/250	2026-03-27 17:50 by kiokin
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4	Ripcord上岸 2026-03-27	4/200	2026-03-27 17:00 by zhyzzh
[考研] 348求調(diào)劑 +4	小懶蟲不懶了 2026-03-27	5/250	2026-03-27 12:47 by 果果媽咪
[考研] 一志愿華東理工大學(xué)081700，初試分數(shù)271 +6	kotoko_ik 2026-03-23	7/350	2026-03-27 12:29 by 惠州彭于晏
[考研] 324求調(diào)劑 +5	hanamiko 2026-03-26	5/250	2026-03-27 10:33 by wangjy2002
[考研] 材料調(diào)劑 +8	匹克i 2026-03-23	8/400	2026-03-27 08:11 by hypershenger
[考研] 機械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3	@357 2026-03-25	3/150	2026-03-26 16:34 by by.MENG
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4	JKSOIID 2026-03-26	4/200	2026-03-26 15:58 by 不吃魚的貓
[考研] 263求調(diào)劑 +6	yqdszhdap－ 2026-03-22	10/500	2026-03-26 13:11 by 公瑾逍遙
[考研] 334分一志愿武理-080500 材料求調(diào)劑 +4	李李不服輸 2026-03-25	4/200	2026-03-25 21:26 by 星空星月
[考研] 【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

zhmdream

[求助] Matlab兩個小問題

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

sudo

sudo

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

zhmdream

sudo

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频