右手握炮

新蟲 (小有名氣)

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

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

求教，我想用matlab求一組數(shù)據(jù)擬合曲線中斜率變化率最大的點，程序如下，出現(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ù)標量。

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復(fù)
一志愿華中科技大學，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學材料學碩321求調(diào)劑已經(jīng)有6人回復(fù)
085410人工智能專碩317求調(diào)劑（0854都可以）已經(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è)求點建議已經(jīng)有6人回復(fù)
求教各位有沒有免費使用 Matlab里的 Partial Differential Equation Toolbox的方法？已經(jīng)有3人回復(fù)
利用matlab求一階變系數(shù)微分方程組的數(shù)值解。已經(jīng)有7人回復(fù)
新人第一次發(fā)帖特此來問個matlab的數(shù)值積分的問題已經(jīng)有3人回復(fù)
matlab編熱泵模型，卡在壓縮機出口的制冷劑狀態(tài)這了已經(jīng)有13人回復(fù)
非常簡單��！MATLAB求解非線性一元方程已經(jīng)有5人回復(fù)
如何繪制向量場（自己編了Matlab程序但結(jié)果好像不正確，求高人指點一下）已經(jīng)有0人回復(fù)
linux下如何設(shè)置腳本文件讓matlab自動運行程序？已經(jīng)有4人回復(fù)
求助大神：要做參數(shù)回歸，fortran和matlab，該選擇哪個？已經(jīng)有14人回復(fù)
誰有成熟的求動力學參數(shù)的MATLAB的程序。已經(jīng)有11人回復(fù)
用matlab來實現(xiàn)分類匯總已經(jīng)有3人回復(fù)
matlab下bvp4c求解非線方程問題已經(jīng)有5人回復(fù)
matlab擬合擬均相動力學參數(shù) 已經(jīng)有6人回復(fù)

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

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

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

hsdwlxhx

金蟲 (正式寫手)

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

【答案】應(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è): 光學

【答案】應(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)驗: +831
數(shù)學EPI: 9
應(yīng)助: 791 (博后)
貴賓: 0.308
金幣: 24609
散金: 310
紅花: 75
帖子: 3004
在線: 881.4小時
蟲號: 1368482
注冊: 2011-08-14
性別: GG
專業(yè): 計算數(shù)學與科學工程計算
管轄: 數(shù)學

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

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

多項式求導命令用polyder函數(shù)!
例如:求x^2+2*x+1在[0.1：0.1：1]處的導數(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è): 燃燒學

引用回帖:

4樓: Originally posted by wurongjun at 2015-07-26 09:32:29
多項式求導命令用polyder函數(shù)!
例如:求x^2+2*x+1在處的導數(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)驗: +831
數(shù)學EPI: 9
應(yīng)助: 791 (博后)
貴賓: 0.308
金幣: 24609
散金: 310
紅花: 75
帖子: 3004
在線: 881.4小時
蟲號: 1368482
注冊: 2011-08-14
性別: GG
專業(yè): 計算數(shù)學與科學工程計算
管轄: 數(shù)學

引用回帖:

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è): 燃燒學

引用回帖:

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è): 燃燒學

引用回帖:

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ù)此樓

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

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

wurongjun

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

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

【答案】應(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è): 燃燒學

引用回帖:

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ā)表

[考研] 085700資源與環(huán)境308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-18 22:35 by bingxueer79
[考研] 求材料調(diào)劑 +8	隔壁陳先生 2026-03-12	8/400	2026-03-18 22:19 by li123456789.
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +5	小材化本科 2026-03-18	5/250	2026-03-18 21:53 by 無懈可擊111
[考研] 321求調(diào)劑 +3	何潤采123 2026-03-18	3/150	2026-03-18 21:27 by li123456789.
[考研] 一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑 +10	yangfz 2026-03-17	10/500	2026-03-18 20:14 by walc
[考研] 0703化學調(diào)劑 +7	妮妮ninicgb 2026-03-15	11/550	2026-03-18 19:48 by macy2011
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	9/450	2026-03-18 15:31 by vgtyfty
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 288求調(diào)劑，一志愿華南理工大學071005 +4	ioodiiij 2026-03-17	4/200	2026-03-18 12:36 by Linda Hu
[考研] 0703化學調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 070300化學學碩求調(diào)劑 +6	太想進步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 本科南京大學一志愿川大藥學327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 081200-11408-276學碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz