clear all; n=0.2; w=[0,0,0]; P=[-9,1,-12,-4,0,5;... 15,-8,4,5,11,9]; d=[0,1,0,0,0,1]; P=[ones(1,6);P]; MAX=20; i=0; while 1 v=w*P; y=hardlim(v); e=(d-y); ee(i+1)=mae(e); if (ee(i+1)<0.001) disp('we have got it:'); disp(w); break; end w=w+n*(d-y)*P'; i=i+1; if (i>=MAX) disp('MAX times loop'); disp(w); disp(ee(i+1)); break; end end figure; subplot(2,1,1); plot([-9,-12,-4,0],[15,4,5,11],'o'); hold on; plot([1,5],[-8,9],'*'); axis([-13,6,-10,16]); legend('第一类','第二类'); title('6个坐标点的二分类'); x=-13:.2:6; y=x*(-w(2)/w(3)-w(1)/w(3)); plot(x,y); hold off; subplot(2,1,2); x=0:i; plot(x,ee,'o-'); s=sprintf('mae的值(迭代次数:%d)',i+1); title(s); 转变成python语言

时间: 2024-04-10 19:33:24 浏览: 21
在Python中,可以将上述MATLAB代码转换为以下代码: ```python import numpy as np import matplotlib.pyplot as plt n = 0.2 w = np.zeros(3) P = np.array([[-9, 1, -12, -4, 0, 5], [15, -8, 4, 5, 11, 9]]) d = np.array([0, 1, 0, 0, 0, 1]) P = np.vstack((np.ones(6), P)) MAX = 20 i = 0 ee = [] while True: v = np.dot(w, P) y = np.heaviside(v, 0) e = d - y ee.append(np.mean(np.abs(e))) if ee[i] < 0.001: print('we have got it:') print(w) break w = w + n * np.dot((d - y), P.T) i += 1 if i >= MAX: print('MAX times loop') print(w) print(ee[i]) break plt.figure() plt.subplot(2, 1, 1) plt.plot([-9, -12, -4, 0], [15, 4, 5, 11], 'o') plt.hold(True) plt.plot([1, 5], [-8, 9], '*') plt.axis([-13, 6, -10, 16]) plt.legend(['Class 1', 'Class 2']) plt.title('Binary Classification of 6 Coordinate Points') x = np.arange(-13, 6, 0.2) y = x * (-w[1] / w[2] - w[0] / w[2]) plt.plot(x, y) plt.hold(False) plt.subplot(2, 1, 2) x = np.arange(0, i + 1) plt.plot(x, ee, 'o-') s = f'mae value (iterations: {i + 1})' plt.title(s) plt.show() ``` 这是将MATLAB代码转换为Python代码的基本过程。请注意,在Python中,需要使用NumPy库进行矩阵运算,并使用Matplotlib库进行绘图。另外,一些函数和语法略有不同,因此可能需要进行适当的调整和修改。

相关推荐

clear all; clc; source_img=imread('C:\Users\LENOVO\Desktop\yes.jpg');%读取图片 [m,n,p]=size(source_img);%计算图片的行数列数层数 %==========从RGB转换到HSV======================= hsv_img=rgb2hsv(source_img); h=hsv_img(:,:,1); s=hsv_img(:,:,2); v=hsv_img(:,:,3); figure; subplot(221);imshow(source_img); subplot(222);imshow(h); subplot(223);imshow(s); subplot(224);imshow(v); %============V分量小波包分解======================================== [cc,ss]=wavedec2(v,1,'haar'); cA=appcoef2(cc,ss,'haar',1); %cc:小波分解的小波系数矩阵;ss:小波分解对应的尺度矩阵;分解的层数为1 cH=detcoef2('h',cc,ss,1); %h:提取水平高频;v:垂直高频;d:对角高频 cV=detcoef2('v',cc,ss,1); cD=detcoef2('d',cc,ss,1); cA1=mapminmax(cA,0,1);%归一化处理 figure; subplot(221);imshow(cA1,[]);title('(a) 近似分量cA'); subplot(222);imshow(cH,[]);title('(b) 细节分量cH'); subplot(223);imshow(cV,[]);title('(c) 细节分量cV'); subplot(224);imshow(cD,[]);title('(d) 细节分量cD'); %=============近似分量cA双边滤波================================== w = 3; % bilateral filter half-width sigma = [3 0.2]; % bilateral filter standard deviations cA2=bfilter2(cA1,w,sigma); %双边滤波 hsize=15; sigma1=15; sigma2=85; sigma3=265; H1=fspecial('gaussian',hsize,sigma1); H2=fspecial('gaussian',hsize,sigma2); H3=fspecial('gaussian',hsize,sigma3); img1=conv2(v,H1,'same'); img2=conv2(v,H2,'same'); img3=conv2(v,H3,'same'); imggaus=1/3*img1+1/3*img2+1/3*img3; alpha5=0.1; k5=alpha5*sum(s(:))/(m*n); vnew5=v*(1+k5)./(max(v,imggaus)+k5); X1=v(:)'; X2=vnew5(:)'; X=[X1 X2]; C=cov(X1,X2); [V,D]=eig(C); diagD=diag(D); if diagD(1)>diagD(2) V1=V(:,1); else V1=V(:,2); end w1=V1(1)/(V1(1)+V1(2)); w2=V1(2)/(V1(1)+V1(2)); recon_set=w1*v+w2*vnew5; %==============HSV转换RGB========================================= hsv(:,:,1)=h; hsv(:,:,2)=s; hsv(:,:,3)=recon_set; rgb_img=hsv2rgb(hsv); figure; subplot(121);imshow(source_img);%显示原图 subplot(122);imshow(rgb_img);%显示增强后的图 根据代码写结果说明及分析

最新推荐

recommend-type

k8s1.16的jenkins部署java项目cicd(cd手动)-kubernetes安装包和详细文档笔记整理

k8s1.16的jenkins部署java项目cicd(cd手动)-kubernetes安装包和详细文档笔记整理
recommend-type

sja1311.x86_64.tar.gz

SQLyong 各个版本,免费下载 SQLyog是业界著名的Webyog公司出品的一款简洁高效、功能强大的图形化MySQL数据库管理工具。使用SQLyog可以快速直观地让您从世界的任何角落通过网络来维护远端的MySQL数据库。
recommend-type

debugpy-1.1.0-cp34-cp34m-manylinux1_x86_64.whl

Python库是一组预先编写的代码模块,旨在帮助开发者实现特定的编程任务,无需从零开始编写代码。这些库可以包括各种功能,如数学运算、文件操作、数据分析和网络编程等。Python社区提供了大量的第三方库,如NumPy、Pandas和Requests,极大地丰富了Python的应用领域,从数据科学到Web开发。Python库的丰富性是Python成为最受欢迎的编程语言之一的关键原因之一。这些库不仅为初学者提供了快速入门的途径,而且为经验丰富的开发者提供了强大的工具,以高效率、高质量地完成复杂任务。例如,Matplotlib和Seaborn库在数据可视化领域内非常受欢迎,它们提供了广泛的工具和技术,可以创建高度定制化的图表和图形,帮助数据科学家和分析师在数据探索和结果展示中更有效地传达信息。
recommend-type

SQLyog-12.5.0-0.x64Community.exe

SQLyog-12.5.0-0.x64Community
recommend-type

debugpy-1.0.0b12-cp35-cp35m-manylinux1_i686.whl

Python库是一组预先编写的代码模块,旨在帮助开发者实现特定的编程任务,无需从零开始编写代码。这些库可以包括各种功能,如数学运算、文件操作、数据分析和网络编程等。Python社区提供了大量的第三方库,如NumPy、Pandas和Requests,极大地丰富了Python的应用领域,从数据科学到Web开发。Python库的丰富性是Python成为最受欢迎的编程语言之一的关键原因之一。这些库不仅为初学者提供了快速入门的途径,而且为经验丰富的开发者提供了强大的工具,以高效率、高质量地完成复杂任务。例如,Matplotlib和Seaborn库在数据可视化领域内非常受欢迎,它们提供了广泛的工具和技术,可以创建高度定制化的图表和图形,帮助数据科学家和分析师在数据探索和结果展示中更有效地传达信息。
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

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

MATLAB柱状图在信号处理中的应用:可视化信号特征和频谱分析

![matlab画柱状图](https://img-blog.csdnimg.cn/3f32348f1c9c4481a6f5931993732f97.png) # 1. MATLAB柱状图概述** MATLAB柱状图是一种图形化工具,用于可视化数据中不同类别或组的分布情况。它通过绘制垂直条形来表示每个类别或组中的数据值。柱状图在信号处理中广泛用于可视化信号特征和进行频谱分析。 柱状图的优点在于其简单易懂,能够直观地展示数据分布。在信号处理中,柱状图可以帮助工程师识别信号中的模式、趋势和异常情况,从而为信号分析和处理提供有价值的见解。 # 2. 柱状图在信号处理中的应用 柱状图在信号处理
recommend-type

用Spring boot和vue写一个登录注册界面

好的,这是一个比较常见的技术问题,我可以为您提供一些参考资料和思路。首先,您可以使用Spring Boot作为后端框架来处理用户认证和注册的逻辑,比如使用Spring Security实现用户登录认证。同时,您还需要设计相应的数据模型和数据库表结构来存储用户信息。在前端方面,您可以使用Vue.js作为框架来构建登录注册页面,使用Axios来发起API请求并和后端进行交互。当然,在实现过程中,还需要考虑一些具体细节,比如数据校验、安全性和用户体验等方面。希望这些信息能够帮助到您。
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。