24小时热门版块排行榜    

查看: 5514  |  回复: 11

右手握炮

新虫 (小有名气)

[求助] 自己编了matlab程序有错误,求教问题出在哪里。。。 已有3人参与

求教,我想用matlab求一组数据拟合曲线中斜率变化率最大的点,程序如下,出现错误,不知道错在哪里。。。求指教。。。感激不尽。。。


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('原始数据','3阶曲线','斜率变化率');
错误使用 diff
差分阶数 N 必须为正整数标量。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

新的风暴已经出现 怎么能够停滞不前 穿越时空竭尽全力 我会来到你身边 微笑面对危险 梦想成真不会遥远 鼓起勇气坚定向前 奇迹一定会出现
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hsdwlxhx

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
上错误代码,我猜可能是第三行

[ 发自小木虫客户端 ]
BringLightonChip
2楼2015-07-26 09:22:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dalerxli

铁虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
同意楼上,第三行看不懂

[ 发自小木虫客户端 ]
3楼2015-07-26 09:25:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
多项式求导命令用polyder函数!
例如:求x^2+2*x+1在[0.1:0.1:1]处的导数值
>> 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
>>
善恶到头终有报,人间正道是沧桑.
4楼2015-07-26 09:32:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

右手握炮

新虫 (小有名气)

引用回帖:
4楼: Originally posted by wurongjun at 2015-07-26 09:32:29
多项式求导命令用polyder函数!
例如:求x^2+2*x+1在处的导数值
>> p=
p =
     1     2     1
>> dp=polyder(p)
dp =
     2     2
>> x=
x =
  Columns 1 through 5
    0.1000    0 ...

多谢,但却出现一下错误是咋回事呢?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('原始数据','3阶曲线','斜率变化率');

dp =

   -0.0000    0.0002   -0.1310

错误使用 plot
矢量长度必须相同。
新的风暴已经出现 怎么能够停滞不前 穿越时空竭尽全力 我会来到你身边 微笑面对危险 梦想成真不会遥远 鼓起勇气坚定向前 奇迹一定会出现
5楼2015-07-26 10:00:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

引用回帖:
5楼: Originally posted by 右手握炮 at 2015-07-26 10:00:30
多谢,但却出现一下错误是咋回事呢?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,'*')
善恶到头终有报,人间正道是沧桑.
6楼2015-07-26 10:05:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

右手握炮

新虫 (小有名气)

引用回帖:
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,'*')...

那请问该怎么改啊。。
新的风暴已经出现 怎么能够停滞不前 穿越时空竭尽全力 我会来到你身边 微笑面对危险 梦想成真不会遥远 鼓起勇气坚定向前 奇迹一定会出现
7楼2015-07-26 10:18:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

右手握炮

新虫 (小有名气)

引用回帖:
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,'*')...

怎样才能输出斜率值最大点的坐标?求指教。。。
新的风暴已经出现 怎么能够停滞不前 穿越时空竭尽全力 我会来到你身边 微笑面对危险 梦想成真不会遥远 鼓起勇气坚定向前 奇迹一定会出现
8楼2015-07-26 10:20:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
右手握炮: 金币+10, ★★★很有帮助, 等下,请帮人帮到底啊。。感激不尽啊。。 2015-07-26 10:44:59
用max找啊!
问题解决了也没有个谢谢呀!
缺乏回帖动力!
再见!
善恶到头终有报,人间正道是沧桑.
9楼2015-07-26 10:43:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

右手握炮

新虫 (小有名气)

引用回帖:
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('原始数据','3阶曲线','斜率变化率');
错误使用 plot
矢量长度必须相同。
新的风暴已经出现 怎么能够停滞不前 穿越时空竭尽全力 我会来到你身边 微笑面对危险 梦想成真不会遥远 鼓起勇气坚定向前 奇迹一定会出现
10楼2015-07-26 10:47:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 右手握炮 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工学调剂 +7 用户573181 2026-03-20 12/600 2026-03-24 11:17 by 用户573181
[考研] 306求0703调剂一志愿华中师范 +8 纸鱼ly 2026-03-21 9/450 2026-03-24 10:24 by 学术搬砖er
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗? +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 一志愿南航材料专317分求调剂 +4 炸呀炸呀炸薯条 2026-03-23 4/200 2026-03-24 07:32 by wangy0907
[考研] 环境学硕288求调剂 +8 皮皮皮123456 2026-03-22 8/400 2026-03-23 23:47 by 热情沙漠
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 070300化学求调剂 +8 苑豆豆 2026-03-20 8/400 2026-03-23 20:57 by baobaoye
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 初始318分求调剂(有工作经验) +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
信息提示
请填处理意见