# 计算带有自定义频率和权重的声压级 fp_tab = [ (25, 1), (31.5, 0.891), (40, 0.794), (50, 0.708), (63, 0.631), (80, 0.562), (100, 0.501), (125, 0.447), (160, 0.398), (200, 0.355), (250, 0.316), (315, 0.282), (400, 0.251), (500, 0.224), (630, 0.200), (800, 0.178), (1000, 0.158), (1250, 0.141), (1600, 0.126), (2000, 0.112), (2500, 0.100), (3150, 0.089), (4000, 0.079), (5000, 0.071), (6300, 0.063), (8000, 0.056), (10000, 0.050), (12500, 0.045), (16000, 0.040) ] spl = [] for f, w in fp_tab: mask = (freqs >= f - 1/6*f) & (freqs < f + 1/6*f) spl.append(np.max(fft_db[mask])) # 绘制图形 plt.stairs(range(len(fp_tab)), spl) plt.xticks(range(len(fp_tab)), [str(f) for f, _ in fp_tab])报错raise ValueError('Size mismatch between "values" and "edges". ' ValueError: Size mismatch between "values" and "edges". Expected `len(values) + 1 == len(edges)`, but `len(values) = 29` and `len(edges) = 29`.帮我修改正确的代码

时间: 2023-09-02 18:14:52 浏览: 37
这段代码是计算带有自定义频率和权重的声压级,并绘制图形。出现错误的原因是在绘制图形时,`plt.stairs`函数的`values`和`edges`参数长度不一致。你需要将`edges`的长度加1,使其等于`values`的长度加1。 下面是修改后的代码: ```python import numpy as np import matplotlib.pyplot as plt freqs = np.arange(1, 20001) # 假设fft_db是一个长度为20000的numpy数组 fft_db = np.random.rand(20000) fp_tab = [ (25, 1), (31.5, 0.891), (40, 0.794), (50, 0.708), (63, 0.631), (80, 0.562), (100, 0.501), (125, 0.447), (160, 0.398), (200, 0.355), (250, 0.316), (315, 0.282), (400, 0.251), (500, 0.224), (630, 0.200), (800, 0.178), (1000, 0.158), (1250, 0.141), (1600, 0.126), (2000, 0.112), (2500, 0.100), (3150, 0.089), (4000, 0.079), (5000, 0.071), (6300, 0.063), (8000, 0.056), (10000, 0.050), (12500, 0.045), (16000, 0.040) ] spl = [] for f, w in fp_tab: mask = (freqs >= f - 1/6*f) & (freqs < f + 1/6*f) spl.append(np.max(fft_db[mask])) # 修改这里 plt.stairs(range(len(fp_tab)+1), [0]+spl+[0]) plt.xticks(range(len(fp_tab)), [str(f) for f, _ in fp_tab]) plt.show() ``` 修改方法是在`plt.stairs`函数的`values`参数前面添加一个0,同时将`edges`参数的长度加1,使其等于`values`的长度加1。这样就可以正常绘制图形了。

相关推荐

这段代码中加一个test loss功能 class LSTM(nn.Module): def __init__(self, input_size, hidden_size, num_layers, output_size, batch_size, device): super().__init__() self.device = device self.input_size = input_size self.hidden_size = hidden_size self.num_layers = num_layers self.output_size = output_size self.num_directions = 1 # 单向LSTM self.batch_size = batch_size self.lstm = nn.LSTM(self.input_size, self.hidden_size, self.num_layers, batch_first=True) self.linear = nn.Linear(65536, self.output_size) def forward(self, input_seq): h_0 = torch.randn(self.num_directions * self.num_layers, self.batch_size, self.hidden_size).to(self.device) c_0 = torch.randn(self.num_directions * self.num_layers, self.batch_size, self.hidden_size).to(self.device) output, _ = self.lstm(input_seq, (h_0, c_0)) pred = self.linear(output.contiguous().view(self.batch_size, -1)) return pred if __name__ == '__main__': # 加载已保存的模型参数 saved_model_path = '/content/drive/MyDrive/危急值/model/dangerous.pth' device = 'cuda:0' lstm_model = LSTM(input_size=1, hidden_size=64, num_layers=1, output_size=3, batch_size=256, device='cuda:0').to(device) state_dict = torch.load(saved_model_path) lstm_model.load_state_dict(state_dict) dataset = ECGDataset(X_train_df.to_numpy()) dataloader = DataLoader(dataset, batch_size=256, shuffle=True, num_workers=0, drop_last=True) loss_fn = nn.CrossEntropyLoss() optimizer = optim.SGD(lstm_model.parameters(), lr=1e-4) for epoch in range(200000): print(f'epoch:{epoch}') lstm_model.train() epoch_bar = tqdm(dataloader) for x, y in epoch_bar: optimizer.zero_grad() x_out = lstm_model(x.to(device).type(torch.cuda.FloatTensor)) loss = loss_fn(x_out, y.long().to(device)) loss.backward() epoch_bar.set_description(f'loss:{loss.item():.4f}') optimizer.step() if epoch % 100 == 0 or epoch == epoch - 1: torch.save(lstm_model.state_dict(), "/content/drive/MyDrive/危急值/model/dangerous.pth") print("权重成功保存一次")

import numpy as np class BPNeuralNetwork: def __init__(self, input_size, hidden_size, output_size): self.input_size = input_size self.hidden_size = hidden_size self.output_size = output_size # 初始化权重和偏置 self.weights_ih = np.random.randn(hidden_size, input_size) self.bias_ih = np.random.randn(hidden_size, 1) self.weights_ho = np.random.randn(output_size, hidden_size) self.bias_ho = np.random.randn(output_size, 1) # 定义激活函数 self.activation = lambda x: 1 / (1 + np.exp(-x)) self.derivative = lambda x: x * (1 - x) def forward(self, inputs): # 计算隐藏层的输出 hidden = self.activation(np.dot(self.weights_ih, inputs) + self.bias_ih) # 计算输出层的输出 output = self.activation(np.dot(self.weights_ho, hidden) + self.bias_ho) return output def backward(self, inputs, targets, output): # 计算输出层的误差 output_error = targets - output output_delta = output_error * self.derivative(output) # 计算隐藏层的误差 hidden_error = np.dot(self.weights_ho.T, output_delta) hidden_delta = hidden_error * self.derivative(hidden) # 更新权重和偏置 self.weights_ho += np.dot(output_delta, hidden.T) self.bias_ho += output_delta self.weights_ih += np.dot(hidden_delta, inputs.T) self.bias_ih += hidden_delta def train(self, inputs, targets, epochs): for i in range(epochs): for j in range(len(inputs)): # 前向传播 output = self.forward(inputs[j].reshape(-1, 1)) # 反向传播 self.backward(inputs[j].reshape(-1, 1), targets[j].reshape(-1, 1), output)

最新推荐

recommend-type

在keras中model.fit_generator()和model.fit()的区别说明

它们都用于更新模型的权重以最小化损失函数,但针对不同类型的输入数据和场景有不同的适用性。 首先,`model.fit()`是Keras中最基础的训练接口,它期望输入数据`x_train`和目标变量`y_train`已经完全加载到内存中。...
recommend-type

对tensorflow中tf.nn.conv1d和layers.conv1d的区别详解

在TensorFlow中,`tf.nn.conv1d`和`layers.conv1d`都是用于执行一维卷积操作的函数,但它们在实现细节和使用上存在一些差异。这篇文章将深入探讨这两个函数的区别,并帮助理解它们在构建一维卷积神经网络(1D CNN)...
recommend-type

美国地图json文件,可以使用arcgis转为spacefile

美国地图json文件,可以使用arcgis转为spacefile
recommend-type

Microsoft Edge 126.0.2592.68 32位离线安装包

Microsoft Edge 126.0.2592.68 32位离线安装包
recommend-type

FLASH源码:读写FLASH内部数据,读取芯片ID

STLINK Utility:读取FLASH的软件
recommend-type

基于Springboot的医院信管系统

"基于Springboot的医院信管系统是一个利用现代信息技术和网络技术改进医院信息管理的创新项目。在信息化时代,传统的管理方式已经难以满足高效和便捷的需求,医院信管系统的出现正是适应了这一趋势。系统采用Java语言和B/S架构,即浏览器/服务器模式,结合MySQL作为后端数据库,旨在提升医院信息管理的效率。 项目开发过程遵循了标准的软件开发流程,包括市场调研以了解需求,需求分析以明确系统功能,概要设计和详细设计阶段用于规划系统架构和模块设计,编码则是将设计转化为实际的代码实现。系统的核心功能模块包括首页展示、个人中心、用户管理、医生管理、科室管理、挂号管理、取消挂号管理、问诊记录管理、病房管理、药房管理和管理员管理等,涵盖了医院运营的各个环节。 医院信管系统的优势主要体现在:快速的信息检索,通过输入相关信息能迅速获取结果;大量信息存储且保证安全,相较于纸质文件,系统节省空间和人力资源;此外,其在线特性使得信息更新和共享更为便捷。开发这个系统对于医院来说,不仅提高了管理效率,还降低了成本,符合现代社会对数字化转型的需求。 本文详细阐述了医院信管系统的发展背景、技术选择和开发流程,以及关键组件如Java语言和MySQL数据库的应用。最后,通过功能测试、单元测试和性能测试验证了系统的有效性,结果显示系统功能完整,性能稳定。这个基于Springboot的医院信管系统是一个实用且先进的解决方案,为医院的信息管理带来了显著的提升。"
recommend-type

管理建模和仿真的文件

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

字符串转Float性能调优:优化Python字符串转Float性能的技巧和工具

![字符串转Float性能调优:优化Python字符串转Float性能的技巧和工具](https://pic1.zhimg.com/80/v2-3fea10875a3656144a598a13c97bb84c_1440w.webp) # 1. 字符串转 Float 性能调优概述 字符串转 Float 是一个常见的操作,在数据处理和科学计算中经常遇到。然而,对于大规模数据集或性能要求较高的应用,字符串转 Float 的效率至关重要。本章概述了字符串转 Float 性能调优的必要性,并介绍了优化方法的分类。 ### 1.1 性能调优的必要性 字符串转 Float 的性能问题主要体现在以下方面
recommend-type

Error: Cannot find module 'gulp-uglify

当你遇到 "Error: Cannot find module 'gulp-uglify'" 这个错误时,它通常意味着Node.js在尝试运行一个依赖了 `gulp-uglify` 模块的Gulp任务时,找不到这个模块。`gulp-uglify` 是一个Gulp插件,用于压缩JavaScript代码以减少文件大小。 解决这个问题的步骤一般包括: 1. **检查安装**:确保你已经全局安装了Gulp(`npm install -g gulp`),然后在你的项目目录下安装 `gulp-uglify`(`npm install --save-dev gulp-uglify`)。 2. **配置
recommend-type

基于Springboot的冬奥会科普平台

"冬奥会科普平台的开发旨在利用现代信息技术,如Java编程语言和MySQL数据库,构建一个高效、安全的信息管理系统,以改善传统科普方式的不足。该平台采用B/S架构,提供包括首页、个人中心、用户管理、项目类型管理、项目管理、视频管理、论坛和系统管理等功能,以提升冬奥会科普的检索速度、信息存储能力和安全性。通过需求分析、设计、编码和测试等步骤,确保了平台的稳定性和功能性。" 在这个基于Springboot的冬奥会科普平台项目中,我们关注以下几个关键知识点: 1. **Springboot框架**: Springboot是Java开发中流行的应用框架,它简化了创建独立的、生产级别的基于Spring的应用程序。Springboot的特点在于其自动配置和起步依赖,使得开发者能快速搭建应用程序,并减少常规配置工作。 2. **B/S架构**: 浏览器/服务器模式(B/S)是一种客户端-服务器架构,用户通过浏览器访问服务器端的应用程序,降低了客户端的维护成本,提高了系统的可访问性。 3. **Java编程语言**: Java是这个项目的主要开发语言,具有跨平台性、面向对象、健壮性等特点,适合开发大型、分布式系统。 4. **MySQL数据库**: MySQL是一个开源的关系型数据库管理系统,因其高效、稳定和易于使用而广泛应用于Web应用程序,为平台提供数据存储和查询服务。 5. **需求分析**: 开发前的市场调研和需求分析是项目成功的关键,它帮助确定平台的功能需求,如用户管理、项目管理等,以便满足不同用户群体的需求。 6. **数据库设计**: 数据库设计包括概念设计、逻辑设计和物理设计,涉及表结构、字段定义、索引设计等,以支持平台的高效数据操作。 7. **模块化设计**: 平台功能模块化有助于代码组织和复用,包括首页模块、个人中心模块、管理系统模块等,每个模块负责特定的功能。 8. **软件开发流程**: 遵循传统的软件生命周期模型,包括市场调研、需求分析、概要设计、详细设计、编码、测试和维护,确保项目的质量和可维护性。 9. **功能测试、单元测试和性能测试**: 在开发过程中,通过这些测试确保平台功能的正确性、模块的独立性和系统的性能,以达到预期的用户体验。 10. **微信小程序、安卓源码**: 虽然主要描述中没有详细说明,但考虑到标签包含这些内容,可能平台还提供了移动端支持,如微信小程序和安卓应用,以便用户通过移动设备访问和交互。 这个基于Springboot的冬奥会科普平台项目结合了现代信息技术和软件工程的最佳实践,旨在通过信息化手段提高科普效率,为用户提供便捷、高效的科普信息管理服务。