matlab 求某函数峰值

matlab 求某函数峰值

x = 0:0.1:100; y = x.*sin(x); plot(x,y); z=diff(y); c=find(z==0) x(c)

问答/197℃/2024-10-25 01:53:18

优质解答：

首先,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

我来回答

猜你喜欢