声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1868|回复: 7

[经典算法] 求最值方法

[复制链接]
发表于 2008-8-14 14:23 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
y=f(x,t),如果我想求当y取得最大值时,所对应的x和t应该怎么弄啊?麻烦高手指教
回复
分享到:

使用道具 举报

发表于 2008-8-14 19:29 | 显示全部楼层

回复 楼主 kexin 的帖子

请问LZ想要用C解决此问题,还是其它的?
发表于 2008-8-15 08:59 | 显示全部楼层
找个合适的优化算法做一下就可以了,论坛有很多相关的实例
 楼主| 发表于 2008-8-15 14:10 | 显示全部楼层
对于优化算法,我不太懂这方面的知识。我想实现如图所示: 图1.jpg
请问怎么做呢,数学方法倒是容易算,可是编程如何实现呢?请求帮助!
发表于 2008-8-17 14:40 | 显示全部楼层
呵呵,挺简单的最值问题,充要条件如下:

(1)    f'(x0)=0;

(2)    f(x0)>=f(x0-);

(3)    f(x0)>=f(x0+);

即可。

具体程序中,你可以:

(1)首先求出f'(x)的解析式,定义为F(x);

(2)使用牛顿法,求出F(x)==0的解x0;

(3)定义条件语句,检验(2)得到的x0,并最终得出需要的解。

解毕

评分

1

查看全部评分

发表于 2008-8-24 10:00 | 显示全部楼层
按照LSD意见办理
 楼主| 发表于 2008-8-25 08:13 | 显示全部楼层
h=hanning(255);
[BT,tt,FreqBins]=tfrstft(y1,1:n,2048,h);
FreqBins = FreqBins(1:1024) * fs * 2;
Pyy=abs(BT).^2;
PP=Pyy(:);   %变成一维数组
[Pmax,Ploc]=max(PP);
fprintf('一维数组中最大值的位置: %10d\n',Ploc);
Tm=fix(Ploc/2048)+1;
Fm=mod(Ploc,2048);
fprintf('相应时间位置: %10d   频率位置: %10d\n',Tm,Fm);
fprintf('最大值的时间: %8.5f   频率: %8.5f\n',t(Tm),FreqBins(Fm));
我想求能量衰减到最大能量的1/3时,所对应的频率和时间,应该怎么弄啊?
 楼主| 发表于 2008-8-25 20:43 | 显示全部楼层
很开心,在自己的摸索下终于作出来了,呵呵分享一下
x=[1:500:2.5e+004];
x=x';
pp=max(x);
M=1/3*pp;
N=min(abs(x-M));
[m n]=find(x==max(M+-N));
Tm=fix(m/2048)+1;
Fm=mod(m,2048);
fprintf('相应时间位置: %10d   频率位置: %10d\n',Tm,Fm);

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-28 11:24 , Processed in 0.096217 second(s), 26 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表