怎样给左图添加标题size_fft1=247; size_fft2=545; S = vortexPatternnew(size_fft1,size_fft2); subplot(121) imagesc(real(fftshift(S))) colormap(jet); colorbar; axis off; caxis([-1,1]); subplot(122) imagesc(imag(fftshift(S))) colormap(jet); colorbar; axis off; caxis([-1,1]); set(gcf,'color','w'); set(gcf, 'Position', [100 200 1200 400]); F_h = I_AC_fft.*S; F_inv = ifft2(F_h); RAC=abs(I_AC+1j.*abs(F_inv(1:size(I11,1),1:size(I11,2)))); imagesc(abs(F_inv(1:size(I11,1),1:size(I11,2))));colorbar; caxis([0,0.8]);colormap(jet);axis off;title('去直流图像Hilbert变换','Fontsize',15); set(gcf,'color','w'); set(gcf, 'Position', [100 200 600 400]);

时间: 2024-04-10 15:35:04 浏览: 231
要给左图添加标题,你可以使用`title`函数在`subplot(121)`之前添加以下代码行: ``` title('实部','Fontsize',15); ``` 同样地,在右图上添加标题,你可以使用`title`函数在`subplot(122)`之前添加以下代码行: ``` title('虚部','Fontsize',15); ``` 这将在左图和右图上分别添加标题为"实部"和"虚部",标题字体大小为15。 请注意,你还需要在第一个`figure`之前添加以下代码行,以确保正确显示图像和标题: ``` clf; ``` 这将清除当前图形窗口并重置图形状态。 最终修改后的代码如下: ```matlab size_fft1=247; size_fft2=545; S = vortexPatternnew(size_fft1,size_fft2); clf; subplot(121) imagesc(real(fftshift(S))) colormap(jet); colorbar; axis off; caxis([-1,1]); title('实部','Fontsize',15); subplot(122) imagesc(imag(fftshift(S))) colormap(jet); colorbar; axis off; caxis([-1,1]); title('虚部','Fontsize',15); set(gcf,'color','w'); set(gcf, 'Position', [100 200 1200 400]); F_h = I_AC_fft.*S; F_inv = ifft2(F_h); RAC=abs(I_AC+1j.*abs(F_inv(1:size(I11,1),1:size(I11,2)))); imagesc(abs(F_inv(1:size(I11,1),1:size(I11,2)))); colorbar; caxis([0,0.8]); colormap(jet); axis off; title('去直流图像Hilbert变换','Fontsize',15); set(gcf,'color','w'); set(gcf, 'Position', [100 200 600 400]); ```
阅读全文

相关推荐

import numpy as npimport cv2# 读取图像img = cv2.imread('lena.png', 0)# 添加高斯噪声mean = 0var = 0.1sigma = var ** 0.5noise = np.random.normal(mean, sigma, img.shape)noisy_img = img + noise# 定义维纳滤波器函数def wiener_filter(img, psf, K=0.01): # 计算傅里叶变换 img_fft = np.fft.fft2(img) psf_fft = np.fft.fft2(psf) # 计算功率谱 img_power = np.abs(img_fft) ** 2 psf_power = np.abs(psf_fft) ** 2 # 计算信噪比 snr = img_power / (psf_power + K) # 计算滤波器 result_fft = img_fft * snr / psf_fft result = np.fft.ifft2(result_fft) # 返回滤波结果 return np.abs(result)# 定义维纳滤波器的卷积核kernel_size = 3kernel = np.ones((kernel_size, kernel_size)) / kernel_size ** 2# 计算图像的自相关函数acf = cv2.calcHist([img], [0], None, [256], [0, 256])# 计算维纳滤波器的卷积核gamma = 0.1alpha = 0.5beta = 1 - alpha - gammapsf = np.zeros((kernel_size, kernel_size))for i in range(kernel_size): for j in range(kernel_size): i_shift = i - kernel_size // 2 j_shift = j - kernel_size // 2 psf[i, j] = np.exp(-np.pi * ((i_shift ** 2 + j_shift ** 2) / (2 * alpha ** 2))) * np.cos(2 * np.pi * (i_shift + j_shift) / (2 * beta))psf = psf / np.sum(psf)# 对带噪声图像进行维纳滤波filtered_img = wiener_filter(noisy_img, psf)# 显示结果cv2.imshow('Original Image', img)cv2.imshow('Noisy Image', noisy_img)cv2.imshow('Filtered Image', filtered_img)cv2.waitKey(0)cv2.destroyAllWindows()这段代码报错为Traceback (most recent call last): File "<input>", line 1, in <module> File "D:\Pycharm\PyCharm 2020.3.5\plugins\python\helpers\pydev\_pydev_bundle\pydev_umd.py", line 197, in runfile pydev_imports.execfile(filename, global_vars, local_vars) # execute the script File "D:\Pycharm\PyCharm 2020.3.5\plugins\python\helpers\pydev\_pydev_imps\_pydev_execfile.py", line 18, in execfile exec(compile(contents+"\n", file, 'exec'), glob, loc) File "E:/Python_project/class_/weinalvboqi.py", line 54, in <module> filtered_img = wiener_filter(noisy_img, psf) File "E:/Python_project/class_/weinalvboqi.py", line 25, in wiener_filter snr = img_power / (psf_power + K) ValueError: operands could not be broadcast together with shapes (1024,2800) (3,3)什么意思,如何修改

% 定义一些常量fft_size = 2048;hop_size = fft_size/4;min_freq = 80;max_freq = 1000;% 读取音频文件filename = 'example.aac';[x, Fs] = audioread(filename);% 计算音高[f0, ~] = yin(x, Fs, fft_size, hop_size, min_freq, max_freq);f0 = medfilt1(f0, 5); % 中值滤波midi = freq2midi(f0);% 计算主音调[~, max_idx] = max(histcounts(midi, 1:128));dominant_note = max_idx - 1;% 输出结果fprintf('主音调:%.2f Hz\n', midi2freq(dominant_note));function [f0, rms_energy] = yin(x, fs, fft_size, hop_size, min_freq, max_freq)% YIN算法计算音频信号的基频% 初始化变量n_frames = floor((length(x)-fft_size)/hop_size) + 1;f0 = zeros(n_frames, 1);rms_energy = zeros(n_frames, 1);% 计算自相关函数x = x(:);x = [x; zeros(fft_size, 1)];acf = xcorr(x, fft_size, 'coeff');acf = acf(ceil(length(acf)/2):end);% 计算差值函数d = zeros(fft_size, n_frames);for i = 1:n_frames frame = x((i-1)*hop_size+1:(i-1)*hop_size+fft_size); for tau = 1:fft_size d(tau, i) = sum((frame(1:end-tau) - frame(1+tau:end)).^2); endend% 计算自相关函数的倒数acf_recip = acf(end:-1:1);acf_recip(1) = acf_recip(2);acf_recip = acf_recip.^(-1);% 计算累积平均能量running_sum = 0;for i = 1:n_frames running_sum = running_sum + sum(x((i-1)*hop_size+1:(i-1)*hop_size+fft_size).^2); rms_energy(i) = sqrt(running_sum/fft_size);end% 计算基频for i = 1:n_frames r = acf_recip./(acf(i:end).*acf_recip(1:end-i+1)); r(1:i) = 0; r(max_freq/fs*fft_size+1:end) = 0; [~, j] = min(r(min_freq/fs*fft_size+1:max_freq/fs*fft_size)); f0(i) = fs/j;endendfunction midi = freq2midi(freq)% 将频率转换为MIDI码midi = 12*log2(freq/440) + 69;endfunction freq = midi2freq(midi)% 将MIDI码转换为频率freq = 440*2^((midi-69)/12);end对于此运算,数组的大小不兼容。 出错 yin (第 26 行) r = acf_recip./(acf(i:end).*acf_recip(1:end-i+1)); 出错 untitled2 (第 8 行) [f0, ~] = yin(x, Fs, fft_size, hop_size, min_freq, max_freq);请修改以上错误

大家在看

recommend-type

计算机图形学-小型图形绘制程序

计算机图形学-小型图形绘制程序
recommend-type

安装验证-浅谈mysql和mariadb区别

3.5 安装验证 客户机上能够启动软件就说明安装成功。 MotorSolve 成功画面 3.6 帮助 MotorSolve 上端的界面中的帮助按钮,点击可以查看详细的说明
recommend-type

基于Python深度学习的目标跟踪系统的设计与实现+全部资料齐全+部署文档.zip

【资源说明】 基于Python深度学习的目标跟踪系统的设计与实现+全部资料齐全+部署文档.zip基于Python深度学习的目标跟踪系统的设计与实现+全部资料齐全+部署文档.zip 【备注】 1、该项目是个人高分项目源码,已获导师指导认可通过,答辩评审分达到95分 2、该资源内项目代码都经过测试运行成功,功能ok的情况下才上传的,请放心下载使用! 3、本项目适合计算机相关专业(人工智能、通信工程、自动化、电子信息、物联网等)的在校学生、老师或者企业员工下载使用,也可作为毕业设计、课程设计、作业、项目初期立项演示等,当然也适合小白学习进阶。 4、如果基础还行,可以在此代码基础上进行修改,以实现其他功能,也可直接用于毕设、课设、作业等。 欢迎下载,沟通交流,互相学习,共同进步!
recommend-type

国密SM4加解密SM2签名验签for delphi等语言.rar

基于C#编写的COM组件DLL,可实现SM2签名验签,SM4加解密,100%适用于黑龙江省国家医保接口中进行应用。 1、调用DLL名称:JQSM2SM4.dll 加解密类名:JQSM2SM4.SM2SM4Util CLSID=5B38DCB3-038C-4992-9FA3-1D697474FC70 2、GetSM2SM4函数说明 函数原型public string GetSM2SM4(string smType, string sM2Prikey, string sM4Key, string sInput) 1)参数一smType:填写固定字符串,识别功能,分别实现SM2签名、SM4解密、SM4加密。SM2签名入参填写“SM2Sign”、SM4解密入参填写“SM4DecryptECB”、SM4加密入参填写“SM4EncryptECB”. 2)参数二sM2Prikey:SM2私钥 3)参数三sM4Key:SM4密钥 4)参数四sInput:当smType=SM2Sign,则sInput入参填写SM4加密串;当smType=SM4DecryptECB,则sInput入参填写待解密SM4密文串;当smType=SM4EncryptECB,则sInput入参填写待加密的明文串; 5)函数返回值:当smType=SM2Sign,则返回SM2签名信息;当smType=SM4DecryptECB,则返回SM4解密信息;当smType=SM4EncryptECB,则返回SM4加密信息;异常时,则返回“加解密异常:详细错误说明” 3、购买下载后,可加QQ65635204、微信feisng,免费提供技术支持。 4、注意事项: 1)基于.NET框架4.0编写,常规win7、win10一般系统都自带无需安装,XP系统则需安装;安装包详见压缩包dotNetFx40_Full_x86_x64.exe 2)C#编写的DLL,需要注册,解压后放入所需位置,使用管理员权限运行“JQSM2SM4注册COM.bat”即可注册成功,然后即可提供给第三方软件进行使用,如delphi等。
recommend-type

基于Android Studio开发的安卓的通讯录管理app

功能包含:新增联系人、编辑联系人、删除联系人、拨打电话、发送短信等相关操作。 资源包含源码:1、apk安装包 2、演示视频 3、 基本安装环境、4、运行文档 5、以及源代码

最新推荐

recommend-type

CSV_MATLAB_FFT

CSV_MATLAB_FFT 该资源主要介绍了如何将示波器波形CSV数据导入到MATLAB,并利用FFTAnalysis进行快速傅里叶变换(FFT)。下面是详细的知识点解释: 1. CSV数据导入MATLAB 将示波器波形CSV数据导入到MATLAB,可以...
recommend-type

(电磁仿真必备)maxwell_电机气隙磁密与用matlab进行fft谐波分析.doc

电磁仿真设计之电机气隙磁密与 MATLAB 进行 FFT 谐波分析 一、电机气隙磁密分析 电磁仿真设计是电机设计中不可或缺的一部分,而电机气隙磁密的分析是电磁仿真设计的重要内容之一。在本文中,我们将使用 Maxwell ...
recommend-type

Matlab的FFT算法程序-MATLAB_FFT.doc

1. **基2时间抽取FFT算法**: 时间抽取FFT,也称为“Radix-2 Decimation in Time”算法,其主要思想是将序列分段,通过级联一系列蝶形运算(Butterfly Operations)来完成DFT。代码中的关键步骤如下: - `M=...
recommend-type

基于OpenCV的人脸识别小程序.zip

【项目资源】: 包含前端、后端、移动开发、操作系统、人工智能、物联网、信息化管理、数据库、硬件开发、大数据、课程资源、音视频、网站开发等各种技术项目的源码。 包括STM32、ESP8266、PHP、QT、Linux、iOS、C++、Java、python、web、C#、EDA、proteus、RTOS等项目的源码。 【项目质量】: 所有源码都经过严格测试,可以直接运行。 功能在确认正常工作后才上传。 【适用人群】: 适用于希望学习不同技术领域的小白或进阶学习者。 可作为毕设项目、课程设计、大作业、工程实训或初期项目立项。 【附加价值】: 项目具有较高的学习借鉴价值,也可直接拿来修改复刻。 对于有一定基础或热衷于研究的人来说,可以在这些基础代码上进行修改和扩展,实现其他功能。 【沟通交流】: 有任何使用上的问题,欢迎随时与博主沟通,博主会及时解答。 鼓励下载和使用,并欢迎大家互相学习,共同进步。。内容来源于网络分享,如有侵权请联系我删除。另外如果没有积分的同学需要下载,请私信我。
recommend-type

免安装JDK 1.8.0_241:即刻配置环境运行

资源摘要信息:"JDK 1.8.0_241 是Java开发工具包(Java Development Kit)的版本号,代表了Java软件开发环境的一个特定发布。它由甲骨文公司(Oracle Corporation)维护,是Java SE(Java Platform, Standard Edition)的一部分,主要用于开发和部署桌面、服务器以及嵌入式环境中的Java应用程序。本版本是JDK 1.8的更新版本,其中的241代表在该版本系列中的具体更新编号。此版本附带了Java源码,方便开发者查看和学习Java内部实现机制。由于是免安装版本,因此不需要复杂的安装过程,解压缩即可使用。用户配置好环境变量之后,即可以开始运行和开发Java程序。" 知识点详细说明: 1. JDK(Java Development Kit):JDK是进行Java编程和开发时所必需的一组工具集合。它包含了Java运行时环境(JRE)、编译器(javac)、调试器以及其他工具,如Java文档生成器(javadoc)和打包工具(jar)。JDK允许开发者创建Java应用程序、小程序以及可以部署在任何平台上的Java组件。 2. Java SE(Java Platform, Standard Edition):Java SE是Java平台的标准版本,它定义了Java编程语言的核心功能和库。Java SE是构建Java EE(企业版)和Java ME(微型版)的基础。Java SE提供了多种Java类库和API,包括集合框架、Java虚拟机(JVM)、网络编程、多线程、IO、数据库连接(JDBC)等。 3. 免安装版:通常情况下,JDK需要进行安装才能使用。但免安装版JDK仅需要解压缩到磁盘上的某个目录,不需要进行安装程序中的任何步骤。用户只需要配置好环境变量(主要是PATH、JAVA_HOME等),就可以直接使用命令行工具来运行Java程序或编译代码。 4. 源码:在软件开发领域,源码指的是程序的原始代码,它是由程序员编写的可读文本,通常是高级编程语言如Java、C++等的代码。本压缩包附带的源码允许开发者阅读和研究Java类库是如何实现的,有助于深入理解Java语言的内部工作原理。源码对于学习、调试和扩展Java平台是非常有价值的资源。 5. 环境变量配置:环境变量是操作系统中用于控制程序执行环境的参数。在JDK中,常见的环境变量包括JAVA_HOME和PATH。JAVA_HOME是JDK安装目录的路径,配置此变量可以让操作系统识别到JDK的位置。PATH变量则用于指定系统命令查找的路径,将JDK的bin目录添加到PATH后,就可以在命令行中的任何目录下执行JDK中的命令,如javac和java。 在实际开发中,了解并正确配置JDK对于Java开发者来说是一个基础且重要的环节。掌握如何安装和配置JDK,以及如何理解JDK中的源码和各种工具,对于进行Java编程和解决问题至关重要。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

【提升效率与稳定性】:深入掌握单相整流器的控制策略

# 摘要 单相整流器作为电力电子技术的基础设备,广泛应用于电力系统和工业自动化领域。本文首先概述了单相整流器的基本概念,详细分析了其工作原理,包括电力电子器件的角色、数学模型以及控制策略。随后,针对控制技术进行了深入探讨,包括传统控制技术如PID控制和高级控制策略如模糊逻辑控制及自适应控制的应用。进一步地,本文提出了一系列提升单相整流器性能的方法,如效率提升策略、系统稳定性措施,以及性能评估与优化方案。最后,文章探讨了单相整流器在不同领域的应用案例,并对新兴技术与市场趋势进行了展望。 # 关键字 单相整流器;电力电子器件;数学模型;控制策略;性能提升;应用案例 参考资源链接:[单相PWM整
recommend-type

你看这是ashx映射的cs文件初始代码,你看这里边根本就没有写对action参数进行任何操作但你.ashx?action=submit这样去做他就能返回出数据这是为什么

在ASP.NET中,.ashx文件(也称为HTTP处理程序)用于处理HTTP请求并返回响应。即使在初始代码中没有对action参数进行任何操作,系统仍然可以通过默认的ProcessRequest方法处理请求并返回数据。 当你在URL中传递参数(如?action=submit)时,这些参数会被包含在请求的查询字符串中。虽然你的代码没有显式地处理这些参数,但默认的ProcessRequest方法会接收这些参数并执行一些默认操作。 以下是一个简单的.ashx文件示例: ```csharp <%@ WebHandler Language="C#" Class="MyHandler" %> us
recommend-type

机器学习预测葡萄酒评分:二值化品尝笔记的应用

资源摘要信息:"wine_reviewer:使用机器学习基于二值化的品尝笔记来预测葡萄酒评论分数" 在当今这个信息爆炸的时代,机器学习技术已经被广泛地应用于各个领域,其中包括食品和饮料行业的质量评估。在本案例中,将探讨一个名为wine_reviewer的项目,该项目的目标是利用机器学习模型,基于二值化的品尝笔记数据来预测葡萄酒评论的分数。这个项目不仅对于葡萄酒爱好者具有极大的吸引力,同时也为数据分析和机器学习的研究人员提供了实践案例。 首先,要理解的关键词是“机器学习”。机器学习是人工智能的一个分支,它让计算机系统能够通过经验自动地改进性能,而无需人类进行明确的编程。在葡萄酒评分预测的场景中,机器学习算法将从大量的葡萄酒品尝笔记数据中学习,发现笔记与葡萄酒最终评分之间的相关性,并利用这种相关性对新的品尝笔记进行评分预测。 接下来是“二值化”处理。在机器学习中,数据预处理是一个重要的步骤,它直接影响模型的性能。二值化是指将数值型数据转换为二进制形式(0和1)的过程,这通常用于简化模型的计算复杂度,或者是数据分类问题中的一种技术。在葡萄酒品尝笔记的上下文中,二值化可能涉及将每种口感、香气和外观等属性的存在与否标记为1(存在)或0(不存在)。这种方法有利于将文本数据转换为机器学习模型可以处理的格式。 葡萄酒评论分数是葡萄酒评估的量化指标,通常由品酒师根据酒的品质、口感、香气、外观等进行评分。在这个项目中,葡萄酒的品尝笔记将被用作特征,而品酒师给出的分数则是目标变量,模型的任务是找出两者之间的关系,并对新的品尝笔记进行分数预测。 在机器学习中,通常会使用多种算法来构建预测模型,如线性回归、决策树、随机森林、梯度提升机等。在wine_reviewer项目中,可能会尝试多种算法,并通过交叉验证等技术来评估模型的性能,最终选择最适合这个任务的模型。 对于这个项目来说,数据集的质量和特征工程将直接影响模型的准确性和可靠性。在准备数据时,可能需要进行数据清洗、缺失值处理、文本规范化、特征选择等步骤。数据集中的标签(目标变量)即为葡萄酒的评分,而特征则来自于品酒师的品尝笔记。 项目还提到了“kaggle”和“R”,这两个都是数据分析和机器学习领域中常见的元素。Kaggle是一个全球性的数据科学竞赛平台,提供各种机器学习挑战和数据集,吸引了来自全球的数据科学家和机器学习专家。通过参与Kaggle竞赛,可以提升个人技能,并有机会接触到最新的机器学习技术和数据处理方法。R是一种用于统计计算和图形的编程语言和软件环境,它在统计分析、数据挖掘、机器学习等领域有广泛的应用。使用R语言可以帮助研究人员进行数据处理、统计分析和模型建立。 至于“压缩包子文件的文件名称列表”,这里可能存在误解或打字错误。通常,这类名称应该表示存储项目相关文件的压缩包,例如“wine_reviewer-master.zip”。这个压缩包可能包含了项目的源代码、数据集、文档和其它相关资源。在开始项目前,研究人员需要解压这个文件包,并且仔细阅读项目文档,以便了解项目的具体要求和数据格式。 总之,wine_reviewer项目是一个结合了机器学习、数据处理和葡萄酒品鉴的有趣尝试,它不仅展示了机器学习在实际生活中的应用潜力,也为研究者提供了丰富的学习资源和实践机会。通过这种跨领域的合作,可以为葡萄酒行业带来更客观、一致的评价标准,并帮助消费者做出更加明智的选择。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依