import numpy as np import cv2 as cv import matplotlib.pyplot as plt a = np.random.randint(90,96,size=(20,2)).astype(np.float32) b = np.random.randint(95,101,size=(20,2)).astype(np.float32) data = np.vstack((a,b)) data = np.array(data,dtype=np.float32) a_label = np.zeros((20,1)) b_label = np.ones((20,1)) label = np.vstack((a_label,b_label)) label = np.array(label,dtype=np.float32) svm = cv.ml.SVM_create() svm.setType(cv.ml.SVM_C_SVC) svm.setKernel(cv.ml.SVM_LINEAR) svm.setC(0.01) result = svm.train(data,cv.ml.ROW_SAMPLE,label) test = np.array([[93.92, 96.97]], dtype=np.float32) p1, p2 = svm.predict(test) plt.scatter(a[:,0],a[:,1],s=80,c="g",marker="o") plt.scatter(b[:,0],b[:,1],s=80,c="b",marker="s") plt.scatter(test[:,0],test[:,1],s=80,c="r",marker="*") plt.show() print(p1, p2)为什么代码会报下列错误n the case of classification problem the responses must be categorical; either specify varType when creating TrainData, or pass integer responses in function 'cv::ml::SVMImpl::train'

时间: 2023-11-29 18:10:17 浏览: 176
这个错误说明在训练SVM模型时,响应变量(即标签)不是分类变量,而是连续变量或其他类型的变量。因此,需要将响应变量转换为分类变量才能训练SVM模型。 为了解决这个问题,您需要将标签变量转换为分类变量。您可以使用OpenCV的ml.TrainData()函数来创建训练数据对象,并将响应变量的类型指定为CV_VAR_CATEGORICAL。具体来说,您可以在创建训练数据对象时指定varType参数,如下所示: ```python train_data = cv.ml.TrainData_create(data, cv.ml.ROW_SAMPLE, label) train_data.setVarType(np.array([cv.ml.VAR_NUMERICAL, cv.ml.VAR_NUMERICAL], dtype=np.uint8)) ``` 在这个例子中,我们将响应变量的类型设置为CV_VAR_CATEGORICAL,将特征变量的类型设置为CV_VAR_NUMERICAL,因为特征变量是连续变量。 下面是修改后的代码: ```python import numpy as np import cv2 as cv import matplotlib.pyplot as plt a = np.random.randint(90,96,size=(20,2)).astype(np.float32) b = np.random.randint(95,101,size=(20,2)).astype(np.float32) data = np.vstack((a,b)) data = np.array(data,dtype=np.float32) a_label = np.zeros((20,1)) b_label = np.ones((20,1)) label = np.vstack((a_label,b_label)) label = np.array(label,dtype=np.float32) train_data = cv.ml.TrainData_create(data, cv.ml.ROW_SAMPLE, label) train_data.setVarType(np.array([cv.ml.VAR_NUMERICAL, cv.ml.VAR_NUMERICAL], dtype=np.uint8)) svm = cv.ml.SVM_create() svm.setType(cv.ml.SVM_C_SVC) svm.setKernel(cv.ml.SVM_LINEAR) svm.setC(0.01) result = svm.train(train_data) test = np.array([[93.92, 96.97]], dtype=np.float32) p1, p2 = svm.predict(test) plt.scatter(a[:,0],a[:,1],s=80,c="g",marker="o") plt.scatter(b[:,0],b[:,1],s=80,c="b",marker="s") plt.scatter(test[:,0],test[:,1],s=80,c="r",marker="*") plt.show() print(p1, p2) ``` 希望这个修改可以帮助您解决问题!
阅读全文

相关推荐

代码改进:import numpy as np import pandas as pd import matplotlib as mpl import matplotlib.pyplot as plt from sklearn.datasets import make_blobs def distEclud(arrA,arrB): #欧氏距离 d = arrA - arrB dist = np.sum(np.power(d,2),axis=1) #差的平方的和 return dist def randCent(dataSet,k): #寻找质心 n = dataSet.shape[1] #列数 data_min = dataSet.min() data_max = dataSet.max() #生成k行n列处于data_min到data_max的质心 data_cent = np.random.uniform(data_min,data_max,(k,n)) return data_cent def kMeans(dataSet,k,distMeans = distEclud, createCent = randCent): x,y = make_blobs(centers=100)#生成k质心的数据 x = pd.DataFrame(x) m,n = dataSet.shape centroids = createCent(dataSet,k) #初始化质心,k即为初始化质心的总个数 clusterAssment = np.zeros((m,3)) #初始化容器 clusterAssment[:,0] = np.inf #第一列设置为无穷大 clusterAssment[:,1:3] = -1 #第二列放本次迭代点的簇编号,第三列存放上次迭代点的簇编号 result_set = pd.concat([pd.DataFrame(dataSet), pd.DataFrame(clusterAssment)],axis = 1,ignore_index = True) #将数据进行拼接,横向拼接,即将该容器放在数据集后面 clusterChanged = True while clusterChanged: clusterChanged = False for i in range(m): dist = distMeans(dataSet.iloc[i,:n].values,centroids) #计算点到质心的距离(即每个值到质心的差的平方和) result_set.iloc[i,n] = dist.min() #放入距离的最小值 result_set.iloc[i,n+1] = np.where(dist == dist.min())[0] #放入距离最小值的质心标号 clusterChanged = not (result_set.iloc[:,-1] == result_set.iloc[:,-2]).all() if clusterChanged: cent_df = result_set.groupby(n+1).mean() #按照当前迭代的数据集的分类,进行计算每一类中各个属性的平均值 centroids = cent_df.iloc[:,:n].values #当前质心 result_set.iloc[:,-1] = result_set.iloc[:,-2] #本次质心放到最后一列里 return centroids, result_set x = np.random.randint(0,100,size=100) y = np.random.randint(0,100,size=100) randintnum=pd.concat([pd.DataFrame(x), pd.DataFrame(y)],axis = 1,ignore_index = True) #randintnum_test, randintnum_test = kMeans(randintnum,3) #plt.scatter(randintnum_test.iloc[:,0],randintnum_test.iloc[:,1],c=randintnum_test.iloc[:,-1]) #result_test,cent_test = kMeans(data, 4) cent_test,result_test = kMeans(randintnum, 3) plt.scatter(result_test.iloc[:,0],result_test.iloc[:,1],c=result_test.iloc[:,-1]) plt.scatter(cent_test[:,0],cent_test[:,1],color = 'red',marker = 'x',s=100)

请在不影响结果的条件下改变代码的样子:import numpy as np import matplotlib.pyplot as plt x1len = 21 x2len = 18 LEN = x1len + x2len POPULATION_SIZE = 100 GENERATIONS = 251 CROSSOVER_RATE = 0.7 MUTATION_RATE = 0.3 pop = np.random.randint(0,2,size=(POPULATION_SIZE,LEN)) def BinToX(pop): x1 = pop[:,0:x1len] x2 = pop[:,x1len:] x1 = x1.dot(2**np.arange(x1len)[::-1]) x2 = x2.dot(2**np.arange(x2len)[::-1]) x1 = -2.9 + x1*(12 + 2.9)/(np.power(2,x1len)-1) x2 = 4.2 + x2*(5.7 - 4.2)/(np.power(2,x2len)-1) return x1,x2 def func(pop): x1,x2 = BinToX(pop) return 21.5 + x1*np.sin(4*np.pi*x1) + x2*np.sin(20*np.pi*x2) def fn(pop): return func(pop); def selection(pop, fitness): idx = np.random.choice(np.arange(pop.shape[0]), size=POPULATION_SIZE, replace=True, p=fitness/fitness.sum()) return pop[idx] def crossover(IdxP1,pop): if np.random.rand() < CROSSOVER_RATE: C = np.zeros((1,LEN)) IdxP2 = np.random.randint(0, POPULATION_SIZE) pt = np.random.randint(0, LEN) C[0,:pt] = pop[IdxP1,:pt] C[0,pt:] = pop[IdxP2, pt:] np.append(pop, C, axis=0) return def mutation(idx,pop): if np.random.rand() < MUTATION_RATE: mut_index = np.random.randint(0, LEN) pop[idx,mut_index] = 1- pop[idx,mut_index] return best_chrom = np.zeros(LEN) best_score = 0 fig = plt.figure() for generation in range(GENERATIONS): fitness = fn(pop) pop = selection(pop, fitness) if generation%50 == 0: ax = fig.add_subplot(2,3,generation//50 +1, projection='3d', title = "generation:"+str(generation)+" best="+str(np.max(fitness))) x1,x2 = BinToX(pop) z = func(pop) ax.scatter(x1,x2,z) for idx in range(POPULATION_SIZE): crossover(idx,pop) mutation(idx,pop) idx = np.argmax(fitness) if best_score < fitness[idx]: best_score = fitness[idx] best_chrom = pop[idx, :] plt.show() print('最优解:', best_chrom, '| best score: %.2f' % best_score)

import matplotlib.pyplot as plt import math import random import numpy as np pop_size = 50 # 种群数量 PC=0.6 # 交叉概率 PM=0.1 #变异概率 X_max=10 #最大值 X_min=0 #最小值 DNA_SIZE=10 #DNA长度与保留位数有关,2**10 当前保留3位小数点 N_GENERATIONS=100 """ 求解的目标表达式为: y = 10 * math.sin(5 * x) + 7 * math.cos(4 * x) x=[0,5] """ def aim(x):return 10*x#np.sin(5*x)+7*np.cos(4*x) def f1(pop): return pop.dot(2 ** np.arange(DNA_SIZE)[::-1]) *(X_max-X_min)/ float(2**DNA_SIZE-1) +X_min def f2(pred): return pred + 1e-3 - np.min(pred) def f3(pop, fitness): idx = np.random.choice(np.arange(pop_size), size=pop_size, replace=True,p=fitness/fitness.sum()) return pop[idx] def f4(parent, pop): if np.random.rand() < PC: i_ = np.random.randint(0, pop_size, size=1) cross_points = np.random.randint(0, 2, size=DNA_SIZE).astype(np.bool) parent[cross_points] = pop[i_, cross_points] return parent def f5(child,pm): for point in range(DNA_SIZE): if np.random.rand() < pm: child[point] = 1 if child[point] == 0 else 0 return child pop = np.random.randint(2, size=(pop_size, DNA_SIZE)) for i in range(N_GENERATIONS): #解码 X_value= ? #获取目标函数值 F_values = ? #获取适应值 fitness = ? if(i==0): max=np.max(F_values) max_DNA = pop[np.argmax(F_values), :] if(max<np.max(F_values)): max=np.max(F_values) max_DNA=pop[np.argmax(F_values), :] if (i % 10 == 0): print("Most fitted value and X: \n", np.max(F_values), decode(pop[np.argmax(F_values), :])) #选择 pop = ? pop_copy = pop.copy() #交叉 变异 for parent in pop: child = ? child = ? parent[:] = child print("目标函数最大值为:",max) print("其DNA值为:",max_DNA) print("其X值为:",decode(max_DNA))

最新推荐

recommend-type

STM32F103单片机连接EC800-4G模块采集GNSS定位数据和多组传感器数据上传到ONENET云平台并接收控制指令.zip

1、嵌入式物联网单片机项目开发例程,简单、方便、好用,节省开发时间。 2、代码使用KEIL 标准库开发,当前在STM32F103运行,如果是STM32F103其他型号芯片,依然适用,请自行更改KEIL芯片型号以及FLASH容量即可。 3、软件下载时,请注意keil选择项是jlink还是stlink。 4、有偿指导v:wulianjishu666; 5、如果接入其他传感器,请查看账号发布的其他资料。 6、单片机与模块的接线,在代码当中均有定义,请自行对照。 7、若硬件有差异,请根据自身情况调整代码,程序仅供参考学习。 8、代码有注释说明,请耐心阅读。
recommend-type

Elasticsearch核心改进:实现Translog与索引线程分离

资源摘要信息:"Elasticsearch是一个基于Lucene构建的开源搜索引擎。它提供了一个分布式多用户能力的全文搜索引擎,基于RESTful web接口。Elasticsearch是用Java语言开发的,并作为Apache许可条款下的开源项目发布,是当前流行的企业级搜索引擎。设计用于云计算中,能够达到实时搜索,稳定,可靠,快速,安装使用方便。" "Elasticsearch的索引线程是处理索引操作的重要部分,负责处理数据的写入、更新和删除等操作。但是,在处理大量数据和高并发请求时,如果索引线程处理速度过慢,就会导致数据处理的延迟,影响整体性能。因此,Elasticsearch采用了事务日志(translog)机制来提高索引操作的效率和可靠性。" "Elasticsearch的事务日志(translog)是一种持久化存储机制,用于记录所有未被持久化到分片中的索引操作。在发生故障或系统崩溃时,事务日志可以确保所有索引操作不会丢失,保证数据的完整性。每个分片都有自己的事务日志文件。" "在Elasticsearch的早期版本中,事务日志的操作和索引线程的操作是在同一个线程中完成的,这可能会导致性能瓶颈。为了解决这个问题,Elasticsearch将事务日志的操作从索引线程中分离出去,使得索引线程可以专注于数据的索引操作,而事务日志的操作可以独立地进行。这样可以大大提高了Elasticsearch的索引性能。" "但是,事务日志的操作是独立于索引操作的,这就需要保证事务日志的操作不会影响到索引操作的性能。因此,在将事务日志从索引线程分离出去的同时,Elasticsearch也引入了一些优化策略,比如批量写入事务日志,减少磁盘I/O操作,以及优化事务日志的数据结构,提高读写效率等。" "需要注意的是,虽然事务日志的分离可以提高索引操作的性能,但是也会增加系统的复杂度和维护难度。因此,开发者在使用这个功能时,需要充分理解其原理和影响,才能确保系统的稳定运行。" "此外,由于这个功能还处于测试和学习阶段,尚未被广泛应用于生产环境,所以开发者在使用时需要谨慎,避免对生产环境造成影响。" "总的来说,Elasticsearch的事务日志的分离是一个重要的优化,可以大大提升索引操作的性能,但是在使用时也需要充分考虑其带来的影响,才能确保系统的稳定运行。"
recommend-type

管理建模和仿真的文件

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

病房呼叫系统设计基础:7个关键架构策略让你一步入门

![病房呼叫系统设计基础:7个关键架构策略让你一步入门](https://zektek.com.mx/wp-content/uploads/2021/03/diagram-enfermeria.jpg) # 摘要 本文对病房呼叫系统进行了深入的概述、需求分析、架构设计、功能实现以及实践应用案例的探讨。通过分析系统架构的重要性、设计原则、模块划分和数据流,确保了系统的高效运行和优化。本文进一步探讨了呼叫信号传输技术、显示与反馈机制、系统安全性与可靠性设计,并分析了系统部署环境、安装调试流程和维护升级策略。最后,文章展望了病房呼叫系统的未来发展趋势,包括智能化、技术融合以及法规遵从与伦理考量,并
recommend-type

Selenium如何获取Shadow DOM下的元素属性?

在Selenium中,获取Shadow DOM下的元素属性通常涉及到两步:首先找到元素,然后访问它的属性。由于Shadow DOM元素默认是不可见的(对于非JavaScript开发者),所以我们需要用JavaScript脚本来获取其内容。 下面是一个示例,展示如何通过Selenium的`execute_script`函数获取Shadow DOM元素的属性: ```python from selenium.webdriver.common.by import By from selenium.webdriver.support.ui import WebDriverWait from sel
recommend-type

分享个人Vim与Git配置文件管理经验

资源摘要信息:"conffiles:我的vim和git配置文件" 在给定的文件信息中,我们可以梳理出一些关键知识点,这些知识点主要涉及到了Vim编辑器和Git版本控制系统,同时涉及到了Linux环境下的一些文件操作知识。 首先,文件标题提到了"conffiles",这通常是指配置文件(configuration files)的缩写。配置文件是软件运行时用于读取用户设置或其他运行参数的文件,它们允许软件按照用户的特定需求进行工作。在本例中,这些配置文件是与Vim编辑器和Git版本控制系统相关的。 Vim是一种流行的文本编辑器,是UNIX系统中vi编辑器的增强版本。Vim不仅支持代码编辑,还支持插件扩展、多种模式(命令模式、插入模式、视觉模式等)和高度可定制化。在这个上下文中,"我的vim"可能指的是使用者为Vim定制的一套配置文件,这些配置文件可能包含键位映射、颜色主题、插件设置、用户界面布局和其他个性化选项。 Git是一个版本控制系统,用于跟踪计算机文件的更改和协作。Git是分布式版本控制,这意味着每个开发者都有一个包含完整项目历史的仓库副本。Git常用于代码的版本控制管理,它允许用户回滚到之前的版本、合并来自不同贡献者的代码,并且有效地管理代码变更。在这个资源中,"git conffiles"可能表示与Git用户相关的配置文件,这可能包括用户凭证、代理设置、别名以及其他一些全局Git配置选项。 描述部分提到了使用者之前使用的编辑器是Vim,但现在转向了Emacs。尽管如此,该用户仍然保留了以前的Vim配置文件。接着,描述中提到了一个安装脚本命令"sh ./.vim/install.sh"。这是一个shell脚本,通常用于自动化安装或配置过程。在这里,这个脚本可能用于创建符号链接(symbolic links),将旧的Vim配置文件链接到当前使用的Emacs配置文件夹中,使用户能够继续使用他们熟悉且习惯的Vim配置。 标签"Vimscript"表明这是一个与Vim脚本相关的资源,Vim脚本是一种专门用于自定义和扩展Vim功能的编程语言。Vimscript可以用于编写宏、自定义函数、插件等。 最后,文件名称列表"conffiles-master"可能表明这个压缩包文件包含了一系列的主配置文件。在Git版本控制的术语中,"master"(现在通常称为"main")分支是项目仓库的默认分支。这暗示了这些配置文件可能是该用户项目的主配置文件,这些配置文件被包含在名为"conffiles-master"的压缩包中。 综上所述,这个资源可能是一个集合了Vim编辑器和Git版本控制系统的个人配置文件的压缩包,附带一个用于符号链接旧Vim配置的安装脚本,它能够帮助用户在转向其他工具时仍然能够使用之前的个性化设置。这个资源对于想要了解如何管理和迁移配置文件的用户具有一定的参考价值。
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

【Genesis 2000教程】:7个技巧助你精通界面布局与操作

![技术专有名词:Genesis 2000](http://intewellos.com/uploadfiles/2023/03/20230321111815165.png) # 摘要 本文全面介绍了Genesis 2000软件的界面布局基础、操作技巧、视觉效果调整、高级功能应用以及综合案例分析,旨在帮助用户高效利用该软件提升工作效率和设计质量。文章首先从界面元素和布局优化入手,讲述了如何定制面板、工具栏以及管理窗口与视图。接着,探讨了通过快捷键和搜索功能实现的高效导航与搜索技巧。第三章强调了视觉效果与图形、文本处理的重要性,并提供实现高级视觉效果的技巧。第四章详细介绍了插件集成、参数化设计
recommend-type

求出所有100到200以内的偶数,并放在数组中,按照每行5个输出

为了求出100到200之间所有的偶数并将它们放入数组中,然后每行打印5个,你可以按以下步骤操作: 1. 创建一个空数组,用于存储偶数。 2. 使用一个for循环,从100开始,每次增加2,直到达到200(包括200)。 - 对于每个数字,检查它是否是偶数(除以2余数为0)。 - 如果是偶数,则添加到数组中。 3. 当找到一个偶数后,检查数组的长度。如果到达或超过5个元素,就开始一个新的行并清空数组,只保留最后一个偶数。 4. 循环结束后,无论数组是否满5个元素,都需要打印剩余的偶数。 这是一个伪代码示例: ```python array = [] num = 100 whil
recommend-type

文本动画新体验:textillate插件功能介绍

资源摘要信息:"textillate是一个文本动画插件,主要应用于前端开发中,为网页上的文字元素添加动态的动画效果。textillate插件通常用于实现文字的淡入淡出、滑动、旋转等多种动画效果,增强用户界面的交互体验。该插件支持HTML5和CSS3,可以轻松集成到现有的网页项目中,无需复杂配置即可使用。textillate插件主要面向设计师和前端开发者,通过简单的JavaScript代码,即可实现复杂且美观的文本动画效果。 textillate插件的使用通常结合了jQuery库,因此在使用前需要确保页面已经加载了jQuery。该插件的安装包中可能包含有多种文件,如JavaScript、CSS和可能的文档说明等。安装后,开发者需要在HTML文档中引入相关文件,然后通过简单的API调用来应用动画效果到指定的文本元素。 textillate插件的核心功能包括但不限于: 1. 支持多种动画类型:包括但不限于淡入淡出、缩放、摇动、翻转、旋转等。 2. 自定义动画序列:开发者可以定义动画的播放顺序,以及每个动画的持续时间等参数。 3. 支持链式调用:可以将多个动画效果串联起来,实现复合的动画序列。 4. 响应式设计支持:确保在不同屏幕尺寸和分辨率的设备上都能保持良好的动画效果。 5. 轻量级性能:虽然提供了丰富的动画效果,但textillate插件经过优化,不会对页面的加载和性能产生明显影响。 开发者在使用textillate时,应阅读相关文档,了解如何正确引用插件文件、如何初始化插件以及如何自定义动画参数等。此外,了解与textillate配合使用的其他前端技术,如CSS预处理器、JavaScript模块加载器等,也有助于更好地发挥textillate插件的性能。 textillate作为一款流行的文本动画插件,非常适合用于制作网页上的宣传动画、产品介绍动画或者页面过渡动画。通过为文本元素添加动画,可以使网页更加生动有趣,提升用户体验。同时,插件提供的灵活性允许设计师和开发者创造出个性化的动画效果,满足不同项目的需求。 在选择使用textillate插件时,开发者应当注意兼容性问题,确保所使用的浏览器版本支持HTML5和CSS3的相关特性。此外,对于动画效果较多或复杂的页面,还应该评估动画对页面性能的影响,避免影响用户访问速度和页面交互流畅度。"