butterworth matlab,Matlab实现Butterworth滤波器

butterworth matlab,Matlab实现Butterworth滤波器下面是用Matlab实现的Butterworth高通、低通滤波器。clc;clearall;closeall;I=imread(‘cameraman.tif’);subplot(3,2,1);imshow(I);title(‘原始图’);f=double(I);%数据类型转换,MATLAB不支持图像的无符号整型的计算g=fft2(f);%傅立叶变换g=fftshi…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46,售后保障稳定

下面是用Matlab实现的Butterworth高通、低通滤波器。

clc;clear all;close all;

I=imread(‘cameraman.tif’);

subplot(3,2,1);imshow(I);

title(‘原始图’);

f=double(I); % 数据类型转换,MATLAB不支持图像的无符号整型的计算

g=fft2(f); % 傅立叶变换

g=fftshift(g); % 转换数据矩阵

subplot(3,2,2);imshow(log(1+abs(g)),[]);

title(‘傅里叶频谱’);

%低通滤波***************************************************************

[M,N]=size(g);

nn=2; % 二阶巴特沃斯(Butterworth)低通滤波器

d0=10;

m=floor(M/2); n=floor(N/2);

for i=1:M

for j=1:N

d=sqrt((i-m)^2+(j-n)^2);

h=1/(1+0.414*(d/d0)^(2*nn)); % 计算低通滤波器传递函数

result(i,j)=h*g(i,j);

end

end

subplot(3,2,3);imshow(log(1+abs(result)),[]);

title(‘butterworth低通滤波’);

result=ifftshift(result);

J2=ifft2(result);

J3=uint8(real(J2));

subplot(3,2,4),imshow(J3); % 显示滤波处理后的图像

title(‘低通滤波图’);

%高通滤波

for i=1:M

for j=1:N

d=sqrt((i-m)^2+(j-n)^2);

h=1/(1+0.414*(d0/d)^(2*nn));% 计算传递函数

result(i,j)=h*g(i,j);

end

end

subplot(3,2,5);imshow(log(1+abs(result)),[]);

title(‘butterworth高通滤波’);

result=ifftshift(result);

J2=ifft2(result);

J3=uint8(real(J2));

subplot(3,2,6),imshow(J3); % 显示滤波处理后的图像

title(‘高通滤波图’);

效果图:

0818b9ca8b590ca3270a3433284dd417.png

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/204561.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • logback-spring.xml配置文件(最佳实践)

    logback-spring.xml配置文件(最佳实践)logback-spring.xml配置文件<?xmlversion="1.0"encoding="UTF-8"?><!–日志级别从低到高分为TRACE<DEBUG&l

    2025年6月30日
    1
  • javascript Date format(js日期格式化)

    javascript Date format(js日期格式化)这个很不错,好像是csdn的Meizz写的://对Date的扩展,将Date转化为指定格式的String//月(M)、日(d)、小时(h)、分(m)、秒(s)、季度(q)可以用1-2个占位符,//年(y)可以用1-4个占位符,毫秒(S)只能用1个占位符(是1-3位的数字)//例子://(newDate()).Format(“yyyy-MM-ddhh:mm:ss.S”)==>2006-07-0208:09:04.423//(newDate())

    2022年4月30日
    36
  • anaconda与pycharm的关系、安装「建议收藏」

    anaconda与pycharm的关系、安装「建议收藏」一、anaconda和pycharm的关系anaconda和miniconda的介绍Anaconda:是一个打包的集合,它里面预装好了conda、某个版本的python、众多packages、科学计算工具等等,就是把很多常用的不常用的库都给你装好了,并且在安装一个库时会询问是否安装与该库相关联的库,它将库之间的关系理的非常清楚,为用户提供极大的便利。Miniconda,顾名思义,它只包含最基本的内容——python与conda,以及相关的必须依赖项,对于空间要求严格的用户,Miniconda是一种

    2022年8月25日
    4
  • sql if语句实例[通俗易懂]

    sql if语句实例[通俗易懂]if函数语法:if(booleantestCondition,value1,value2)返回值:当判断条件testCondition为TRUE时,返回value1;否则返回value2示例:原始数据转化:SELECTapp_list,if(account_id>20,1,0)FROMtemp.bdl_op_crawler_loan0902…

    2022年7月17日
    12
  • 开发Windows Mobile今日插件 — 内存电量,桌面便笺,桌面记单词

    开发Windows Mobile今日插件 — 内存电量,桌面便笺,桌面记单词

    2021年7月31日
    78
  • pychame激活码(破解版激活)[通俗易懂]

    pychame激活码(破解版激活),https://javaforall.net/100143.html。详细ieda激活码不妨到全栈程序员必看教程网一起来了解一下吧!

    2022年3月16日
    55

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号