版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

lmszz

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 95
帖子: 1
在線: 1.4小時(shí)
蟲號(hào): 3192658
注冊(cè): 2014-05-09
性別: MM

[求助] 程序求助

幫忙看看matlab程序兩列隨機(jī)數(shù)據(jù)經(jīng)過互相關(guān)計(jì)算，為什么不能畫出正確的互相關(guān)圖并得到正確的延遲點(diǎn)數(shù)？新蟲存糧是12個(gè)金幣，不好意思，謝謝！
clear all;
Fs=1468;       %采樣頻率間序列
T=1/Fs;       %采樣周期
W=20479;       %信號(hào)長度(20479)
t=(0:W-1)*T; %時(shí)
A=xlsread('E:\0.xls','sheet1','A1:B2000');  %以矩陣形式讀進(jìn)數(shù)據(jù)
L=10;          %電極間距（可變）
x=1:1:2000;
N=2000;       %采樣信號(hào)長度
y1=A(:,1);
y2=A(:,2);
figure(1);
subplot(3,1,1),
plot(x,y1); %第一組數(shù)據(jù)畫圖
title('圖1');
subplot(3,1,2),
plot(x,y2); %第二組數(shù)據(jù)畫圖
title('圖2');
%設(shè)計(jì)低通濾波器
Wp =100/Fs;
Ws = 200/Fs; %截止頻率100Hz,阻帶截止頻率200Hz
[n,Wn] = buttord(Wp,Ws,1,50); %阻帶衰減大于50db,通帶紋波小于1db
%估算得到Butterworth低通濾波器的最小階數(shù)N(10)和3dB截止頻率Wn(0.0774)
[a,b]=butter(n,Wn);    %設(shè)計(jì)Butterworth低通濾波器，確定傳遞函數(shù)的分子分母系數(shù)
[h,f]=freqz(a,b,'whole',Fs);       %求數(shù)字低通濾波器的頻率響應(yīng)
f=(0:length(f)-1)'*Fs/length(f); %進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換
subplot(3,1,3),
plot(f,abs(h));    %繪制Butterworth低通濾波器的幅頻響應(yīng)圖
title('巴氏低通濾波器');
grid;
%低通濾波
p1=filter(a,b,y1);    %y1經(jīng)過低通濾波器以后的新數(shù)據(jù)
p2=filter(a,b,y2);    %y2經(jīng)過低通濾波器以后的新數(shù)據(jù)
figure(2);
subplot(2,1,1),
plot(x,p1);    %繪制y1經(jīng)過低通濾波器以后的時(shí)域圖形
title('除燥數(shù)據(jù)1');
subplot(2,1,2),
plot(x,p2);    %繪制y2經(jīng)過低通濾波器以后的時(shí)域圖形
title('除燥數(shù)據(jù)2');
figure(3);
plot(x,p1,x,p2);
%相關(guān)運(yùn)算
[a,b]=xcorr(p1,p2,'coeff');  %互相關(guān)計(jì)算, 返回矢量長度為2*N-1互相關(guān)函數(shù)序列,,當(dāng)b沒有指定時(shí)，其范圍為[-N+1,N-1]
dt=.1;
figure(4);
plot(b*dt,a);    %絕對(duì)互相關(guān)圖
title('互相關(guān)圖');
a=abs(a);
[amax,bmax]=max(a);
sample_delay=bmax-N;             %計(jì)算與中心點(diǎn)相差的樣點(diǎn)數(shù)
t0=sample_delay/Fs;
v=L/t0;

回復(fù)此樓

» 猜你喜歡

材料科學(xué)與工程求調(diào)劑已經(jīng)有8人回復(fù)
一志愿大連理工大學(xué)材料求調(diào)劑已經(jīng)有3人回復(fù)
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑已經(jīng)有9人回復(fù)
材料求調(diào)劑一志愿哈工大總分298分，前三科223分已經(jīng)有8人回復(fù)
274求調(diào)劑已經(jīng)有3人回復(fù)
合肥區(qū)域性重點(diǎn)一本招收調(diào)劑已經(jīng)有4人回復(fù)
085600材料與化工調(diào)劑已經(jīng)有10人回復(fù)
22408 359分調(diào)劑已經(jīng)有7人回復(fù)
抱歉已經(jīng)有4人回復(fù)
085602化工求調(diào)劑（331分）已經(jīng)有8人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2014-05-19 09:44:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 lmszz 的主題更新