右手握炮

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 4734.4
帖子: 119
在線: 179.8小時
蟲號: 2536380
注冊: 2013-07-07
性別: GG
專業(yè): 燃燒學(xué)

[求助] 自己編了matlab程序有錯誤，求教問題出在哪里。。。已有3人參與

求教，我想用matlab求一組數(shù)據(jù)擬合曲線中斜率變化率最大的點(diǎn)，程序如下，出現(xiàn)錯誤，不知道錯在哪里。。。求指教。。。感激不盡。。。

x=0:2:450;
data=xlsread('C:\Users\dell\Desktop\1.xlsx',1,'B2:B227');
y=[data];
y=y';
n=3;
p=polyfit(x,y,n);
dy=diff(y,x,2);
plot(x,y,'o',x,y,'b',x,dy,'*')
legend('原始數(shù)據(jù)','3階曲線','斜率變化率');
錯誤使用 diff
差分階數(shù) N 必須為正整數(shù)標(biāo)量。

回復(fù)此樓

» 猜你喜歡

321求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中海洋材料工程專碩330分求調(diào)劑已經(jīng)有7人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有5人回復(fù)
288求調(diào)劑已經(jīng)有3人回復(fù)
材料080500調(diào)劑求收留已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有9人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有8人回復(fù)
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑已經(jīng)有10人回復(fù)
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑已經(jīng)有3人回復(fù)
332求調(diào)劑已經(jīng)有3人回復(fù)

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

用Matlab編程序來解決一個小問題，但結(jié)果不太滿意，哪位蟲友大神能給些建議哈已經(jīng)有6人回復(fù)
自己編了個matlab數(shù)據(jù)擬合公式，擬合成功但不顯示參數(shù)，請問問題出在哪了？已經(jīng)有11人回復(fù)
大家是如何管理自己的Matlab仿真程序的呢？已經(jīng)有8人回復(fù)
圖像分割的超像素SLIC算法有沒有MATLAB語言的程序已經(jīng)有10人回復(fù)
MATLAB自編工具箱已經(jīng)有0人回復(fù)
請用matlab編下代碼，SO EASY只有一個未知數(shù) 已經(jīng)有4人回復(fù)
自己的寫的matlab跑不動，哪位大神幫我用牛頓迭代編一下已經(jīng)有4人回復(fù)
matlab小作業(yè)求點(diǎn)建議已經(jīng)有6人回復(fù)
求教各位有沒有免費(fèi)使用 Matlab里的 Partial Differential Equation Toolbox的方法？已經(jīng)有3人回復(fù)
利用matlab求一階變系數(shù)微分方程組的數(shù)值解。已經(jīng)有7人回復(fù)
新人第一次發(fā)帖特此來問個matlab的數(shù)值積分的問題已經(jīng)有3人回復(fù)
matlab編熱泵模型，卡在壓縮機(jī)出口的制冷劑狀態(tài)這了已經(jīng)有13人回復(fù)
非常簡單��！MATLAB求解非線性一元方程已經(jīng)有5人回復(fù)
如何繪制向量場（自己編了Matlab程序但結(jié)果好像不正確，求高人指點(diǎn)一下）已經(jīng)有0人回復(fù)
linux下如何設(shè)置腳本文件讓matlab自動運(yùn)行程序？已經(jīng)有4人回復(fù)
求助大神：要做參數(shù)回歸，fortran和matlab，該選擇哪個？已經(jīng)有14人回復(fù)
誰有成熟的求動力學(xué)參數(shù)的MATLAB的程序。已經(jīng)有11人回復(fù)
用matlab來實(shí)現(xiàn)分類匯總已經(jīng)有3人回復(fù)
matlab下bvp4c求解非線方程問題已經(jīng)有5人回復(fù)
matlab擬合擬均相動力學(xué)參數(shù) 已經(jīng)有6人回復(fù)

新的風(fēng)暴已經(jīng)出現(xiàn) 怎么能夠停滯不前穿越時空竭盡全力我會來到你身邊微笑面對危險夢想成真不會遙遠(yuǎn) 鼓起勇氣堅定向前奇跡一定會出現(xiàn)

1樓 2015-07-26 09:09:47

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

hsdwlxhx

金蟲 (正式寫手)

應(yīng)助: 20 (小學(xué)生)
金幣: 2986.8
散金: 204
紅花: 6
帖子: 430
在線: 318.5小時
蟲號: 1537165
注冊: 2011-12-14
性別: GG
專業(yè): 配位化學(xué)

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

感謝參與，應(yīng)助指數(shù) +1

上錯誤代碼，我猜可能是第三行

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

BringLightonChip

2樓2015-07-26 09:22:12

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

dalerxli

鐵蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 2413
帖子: 242
在線: 176.3小時
蟲號: 1864211
注冊: 2012-06-19
性別: GG
專業(yè): 光學(xué)

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

感謝參與，應(yīng)助指數(shù) +1

同意樓上，第三行看不懂

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

3樓2015-07-26 09:25:04

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

wurongjun

專家顧問 (職業(yè)作家)

專家經(jīng)驗(yàn): +831
數(shù)學(xué)EPI: 9
應(yīng)助: 791 (博后)
貴賓: 0.308
金幣: 24609
散金: 310
紅花: 75
帖子: 3004
在線: 881.4小時
蟲號: 1368482
注冊: 2011-08-14
性別: GG
專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算
管轄: 數(shù)學(xué)

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

感謝參與，應(yīng)助指數(shù) +1

多項式求導(dǎo)命令用polyder函數(shù)!
例如:求x^2+2*x+1在[0.1：0.1：1]處的導(dǎo)數(shù)值
>> p=[1 2 1]
p =
   1    2    1
>> dp=polyder(p)
dp =
   2    2
>> x=[0.1 :0.1:1]
x =
  Columns 1 through 5
0.1000 0.2000 0.3000 0.4000 0.5000
  Columns 6 through 10
0.6000 0.7000 0.8000 0.9000 1.0000
>> polyval(dp,x)
ans =
  Columns 1 through 5
2.2000 2.4000 2.6000 2.8000 3.0000
  Columns 6 through 10
3.2000 3.4000 3.6000 3.8000 4.0000
>>

贊一下

回復(fù)此樓

善惡到頭終有報,人間正道是滄桑.

4樓2015-07-26 09:32:29

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

右手握炮

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 4734.4
帖子: 119
在線: 179.8小時
蟲號: 2536380
注冊: 2013-07-07
性別: GG
專業(yè): 燃燒學(xué)

引用回帖:

4樓: Originally posted by wurongjun at 2015-07-26 09:32:29
多項式求導(dǎo)命令用polyder函數(shù)!
例如:求x^2+2*x+1在處的導(dǎo)數(shù)值
>> p=
p =
   1    2    1
>> dp=polyder(p)
dp =
   2    2
>> x=
x =
  Columns 1 through 5
0.1000 0 ...

多謝，但卻出現(xiàn)一下錯誤是咋回事呢？x=0:2:450;
data=xlsread('C:\Users\dell\Desktop\1.xlsx',1,'B2:B227');
y=[data];
y=y';
n=3;
p=polyfit(x,y,n);
dp=polyder(p)
plot(x,y,'o',x,y,'b',x,dp,'*')
legend('原始數(shù)據(jù)','3階曲線','斜率變化率');

dp =

-0.0000 0.0002 -0.1310

錯誤使用 plot
矢量長度必須相同。

贊一下

回復(fù)此樓

5樓2015-07-26 10:00:30

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

wurongjun

專家顧問 (職業(yè)作家)

專家經(jīng)驗(yàn): +831
數(shù)學(xué)EPI: 9
應(yīng)助: 791 (博后)
貴賓: 0.308
金幣: 24609
散金: 310
紅花: 75
帖子: 3004
在線: 881.4小時
蟲號: 1368482
注冊: 2011-08-14
性別: GG
專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算
管轄: 數(shù)學(xué)

引用回帖:

5樓: Originally posted by 右手握炮 at 2015-07-26 10:00:30
多謝，但卻出現(xiàn)一下錯誤是咋回事呢？x=0:2:450;
data=xlsread('C:\Users\dell\Desktop\1.xlsx',1,'B2:B227');
y=;
y=y';
n=3;
p=polyfit(x,y,n);
dp=polyder(p)
plot(x,y,'o',x,y,'b',x,dp,'*')
legend('原 ...

你沒有代值啊!
dy=polyval(dp,x);
plot(x,y,'o',x,y,'b',x,dy,'*')

贊一下

回復(fù)此樓

善惡到頭終有報,人間正道是滄桑.

6樓2015-07-26 10:05:12

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

右手握炮

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 4734.4
帖子: 119
在線: 179.8小時
蟲號: 2536380
注冊: 2013-07-07
性別: GG
專業(yè): 燃燒學(xué)

引用回帖:

6樓: Originally posted by wurongjun at 2015-07-26 10:05:12
你沒有代值啊!
dy=polyval(dp,x);
plot(x,y,'o',x,y,'b',x,dy,'*')...

那請問該怎么改啊。。

贊一下

回復(fù)此樓

7樓2015-07-26 10:18:17

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

右手握炮

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 4734.4
帖子: 119
在線: 179.8小時
蟲號: 2536380
注冊: 2013-07-07
性別: GG
專業(yè): 燃燒學(xué)

引用回帖:

6樓: Originally posted by wurongjun at 2015-07-26 10:05:12
你沒有代值啊!
dy=polyval(dp,x);
plot(x,y,'o',x,y,'b',x,dy,'*')...

怎樣才能輸出斜率值最大點(diǎn)的坐標(biāo)？求指教。。。

贊一下

回復(fù)此樓

8樓2015-07-26 10:20:08

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

wurongjun

專家顧問 (職業(yè)作家)

專家經(jīng)驗(yàn): +831
數(shù)學(xué)EPI: 9
應(yīng)助: 791 (博后)
貴賓: 0.308
金幣: 24609
散金: 310
紅花: 75
帖子: 3004
在線: 881.4小時
蟲號: 1368482
注冊: 2011-08-14
性別: GG
專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算
管轄: 數(shù)學(xué)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
右手握炮: 金幣+10, ★★★很有幫助, 等下，請幫人幫到底啊。。感激不盡啊。。 2015-07-26 10:44:59

用max找啊!
問題解決了也沒有個謝謝呀!
缺乏回帖動力!
再見!

贊一下

回復(fù)此樓

善惡到頭終有報,人間正道是滄桑.

9樓2015-07-26 10:43:27

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

右手握炮

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 4734.4
帖子: 119
在線: 179.8小時
蟲號: 2536380
注冊: 2013-07-07
性別: GG
專業(yè): 燃燒學(xué)

引用回帖:

9樓: Originally posted by wurongjun at 2015-07-26 10:43:27
用max找啊!
問題解決了也沒有個謝謝呀!
缺乏回帖動力!
再見!

哥們。。怎么賦值之后還是有錯。。。求指教啊。。只要解決問題金幣可以追加啊。。。謝謝拉
x=0:2:450;
data=xlsread('C:\Users\dell\Desktop\1.xlsx',1,'B2:B227');
y=[data];
y=y';
n=3;
p=polyfit(x,y,n);
dp=polyder(p);
x=[0:2:450];
plot(x,y,'o',x,y,'b',x,dp,'*')
legend('原始數(shù)據(jù)','3階曲線','斜率變化率');
錯誤使用 plot
矢量長度必須相同。

贊一下

回復(fù)此樓

10樓2015-07-26 10:47:54

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

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

[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +10	SUICHILD 2026-03-12	10/500	2026-03-19 09:51 by Delta2012
[教師之家] 焦慮 +9	水冰月月野兔 2026-03-13	13/650	2026-03-19 09:50 by otani
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無懈可擊111
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 085600材料與化工求調(diào)劑 +5	緒幸與子 2026-03-17	5/250	2026-03-17 16:40 by laoshidan
[考研] 274求調(diào)劑 +5	時間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎覬?/a>
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz