峰值演算法
峰值就是一堆數值中的最大值或者最小值吧,你可以使用max和min函數。
l=[1,4,5,6,74,23,2,1,5,7]
print max(l)
print min(l)
『貳』 快速查找二維數組的所有峰值,c語言實現最好,python也可以,最好能實現濾波。
#include<stdio.h>
#include<stdlib.h>
#include<limits.h>
#include<malloc.h>
#definemaxsize6//每個波形數據最大採集個數默認6
typedefstructarray1
{
int*data;
structarray1*next;
}ARR1;
ARR1*addNewArr(ARR1*arrHead,ARR1*arrTail);//插入一組波形數據節點返回尾節點
intfindMaxF(ARR1*arrHead);//查找最大峰值
intremoveMinF(ARR1*arrHead);//移除最小峰值數據組返回最小峰值
voidprintfArr(ARR1*arrHead);//列印數據鏈表
intmain()
{
ARR1*arrHead=(ARR1*)malloc(sizeof(ARR1));
arrHead->next=NULL;
ARR1*arrTail=NULL;
arrTail=addNewArr(arrHead,arrTail);//想插入幾組,就調用幾次我就測試6組
arrTail=addNewArr(arrHead,arrTail);
arrTail=addNewArr(arrHead,arrTail);
arrTail=addNewArr(arrHead,arrTail);
arrTail=addNewArr(arrHead,arrTail);
arrTail=addNewArr(arrHead,arrTail);
printf("採集的原數據組為:
");
printfArr(arrHead);
printf("最大峰值:%d,最小峰值值:%d
",findMaxF(arrHead),removeMinF(arrHead));
printf("刪除最小峰值數據組後的數據為:
");
printfArr(arrHead);
return0;
}
voidprintfArr(ARR1*arrHead)//列印數據鏈表
{
while(arrHead->next!=NULL)
{
printf("%d,%d,%d,%d,%d,%d
",arrHead->next->data[0],arrHead->next->data[1],arrHead->next->data[2],arrHead->next->data[3],arrHead->next->data[4],arrHead->next->data[5]);
arrHead=arrHead->next;
}
}
ARR1*addNewArr(ARR1*arrHead,ARR1*arrTail)//插入一組波形數據返回尾節點
{
int*data=(int*)malloc(sizeof(int)*maxsize),i;
printf("採集一組波形數據(最大採集個數%d):",maxsize);
for(i=0;i<maxsize;i++)
scanf("%d",&data[i]);
ARR1*arrNew=(ARR1*)malloc(sizeof(ARR1));
arrNew->data=data;
arrNew->next=NULL;
if(arrHead->next==NULL)//插入第一組數組作為首節點
arrHead->next=arrNew;
else
arrTail->next=arrNew;
arrTail=arrNew;
returnarrTail;
}
intfindMaxF(ARR1*arrHead)//查找最大峰值
{
inti,maxNum=0,maxf=0;
while(arrHead->next!=NULL)
{
maxNum=0;
for(i=0;i<maxsize;i++)//獲取每組數據的峰值
{
if(arrHead->next->data[i]>maxNum)
maxNum=arrHead->next->data[i];
}
if(maxNum>maxf)//獲取最大峰值
maxf=maxNum;
arrHead=arrHead->next;
}
returnmaxf;
}
intremoveMinF(ARR1*arrHead)//移除最小峰值數據組返回最小峰值
{
inti,maxNum=0,minf=INT_MAX;
ARR1*minDataSave=NULL,*arrHeadSave=arrHead;
while(arrHead->next!=NULL)
{
maxNum=0;
for(i=0;i<maxsize;i++)//獲取每組數據的峰值
{
if(arrHead->next->data[i]>maxNum)
maxNum=arrHead->next->data[i];
}
if(maxNum<minf)//獲取最小峰值
{
minDataSave=arrHead->next;
minf=maxNum;
}
arrHead=arrHead->next;
}
arrHead=arrHeadSave;
//移除最小峰值數據組
while(arrHead->next!=NULL)
{
if(arrHead->next==minDataSave)//刪除節點重組鏈表
{
arrHead->next=minDataSave->next;
minDataSave->next=NULL;
free(minDataSave->data);//釋放節點內存
free(minDataSave);
break;
}
arrHead=arrHead->next;
}
returnminf;
}
『叄』 matlab 求某函數峰值
首先,diff(y)算出來的z並不是嚴格意義上的y的微分,而是有一定誤差的,所以要它嚴格等於0除非是巧合才能碰到。
但是即使找絕對值z小於error(error是一個自己設定的比較小的數,作為誤差)的話也會有問題,主要是error的值不好確定,和x的取值有關,而且會出現如圖1所示的問題,有些點能探測到,有些又探測過多。
clearall;closeall;clc;
x=0:0.01:100;
y=x.*sin(x);
error=1e-2
z=diff(y);
length(z)
plot(x,y);holdon;
c=find(z<=error);
c=find(z(c)>=-error);
plot(x(c),y(c),'r*')
圖2
『肆』 給定一個二維數組.用Python設計實現演算法找到一個峰值
#include<stdio.h>
int main()
{int i,j,t,a[3][3]={{1,2,3},{4,5,6},{7,8,9}};
printf("原來的數組:\n");
for(i=0;i<3;i++)
{
for(j=0;j<3;j++)
printf("%d ",a[i][j]);
printf("\n");
}
for(i=0;i<2;i++)
for(j=i+1;j<3;j++)
{t=a[i][j];a[i][j]=a[j][i];a[j][i]=t;}
printf("轉置後的數組:\n");
for(i=0;i<3;i++)
{
『伍』 測量三極體正玄波電路示波器怎麼算測量頻率的峰值丶峰峰值
首先先先自矯,示波器本身輸出1KHZ信號,自矯完成會出現一個矩形波,橫2.5格,豎4格,2.5+2.5=5,4x5=2V(這是信號輸入電壓)
然後算輸入頻率,左邊的mvx右邊的mv=0.2x0.5=1vs,1處於1vs=1kHZ,1s=1000ms,1ms=1000vs
自矯部分算完了,接下來算輸入波形的峰值,封峰值:我們任意輸入一個正玄波,得出圖中的波形,假設mv檔調到0.2,ms檔調到0.5,先算周期的ms(毫秒):橫著數一個完整的正玄波為8格,8x0.5=4ms(這里算出一個周期4毫秒),然後用1處於4=0.25=250HZ(輸入頻率為250HZ,一個周期4毫秒)峰峰值演算法如下:從正玄波峰頂到峰底數為7格,7x0.2(ms取左邊的值,你打到多少就x多少)=1.4(V這是峰峰值),峰值就是峰峰值除於2,ok,教程到此結束!
有問題請用QQ搜索群雨道~
『陸』 電流電壓有效值,平均值
單相正弦交流電的有效值為220V、峰值為310V(有效值乘以1.414)、算術平均值為0V、整流平均值(絕對值的平均值)為198V(有效值除以1.1107)。
三相電峰值、算術平均值、整流平均值、峰值的演算法相同。
三相電的相電壓有效值為220V。
三相電的線電壓有效值為380V。
電流有效值並不等同於電流的平均值,交變電流的有效值是根據電流的熱效應定義的,交變電流的平均值是指在某段時間內平均電流的大小。這是兩個不同的物理量。
有效值在相同的電阻上分別通過直流電流和交流電流,經過一個交流周期的時間,如果它們在電阻上所消耗的電能相等的話,則把該直流電流(電壓)的大小作為交流電流(電壓)的有效值,正弦電流(電壓)的有效值等於其最大值(幅值)的1/√2,約0.707倍。
在正弦交流電流電中根據熱等效原理,定義電流和電壓的有效值為其瞬時值在一個周期內的均方根值。
有效值是根據電流熱效應來規定的,讓一個交流電流和一個直流電流分別通過阻值相同的電阻,如果在相同時間內產生的熱量相等,那麼就把這一直流電的數值叫做這一交流電的有效值。
在電工技術中,有時並不需要知道交流電的瞬時值,而規定一個能夠表徵其大小的特定值——有效值,其依據是交流電流和直流電流通過電阻時,電阻都要消耗電能(熱效應)。
『柒』 現在無錫電費如何演算法,什麼時間是峰值什麼時候是谷值
峰:0.56元(電費單顯示是四位小數,保留兩位是這么多)每度電,谷:0.36(也是保留兩位數)元每度電。
峰時間:早八點到晚九點;其他為谷
『捌』 求C#峰值演算法公式,在線等,急
能否說的詳細點...
這個峰值指什麼?
實時讀取波形計算峰值? 還是什麼....
『玖』 利用C8051F020實現音頻的峰值檢測 求峰值檢測軟體演算法
分析法
『拾』 python分治法求二維數組局部峰值方法
python分治法求二維數組局部峰值方法
下面小編就為大家分享一篇python分治法求二維數組局部峰值方法,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧
題目的意思大致是在一個n*m的二維數組中,找到一個局部峰值。峰值要求大於相鄰的四個元素(數組邊界以外視為負無窮),比如最後我們找到峰值A[j][i],則有A[j][i] > A[j+1][i] && A[j][i] > A[j-1][i] && A[j][i] > A[j][i+1] && A[j][i] > A[j][i-1]。返回該峰值的坐標和值。
當然,最簡單直接的方法就是遍歷所有數組元素,判斷是否為峰值,時間復雜度為O(n^2)
再優化一點求每一行(列)的最大值,再通過二分法找最大值列的峰值(具體方法可見一維數組求峰值),這種演算法時間復雜度為O(logn)
這里討論的是一種復雜度為O(n)的演算法,演算法思路分為以下幾步:
1、找「田」字。包括外圍的四條邊和中間橫豎兩條邊(圖中綠色部分),比較其大小,找到最大值的位置。(圖中的7)
2、找到田字中最大值後,判斷它是不是局部峰值,如果是返回該坐標,如果不是,記錄找到相鄰四個點中最大值坐標。通過該坐標所在的象限縮小范圍,繼續比較下一個田字
3、當范圍縮小到3*3時必定會找到局部峰值(也可能之前就找到了)
關於為什麼我們選擇的范圍內一定存在峰值,大家可以這樣想,首先我們有一個圈,我們已知有圈內至少有一個元素大於這個圈所有的元素,那麼,是不是這個圈中一定有一個最大值?
可能說得有點繞,但是多想想應該能夠理解,也可以用數學的反證法來證明。
演算法我們理解後接下來就是代碼實現了,這里我用的語言是python(初學python,可能有些用法上不夠簡潔請見諒),先上代碼:
import numpy as np
def max_sit(*n): #返回最大元素的位置
temp = 0
sit = 0
for i in range(len(n)):
if(n[i]>temp):
temp = n[i]
sit = i
return sit
def dp(s1,s2,e1,e2):
m1 = int((e1-s1)/2)+s1 #row
m2 = int((e2-s1)/2)+s2 #col
nub = e1-s1
temp = 0
sit_row = 0
sit_col = 0
for i in range(nub):
t = max_sit(list[s1][s2+i], #第一排
list[m1][s2+i], #中間排
list[e1][s2+i], #最後排
list[s1+i][s2], #第一列
list[s1+i][m2], #中間列
list[s1+i][e2], #最後列
temp)
if(t==6):
pass
elif(t==0):
temp = list[s1][s2+i]
sit_row = s1
sit_col = s2+i
elif(t==1):
temp = list[m1][s2+i]
sit_row = m1
sit_col = s2+i
elif(t==2):
temp = list[e1][s2+i]
sit_row = e1
sit_col = s2+i
elif(t==3):
temp = list[s1+i][s2]
sit_row = s1+i
sit_row = s2
elif(t==4):
temp = list[s1+i][m2]
sit_row = s1+i
sit_col = m2
elif(t==5):
temp = list[s1+i][e2]
sit_row = s1+i
sit_col = m2
t = max_sit(list[sit_row][sit_col], #中
list[sit_row-1][sit_col], #上
list[sit_row+1][sit_col], #下
list[sit_row][sit_col-1], #左
list[sit_row][sit_col+1]) #右
if(t==0):
return [sit_row-1,sit_col-1]
elif(t==1):
sit_row-=1
elif(t==2):
sit_row+=1
elif(t==3):
sit_col-=1
elif(t==4):
sit_col+=1
if(sit_row<m1):
e1 = m1
else:
s1 = m1
if(sit_col<m2):
e2 = m2
else:
s2 = m2
return dp(s1,s2,e1,e2)
f = open("demo.txt","r")
list = f.read()
list = list.split("n") #對行進行切片
list = ["0 "*len(list)]+list+["0 "*len(list)] #加上下的圍牆
for i in range(len(list)): #對列進行切片
list[i] = list[i].split()
list[i] = ["0"]+list[i]+["0"] #加左右的圍牆
list = np.array(list).astype(np.int32)
row_n = len(list)
col_n = len(list[0])
ans_sit = dp(0,0,row_n-1,col_n-1)
print("找到峰值點位於:",ans_sit)
print("該峰值點大小為:",list[ans_sit[0]+1,ans_sit[1]+1])
f.close()
首先我的輸入寫在txt文本文件里,通過字元串轉換變為二維數組,具體轉換過程可以看我上一篇博客——python中字元串轉換為二維數組。(需要注意的是如果在windows環境中split後的列表沒有空尾巴,所以不用加list.pop()這句話)。有的變動是我在二維數組四周加了「0」的圍牆。加圍牆可以再我們判斷峰值的時候不用考慮邊界問題。
max_sit(*n)函數用於找到多個值中最大值的位置,返回其位置,python的內構的max函數只能返回最大值,所以還是需要自己寫,*n表示不定長參數,因為我需要在比較田和十(判斷峰值)都用到這個函數
def max_sit(*n): #返回最大元素的位置
temp = 0
sit = 0
for i in range(len(n)):
if(n[i]>temp):
temp = n[i]
sit = i
return sit
dp(s1,s2,e1,e2)函數中四個參數的分別可看為startx,starty,endx,endy。即我們查找范圍左上角和右下角的坐標值。
m1,m2分別是row 和col的中間值,也就是田字的中間。
def dp(s1,s2,e1,e2):
m1 = int((e1-s1)/2)+s1 #row
m2 = int((e2-s1)/2)+s2 #col
依次比較3行3列中的值找到最大值,注意這里要求二維數組為正方形,如果為矩形需要做調整
for i in range(nub):
t = max_sit(list[s1][s2+i], #第一排
list[m1][s2+i], #中間排
list[e1][s2+i], #最後排
list[s1+i][s2], #第一列
list[s1+i][m2], #中間列
list[s1+i][e2], #最後列
temp)
if(t==6):
pass
elif(t==0):
temp = list[s1][s2+i]
sit_row = s1
sit_col = s2+i
elif(t==1):
temp = list[m1][s2+i]
sit_row = m1
sit_col = s2+i
elif(t==2):
temp = list[e1][s2+i]
sit_row = e1
sit_col = s2+i
elif(t==3):
temp = list[s1+i][s2]
sit_row = s1+i
sit_row = s2
elif(t==4):
temp = list[s1+i][m2]
sit_row = s1+i
sit_row = m2
elif(t==5):
temp = list[s1+i][e2]
sit_row = s1+i
sit_row = m2
判斷田字中最大值是不是峰值,並找不出相鄰最大值
t = max_sit(list[sit_row][sit_col], #中
list[sit_row-1][sit_col], #上
list[sit_row+1][sit_col], #下
list[sit_row][sit_col-1], #左
list[sit_row][sit_col+1]) #右
if(t==0):
return [sit_row-1,sit_col-1]
elif(t==1):
sit_row-=1
elif(t==2):
sit_row+=1
elif(t==3):
sit_col-=1
elif(t==4):
sit_col+=1
縮小范圍,遞歸求解
if(sit_row<m1):
e1 = m1
else:
s1 = m1
if(sit_col<m2):
e2 = m2
else:
s2 = m2
return dp(s1,s2,e1,e2)
好了,到這里代碼基本分析完了。如果還有不清楚的地方歡迎下方留言。
除了這種演算法外,我也寫一種貪心演算法來求解這道題,只可惜最壞的情況下演算法復雜度還是O(n^2),QAQ。
大體的思路就是從中間位置起找相鄰4個點中最大的點,繼續把該點來找相鄰最大點,最後一定會找到一個峰值點,有興趣的可以看一下,上代碼:
#!/usr/bin/python3
def dp(n):
temp = (str[n],str[n-9],str[n-1],str[n+1],str[n+9]) #中 上 左 右 下
sit = temp.index(max(temp))
if(sit==0):
return str[n]
elif(sit==1):
return dp(n-9)
elif(sit==2):
return dp(n-1)
elif(sit==3):
return dp(n+1)
else:
return dp(n+9)
f = open("/home/nancy/桌面/demo.txt","r")
list = f.read()
list = list.replace(" ","").split() #轉換為列表
row = len(list)
col = len(list[0])
str="0"*(col+3)
for x in list: #加圍牆 二維變一維
str+=x+"00"
str+="0"*(col+1)
mid = int(len(str)/2)
print(str,mid)
p = dp(mid)
print (p)
f.close()
以上這篇python分治法求二維數組局部峰值方法就是小編分享給大家的全部內容了,希望能給大家一個參考