检查下面代码中的错误并修改# -- coding: UTF-8 -- import numpy as np import matplotlib.pyplot as plt # 设置信道错误概率p和码率R p = 0.01 R = 1000 # 设置典型集参数epsilon和n epsilon = 0.01 n = 10 # 生成随机序列 np_random_sequence = np.random.randint(0, 10000, size=R) # 接收比特 data = np.zeros((R, n)) data[np.where(np.random.rand(R, n) < p)] = 1 # 进行解码 decoded = np.zeros((R, n)) decoded[np.where(np.random.rand(R, n) < p)] = 1 # 仿真性能 P_a = np.sum(data != decoded, axis=1) / n * 100 P_b = np.sum(data != decoded, axis=1) / n * 100 # 绘制性能曲线 plt.plot(np.arange(1, R+1, 1) * R / n, P_a, label='A') plt.xlabel('Time (s)') plt.ylabel('错误率') plt.title('Random 编码 (p = 0.01, R = 1000)') plt.show() # 绘制性能曲线 plt.plot(np.arange(1, R+1, 1) * R / n, P_b, label='B') plt.xlabel('Time (s)') plt.ylabel('错误率') plt.title('Random 编码 (p = 0.01, R = 1000)') plt.show()

时间: 2024-04-28 19:19:20 浏览: 17
# -- coding: UTF-8 -- import numpy as np import matplotlib.pyplot as plt # 设置信道错误概率p和码率R p = 0.01 R = 1000 # 设置典型集参数epsilon和n epsilon = 0.01 n = 10 # 生成随机序列 np_random_sequence = np.random.randint(0, 10000, size=R) # 接收比特 data = np.zeros((R, n)) data[np.where(np.random.rand(R, n) < p)] = 1 # 进行解码 decoded = np.zeros((R, n)) decoded[np.where(np.random.rand(R, n) < p)] = 1 # 仿真性能 P_a = np.sum(data != decoded, axis=1) / n * 100 P_b = np.sum(data != decoded, axis=1) / R * 100 # 绘制性能曲线 plt.plot(np.arange(1, R+1, 1) * n / R, P_a, label='A') plt.xlabel('Time (s)') plt.ylabel('错误率') plt.title('Random 编码 (p = 0.01, R = 1000)') plt.show() # 绘制性能曲线 plt.plot(np.arange(1, R+1, 1) * n / R, P_b, label='B') plt.xlabel('Time (s)') plt.ylabel('错误率') plt.title('Random 编码 (p = 0.01, R = 1000)') plt.show() 修改后代码中的错误为: 1. 第23行中应该除以R而不是n,因为每次接收的data都是R个比特,而不是n个; 2. 第30行和第36行中应该用n除以R而不是R除以n,因为时间应该是按照接收比特的数量来计算的。

相关推荐

#coding:utf8 import pandas as pd import matplotlib.pyplot as plt import seaborn as sns import numpy as np sns.set_style('darkgrid') import matplotlib.font_manager as fm myfont=fm.FontProperties(fname=r'./data/simhei.ttf') #请完善下面的函数 def push_week(new_data): ############ Begin ############ new_data=new_data[new_data['type']==4].copy() #选取样本 new_data['weekdays'] = pd.to_datetime(new_data['time']).apply(lambda x: x.weekday()+1) #时间转化 week_days = new_data.groupby('weekdays')['user_id'].count() #统计购买次数 fig=plt.figure(figsize=(8,6)) #设置大小 bar_width = 0.33 # 设置宽度 plt.bar(week_days.index.values , week_days.values, bar_width, label='下单的次数') plt.xlabel('时间',fontproperties=myfont,fontsize=9) plt.ylabel('数量',fontproperties=myfont,fontsize=9) plt.title('一周内每天的下单情况',fontproperties=myfont,fontsize=12) plt.xticks(week_days.index.values, ('周一', '周二', '周三', '周四', '周五', '周六', '周日'),fontproperties=myfont,fontsize=9) plt.ylim(0,300) plt.legend(prop=myfont) ############ End ############ plt.savefig('./task2/task2_week.png') plt.close(fig) def push_date(new_data): new_data = new_data[(new_data['type'] == 4) & (pd.to_datetime(new_data['time']) < pd.to_datetime('2016-03-01'))].copy() #选出2016年数据 new_data['days'] = [x.day for x in pd.to_datetime(new_data['time'])] #选出天数 renew=new_data.groupby('days')['sku_id'].count() fig = plt.figure(figsize=(8, 6)) plt.plot(renew.index.values,renew.values,label='购买次数') plt.xlabel('天数',fontproperties=myfont,fontsize=9) plt.ylabel('次数',fontproperties=myfont,fontsize=9) plt.title('购买量和月内日期的关系',fontproperties=myfont,fontsize=12) plt.legend(prop=myfont) ############ End ############ plt.savefig('./task2/task2_date.png') plt.close(fig) 报错src/task2_test.py:22: FutureWarning: The pandas.datetime class is deprecated and will be removed from pandas in a future version. Import from datetime instead. data['weekdays'] = pd.to_datetime(data['time']).apply(pd.datetime.weekday) + 1 购买意愿与星期之间的关系图完成! 购买意愿与日期之间的关系图完成!

# -*- coding: utf-8 -*- """ Created on Fri Apr 23 21:10:25 2021 例题:我们把(2,0),(0,2),(0,0)这三个点当作类别1; (3,0),(0,3),(3,3)这三个点当作类别2, 训练好SVM分类器之后,我们预测(-1,-1),(4,4)这两个点所属的类别。 @author: Administrator """ import numpy as np from sklearn.svm import SVC import matplotlib.pyplot as plt data = np.array([[2,0,1],[0,2,1],[0,0,1],[3,0,2],[0,3,2],[3,3,2]]) x = np.array(data[:, 0:2]) y = np.array(data[:,2]) model = SVC(kernel='linear') model.fit(x,y) # ============================================================================= # print(model.dual_coef_) #决策函数中支持向量的系数 # print(model.coef_) #赋予特征的权重(原始问题中的系数)。这仅适用于线性内核 # print(model.intercept_) # 决策函数中的常量 # print(model.support_) #支持向量索引 # print(model.n_support_) #每一类的支持向量数目 print(model.support_vectors_) #支持向量 # ============================================================================= Cp = [[-1,-1],[4,4]] pre = model.predict(Cp) #对Cp中的点进行类别预测 print(pre) plt.scatter(x[:, 0], x[:, 1], c=y, s=30, cmap=plt.cm.Paired) # plot the decision function ax = plt.gca() xlim = ax.get_xlim() ylim = ax.get_ylim() # create grid to evaluate model xx = np.linspace(xlim[0], xlim[1], 30) yy = np.linspace(ylim[0], ylim[1], 30) YY, XX = np.meshgrid(yy, xx) xy = np.vstack([XX.ravel(), YY.ravel()]).T Z = model.decision_function(xy).reshape(XX.shape) # plot decision boundary and margins ax.contour(XX, YY, Z, colors='k', levels=[0], alpha=1, linestyles=['-']) # plot support vectors ax.scatter(model.support_vectors_[:, 0], model.support_vectors_[:, 1], s=100, linewidth=1, facecolors='none', edgecolors='k') plt.show()代码解释

请将下面代码中的两张波动率微笑的图绘制到一张图中# -- coding: utf-8 -- """ Created on Sun May 28 18:08:36 2023 @author: lll """ import numpy as np import matplotlib.pyplot as plt from scipy.optimize import brentq from scipy.stats import norm # 定义BS模型计算期权价格的函数 def bs_price(S, K, r, T, sigma, option='call'): d1 = (np.log(S/K) + (r + 0.5sigma**2)T) / (sigmanp.sqrt(T)) d2 = d1 - sigmanp.sqrt(T) if option == 'call': price = Snorm.cdf(d1) - Knp.exp(-rT)norm.cdf(d2) else: price = Knp.exp(-rT)norm.cdf(-d2) - Snorm.cdf(-d1) return price # 定义计算隐含波动率的函数 def implied_vol(S, K, r, T, price, option='call'): def f(sigma): return bs_price(S, K, r, T, sigma, option) - price return brentq(f, 0.001, 10) # 定义计算波动率微笑图形的函数 def smile_vol(S, r, T, vol_list, K_list, option='call'): if len(vol_list) != len(K_list): raise ValueError("Lengths of vol_list and K_list must be equal.") implied_vol_list = [] for K, vol in zip(K_list, vol_list): price = bs_price(S, K, r, T, vol, option) implied_vol_list.append(implied_vol(S, K, r, T, price, option)) plt.plot(K_list, implied_vol_list) plt.xlabel('Strike') plt.ylabel('Implied Volatility') plt.title('Volatility Smile') plt.show() # 示例代码 S = 100 r = 0.05 T = 1 K_list = np.arange(80, 121, 9) vol_call_list = [0.2, 0.3, 0.4, 0.5, 0.6] vol_put_list = [0.6, 0.5, 0.4, 0.3, 0.2] smile_vol(S, r, T, vol_call_list, K_list, option='call') smile_vol(S, r, T, vol_put_list, K_list, option='put')

最新推荐

recommend-type

Python基于matplotlib实现绘制三维图形功能示例

import matplotlib.pyplot as plt import mpl_toolkits.mplot3d x, y = np.mgrid[-2:2:20j,-2:2:20j] z = x*np.exp(-x**2-y**2) ax = plt.subplot(111, projection='3d') ax.set_title('www.jb51.net - matplotlib...
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

基于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的冬奥会科普平台项目结合了现代信息技术和软件工程的最佳实践,旨在通过信息化手段提高科普效率,为用户提供便捷、高效的科普信息管理服务。
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。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依
recommend-type

Python字符串转Float最佳实践:从初学者到专家的进阶指南

![Python字符串转Float最佳实践:从初学者到专家的进阶指南](https://img-blog.csdnimg.cn/img_convert/1678da8423d7b3a1544fd4e6457be4d1.png) # 1. Python字符串转Float基础** Python中字符串转Float的本质是将文本表示的数字转换为浮点数。这在数据处理、科学计算和许多其他应用中至关重要。本章将介绍字符串转Float的基础知识,包括: * **字符串转Float的意义:**理解字符串和浮点数之间的差异,以及为什么需要进行转换。 * **内置函数:**探索float()函数和decima