zhitong.m

来自「数字信号处理课程的有关实验」· M 代码 · 共 42 行

val0=[];val=[];mean0=[];val1=[];valh=[];vallvbo=[];
for i=35:213;
    file=['06.2.23_1xshiyan\\' num2str(i)];
    fid=fopen(file,'rb');
    if(fid~=-1)
        val=fread(fid,20000,'uint8');
        fclose(fid);
    end       
    val0=[val0 val];%
    mean0(i)=mean(val);
end
figure(1),subplot(211),imshow(val0',[20,150]);
h0=fir1(64,0.33,'low');
for i=1:480
    xx=zeros(1,64);
    xx(13:64)=val0(i,1:52);
    ff=conv(xx,h0);
    ff=ff(33:97);
    gg(i,1:40)=ff(1:40);
end
for j=0:2
for i=1:480
    xx=zeros(64,1);
    xx(1:64)=val0(i,29+40*j:92+40*j);
    ff=conv(xx,h0);
    ff=ff(33:97);
    gg(i,41+40*j:80+40*j)=ff(1:40);
end
end
for i=1:480
    xx=zeros(1,64);
    xx(1:31)=val0(i,149:179);
    ff=conv(xx,h0);
    ff=ff(33:97);
    gg(i,149:179)=ff(1:31);
end
mean(mean(gg))
subplot(212),imshow(gg',[20,150]);