作业帮matlab 求某函数峰值
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/05/28 04:56:26
matlab 求某函数峰值
x = 0:0.1:100; y = x.*sin(x); plot(x,y); z=diff(y); c=find(z==0) x(c)
x = 0:0.1:100; y = x.*sin(x); plot(x,y); z=diff(y); c=find(z==0) x(c)
首先,diff(y)算出来的z并不是严格意义上的y的微分,而是有一定误差的,所以要它严格等于0除非是巧合才能碰到.但是即使找绝对值z小于error (error是一个自己设定的比较小的数,作为误差)的话也会有问题,主要是error的值不好确定,和x的取值有关,而且会出现如图1所示的问题,有些点能探测到,有些又探测过多.clear all;close all;clc;x = 0:0.01:100;y = x.*sin(x);error=1e-2z=diff(y);length(z)plot(x,y);hold on;c=find(z<=error);c=find(z(c)>=-error);plot(x(c),y(c),'r*')
图1
因此不建议用find函数,另一种方法是探测z的绝对值最小的点作为近似的峰值点,代码如下clear all;close all;clc;x = 0:0.01:100;y = x.*sin(x);z=diff(y);z=abs(z);j=1;for i=2:length(z)-1 if z(i)<z(i+1) && z(i)<z(i-1) % 如果微分的绝对值最小则记录下来 c(j)=i;j=j+1; endendfigure(1);hold on;plot(x,y);plot(x(c),y(c),'r*')效果很好,见图2.但不排除有更简洁的算法.
图2
图1
因此不建议用find函数,另一种方法是探测z的绝对值最小的点作为近似的峰值点,代码如下clear all;close all;clc;x = 0:0.01:100;y = x.*sin(x);z=diff(y);z=abs(z);j=1;for i=2:length(z)-1 if z(i)<z(i+1) && z(i)<z(i-1) % 如果微分的绝对值最小则记录下来 c(j)=i;j=j+1; endendfigure(1);hold on;plot(x,y);plot(x(c),y(c),'r*')效果很好,见图2.但不排除有更简洁的算法.
图2