for c = 1 : length(opt.C) for n = 1 : length(opt.nu) cur_C = opt.C(c); cur_nu = opt.nu(n); mean_Xbase_rec = zeros(size(Sig_Ybase, 1), opt.pca_d(d)); mean_Xval_rec = zeros(size(Sig_Yval, 1), opt.pca_d(d)); mean_Xval = zeros(size(Sig_Yval, 1), opt.pca_d(d)); parfor j = 1 : opt.pca_d(d) % SVR learning and testing regressor = svmtrain(mean_Xbase(:, j), Ker_base, ['-s 4 -t 4 -c ' num2str(cur_C) ' -n ' num2str(cur_nu) ' -m 10000']); mean_Xbase_rec(:, j) = svmpredict(mean_Xbase(:, j), Ker_base, regressor); mean_Xval_rec(:, j) = svmpredict(mean_Xval(:, j), Ker_val, regressor); end [dis_eu, dis_seu, acc_eu, acc_seu] = test_EXEM(mapped_Xval, Yval, mean_Xval_rec, avg_std_Xbase); val_dis_eu(c, n, g, d) = val_dis_eu(c, n, g, d) + dis_eu / nr_fold; val_dis_seu(c, n, g, d) = val_dis_seu(c, n, g, d) + dis_seu / nr_fold; val_acc_eu(c, n, g, d) = val_acc_eu(c, n, g, d) + acc_eu / nr_fold; val_acc_seu(c, n, g, d) = val_acc_seu(c, n, g, d) + acc_seu / nr_fold; end end

时间: 2024-02-14 18:27:38 浏览: 214
这段代码是用于在EXEM训练中进行SVR(Support Vector Regression)学习和评估的部分。在C和nu的循环中,对每个组合进行以下操作: 1. 设置当前的C和nu参数值为`cur_C`和`cur_nu`。 2. 初始化用于存储重构结果的矩阵`mean_Xbase_rec`和`mean_Xval_rec`,以及用于存储验证数据集平均值的矩阵`mean_Xval`。 3. 并行地对每个PCA维度进行以下操作: a. 使用`svmtrain`函数,基于训练数据集的PCA降维结果`mean_Xbase(:, j)`和训练数据集的核矩阵`Ker_base`进行SVR学习。使用给定的C和nu参数进行模型训练。 b. 使用`svmpredict`函数,基于训练得到的SVR模型,对训练数据集和验证数据集中的PCA降维结果进行预测,得到重构后的结果`mean_Xbase_rec(:, j)`和`mean_Xval_rec(:, j)`。 4. 调用`test_EXEM`函数,对重构后的验证数据集、真实类别标签和平均标准差进行评估。返回值包括欧氏距离(dis_eu)、标准化欧氏距离(dis_seu)、准确率(acc_eu)和标准化准确率(acc_seu)。 5. 将评估结果加权平均到相应的结果矩阵中。`val_dis_eu`、`val_dis_seu`、`val_acc_eu`和`val_acc_seu`是用于存储不同参数组合下的评估结果的多维矩阵。 这段代码的目的是在EXEM训练过程中,对不同参数组合下的SVR模型进行学习和评估,并将评估结果存储在相应的结果矩阵中,以备后续使用。
阅读全文

相关推荐

用C++编写程序,要求如下: ①输入多组数据,总计n*( a+b+2)+1行。其中,第一行整数n代表总计有n组数据,之后依次输入n组数据。每组数据包括a+b+2行,其中第一行是两个整数a和b,分别代表A(x)与B(x)的项数。之后紧跟a行,每行两个整数a1和a2,分别代表A(x)每项的系数和指数,再之后紧跟b行,每行两个整数b1和b2,分别代表B(x)每项的系数和指数,每组数据最后一行为一个字符(+、-、*、'),分别代表多项式的加法、减法、乘法和求导运算。所有数据的绝对值小于100,指数大于等于0。 ②编写的程序在我给出的代码上进行补充 ③当用户输入: 4 1 1 1 0 1 1 + 4 3 7 0 3 1 9 8 5 17 8 1 22 7 -9 8 + 1 1 1 1 1 1 - 1 1 1 1 1 1 ' 输出: 1x^1+1 5x^17+22x^7+11x^1+7 0 1 1 #include <iostream>#include <string> using namespace std; typedef struct LNode{ int coe;int exp;struct LNode *next; }LNode,*LinkList; void CreatePolynomial(LinkList &L,int n){ L=new LNode;L->next=NULL; for(int i=0;i<n;i++){ LinkList p=new LNode;cin>>p->coe>>p->exp; LinkList pre=L,cur=L->next; while(cur&&p->exp<cur->exp){ pre=cur;cur=cur->next;} p->next=cur;pre->next=p;} } void OutputPolynomial(LinkList L){ if(!L||!L->next) cout<<0;LinkList p=L->next; while(p){ if(p==L->next){ if (p->exp!=0) cout<coe<<"x^"<exp; else cout<coe;} else{ if(p->coe>0) cout<<"+"; if(p->exp!=0) cout<coe<<"x^"<exp; else cout<coe;} p=p->next;} cout<<endl;} LinkList Add(LinkList LA,LinkList LB){} void Minus(LinkList LA,LinkList LB){} void Mul(LinkList LA,LinkList LB){} void Diff(LinkList L){ LinkList p=L->next;LinkList r=NULL; while(p){ p->coe*=p->exp;p->exp--; if(p->exp<0){ r=p;p=p->next;delete r;} else{ p=p->next;} } OutputPolynomial(L);} void Opt(LinkList &LA,LinkList &LB,string s){ if(s=="+") OutputPolynomial(Add(LA, LB));if(s=="-") Minus(LA, LB); if(s=="*") Mul(LA, LB);if(s=="'"){ Diff(LA);Diff(LB);} } int main(){ int n;cin>>n; while(n--){ int a,b;cin>>a>>b;LinkList LA,LB;CreatePolynomial(LA,a); CreatePolynomial(LB,b);string s;cin>>s;Opt(LA,LB,s);} return 0;}

请将如下的matlab代码转为python代码,注意使用pytorch框架实现,并对代码做出相应的解释:function [nets,errors]=BPMLL_train(train_data,train_target,hidden_neuron,alpha,epochs,intype,outtype,Cost,min_max) rand('state',sum(100clock)); if(nargin<9) min_max=minmax(train_data'); end if(nargin<8) Cost=0.1; end if(nargin<7) outtype=2; end if(nargin<6) intype=2; end if(nargin<5) epochs=100; end if(nargin<4) alpha=0.05; end if(intype==1) in='logsig'; else in='tansig'; end if(outtype==1) out='logsig'; else out='tansig'; end [num_class,num_training]=size(train_target); [num_training,Dim]=size(train_data); Label=cell(num_training,1); not_Label=cell(num_training,1); Label_size=zeros(1,num_training); for i=1:num_training temp=train_target(:,i); Label_size(1,i)=sum(temp==ones(num_class,1)); for j=1:num_class if(temp(j)==1) Label{i,1}=[Label{i,1},j]; else not_Label{i,1}=[not_Label{i,1},j]; end end end Cost=Cost2; %Initialize multi-label neural network incremental=ceil(rand100); for randpos=1:incremental net=newff(min_max,[hidden_neuron,num_class],{in,out}); end old_goal=realmax; %Training phase for iter=1:epochs disp(strcat('training epochs: ',num2str(iter))); tic; for i=1:num_training net=update_net_ml(net,train_data(i,:)',train_target(:,i),alpha,Cost/num_training,in,out); end cur_goal=0; for i=1:num_training if((Label_size(i)~=0)&(Label_size(i)~=num_class)) output=sim(net,train_data(i,:)'); temp_goal=0; for m=1:Label_size(i) for n=1:(num_class-Label_size(i)) temp_goal=temp_goal+exp(-(output(Label{i,1}(m))-output(not_Label{i,1}(n)))); end end temp_goal=temp_goal/(mn); cur_goal=cur_goal+temp_goal; end end cur_goal=cur_goal+Cost0.5(sum(sum(net.IW{1}.*net.IW{1}))+sum(sum(net.LW{2,1}.*net.LW{2,1}))+sum(net.b{1}.*net.b{1})+sum(net.b{2}.*net.b{2})); disp(strcat('Global error after ',num2str(iter),' epochs is: ',num2str(cur_goal))); old_goal=cur_goal; nets{iter,1}=net; errors{iter,1}=old_goal; toc; end disp('Maximum number of epochs reached, training process completed');

void ifOptLaserTmpGet15minHistory(net_if *pNetIf,UINT8 portNum) { char objectName[MAX_IFNAME_LEN] = {0}; char pmParaName[MAX_COMMON_LEN] = {0}; char granularity[MAX_COMMON_LEN] = {0}; char startTime[MAX_ALARM_TIME_LENGTH] = {0}; char objectType[MAX_COMMON_LEN] = {0}; char maxValue[MAX_COMMON_LEN] = {0}; char minValue[MAX_COMMON_LEN] = {0}; char aveValue[MAX_COMMON_LEN] = {0}; char curValue[MAX_COMMON_LEN] = {0}; char *pTime = NULL; char *pStartTime = NULL; time_t etime; UINT32 length = 0; char timechange[32]={0}; etime = time(NULL); if (etime == -1) { ERROR("ifAnalogGet get start or end time failed (%s).", strerror(errno)); return; } if ((pTime = nc_time2datetime(etime, NULL)) == NULL) { ERROR("ifAnalogGet Internal error when converting time formats."); return; } pStartTime = pTime; timechange15Min(pStartTime,timechange); snprintf(objectName,MAX_IFNAME_LEN,"PTP=/shelf=1/slot=1/subslot=1/port=%u",portNum); snprintf(pmParaName,MAX_COMMON_LEN,"LASER_TMP"); snprintf(granularity,MAX_COMMON_LEN,"%s","15min"); snprintf(startTime,MAX_ALARM_TIME_LENGTH,"%s",timechange); snprintf(objectType,MAX_COMMON_LEN,"%s","PTP"); ncGetPerfInfo(portNum,NC_LASER_TMP,NC_MAX_VALUE,maxValue); ncGetPerfInfo(portNum,NC_LASER_TMP,NC_MIN_VALUE,minValue); ncGetPerfInfo(portNum,NC_LASER_TMP,NC_AVE_VALUE,aveValue); ifGetCurLaserTmp(pNetIf,curValue); //ncGetPerfInfo(portNum - 1,NC_CUR_VALUE,NC_LASER_TMP,curValue); memset(g_15minhistoryperf[portNum-1].parame[Opt_LASERTMP].pkg,0,HISTORY_PERF_PARAM_BUFF); length = snprintf(g_15minhistoryperf[portNum-1].parame[Opt_LASERTMP].pkg, HISTORY_PERF_PARAM_BUFF, PM_MODE_ANALOG_INIT_XML, objectName, pmParaName, granularity, startTime, objectType, maxValue, minValue, aveValue, curValue); DBG("%s/%d xml_send_buffer %s length %d",__func__,__LINE__,g_15minhistoryperf[portNum-1].parame[Opt_LASERTMP].pkg,length); // ncds_apply_pkg_info(pXmlSendBuf,NC_DATASTORE_RUNNING,"acc-performance"); }

zip

大家在看

recommend-type

STM32的FOC库教程

内容如下: 1、STM32_FOC _library_v2.0新功能 2、STM32F103_永磁同步电机_PMSM_FOC软件库_用户手册_中文版 3、STM32F103xx-PMSM-FOC-software-library-UM 4、基于STM32的PMSM FOC软件库(一) 5、基于STM32的PMSM FOC软件库(二) 6、基于STM32的PMSM FOC软件库(三) 7、基于STM32的PMSM FOC软件库(四)
recommend-type

2000-2022年 上市公司-股价崩盘风险相关数据(数据共52234个样本,包含do文件、excel数据和参考文献).zip

上市公司股价崩盘风险是指股价突然大幅下跌的可能性。这种风险可能由多种因素引起,包括公司的财务状况、市场环境、政策变化、投资者情绪等。 测算方式:参考《管理世界》许年行老师和《中国工业经济》吴晓晖老师的做法,使用负收益偏态系数(NCSKEW)和股票收益上下波动比率(DUVOL)度量股价崩盘风险。 数据共52234个样本,包含do文件、excel数据和参考文献。 相关数据指标 stkcd、证券代码、year、NCSKEW、DUVOL、Crash、Ret、Sigma、证券代码、交易周份、周个股交易金额、周个股流通市值、周个股总市值、周交易天数、考虑现金红利再投资的周个股回报率、市场类型、周市场交易总股数、周市场交易总金额、考虑现金红利再投资的周市场回报率(等权平均法)、不考虑现金红利再投资的周市场回报率(等权平均法)、考虑现金红利再投资的周市场回报率(流通市值加权平均法)、不考虑现金红利再投资的周市场回报率(流通市值加权平均法)、考虑现金红利再投资的周市场回报率(总市值加权平均法)、不考虑现金红利再投资的周市场回报率(总市值加权平均法)、计算周市场回报率的有效公司数量、周市场流通市值、周
recommend-type

Mac OS X10.6.3 Snow Leopard系统 中文版完整安装盘 下载地址连接

Mac OS X10.6.3 Snow Leopard系统 中文版完整安装盘 下载链接,速度稳定。 Mac OS X10.6.3 Snow Leopard系统 中文版完整安装盘 下载链接,速度稳定。
recommend-type

SigmaStudioHelp_3.0(中文)

关于DSP 的技术文档,留住入门DSP 控制用作备份;DSP核心技术都在里面了解;
recommend-type

涉密网络建设方案模板.doc

涉密网络建设方案模板.doc

最新推荐

recommend-type

微软内部资料-SQL性能优化2

Lesson 1: Memory 3 Lesson 2: I/O 73 Lesson 3: CPU 111 Module 3: Troubleshooting Server Performance Module Overview Troubleshooting server performance-based support calls requires product...
recommend-type

java计算器源码.zip

java毕业设计源码,可供参考
recommend-type

FRP Manager-V1.19.2

Windows下的FRP图形化客户端,对应FRP版本0.61.1,需要64位操作系统
recommend-type

基于优化EKF的PMSM无位置传感器矢量控制研究_崔鹏龙.pdf

基于优化EKF的PMSM无位置传感器矢量控制研究_崔鹏龙.pdf
recommend-type

旧物置换网站(基于springboot,mysql,java).zip

旧物置换网站的开发过程中,采用B / S架构,主要使用Java技术进行开发,结合最新流行的springboot框架。中间件服务器是Tomcat服务器,使用Mysql数据库和Eclipse开发 环境。该旧物置换网站包括管理员、用户、卖家。其主要功能包括管理员:首页、个人中心、用户管理、卖家管理、旧物类型管理、旧物信息管理、置换交易管理、系统管理等,卖家后台:首页、个人中心、旧物类型管理、旧物信息管理、置换交易管理。前台首页;首页、旧物信息、网站公告、个人中心、后台管理等,用户后台:首页、个人中心、旧物信息管理、置换交易管理、用户可根据关键字进行信息的查找自己心仪的信息等。 (1)用户功能需求 用户进入前台系统可以查看首页、旧物信息、网站公告、个人中心、后台管理等操作。前台首页用例如图3-1所示。 (2)管理员功能需求 管理员登陆后,主要功能模块包括首页、个人中心、用户管理、卖家管理、旧物类型管理、旧物信息管理、置换交易管理、系统管理等功能。 关键词:旧物置换网站,Mysql数据库,Java技术 springboot框架
recommend-type

PHP集成Autoprefixer让CSS自动添加供应商前缀

标题和描述中提到的知识点主要包括:Autoprefixer、CSS预处理器、Node.js 应用程序、PHP 集成以及开源。 首先,让我们来详细解析 Autoprefixer。 Autoprefixer 是一个流行的 CSS 预处理器工具,它能够自动将 CSS3 属性添加浏览器特定的前缀。开发者在编写样式表时,不再需要手动添加如 -webkit-, -moz-, -ms- 等前缀,因为 Autoprefixer 能够根据各种浏览器的使用情况以及官方的浏览器版本兼容性数据来添加相应的前缀。这样可以大大减少开发和维护的工作量,并保证样式在不同浏览器中的一致性。 Autoprefixer 的核心功能是读取 CSS 并分析 CSS 规则,找到需要添加前缀的属性。它依赖于浏览器的兼容性数据,这一数据通常来源于 Can I Use 网站。开发者可以通过配置文件来指定哪些浏览器版本需要支持,Autoprefixer 就会自动添加这些浏览器的前缀。 接下来,我们看看 PHP 与 Node.js 应用程序的集成。 Node.js 是一个基于 Chrome V8 引擎的 JavaScript 运行时环境,它使得 JavaScript 可以在服务器端运行。Node.js 的主要特点是高性能、异步事件驱动的架构,这使得它非常适合处理高并发的网络应用,比如实时通讯应用和 Web 应用。 而 PHP 是一种广泛用于服务器端编程的脚本语言,它的优势在于简单易学,且与 HTML 集成度高,非常适合快速开发动态网站和网页应用。 在一些项目中,开发者可能会根据需求,希望把 Node.js 和 PHP 集成在一起使用。比如,可能使用 Node.js 处理某些实时或者异步任务,同时又依赖 PHP 来处理后端的业务逻辑。要实现这种集成,通常需要借助一些工具或者中间件来桥接两者之间的通信。 在这个标题中提到的 "autoprefixer-php",可能是一个 PHP 库或工具,它的作用是把 Autoprefixer 功能集成到 PHP 环境中,从而使得在使用 PHP 开发的 Node.js 应用程序时,能够利用 Autoprefixer 自动处理 CSS 前缀的功能。 关于开源,它指的是一个项目或软件的源代码是开放的,允许任何个人或组织查看、修改和分发原始代码。开源项目的好处在于社区可以一起参与项目的改进和维护,这样可以加速创新和解决问题的速度,也有助于提高软件的可靠性和安全性。开源项目通常遵循特定的开源许可证,比如 MIT 许可证、GNU 通用公共许可证等。 最后,我们看到提到的文件名称 "autoprefixer-php-master"。这个文件名表明,该压缩包可能包含一个 PHP 项目或库的主分支的源代码。"master" 通常是源代码管理系统(如 Git)中默认的主要分支名称,它代表项目的稳定版本或开发的主线。 综上所述,我们可以得知,这个 "autoprefixer-php" 工具允许开发者在 PHP 环境中使用 Node.js 的 Autoprefixer 功能,自动为 CSS 规则添加浏览器特定的前缀,从而使得开发者可以更专注于内容的编写而不必担心浏览器兼容性问题。
recommend-type

揭秘数字音频编码的奥秘:非均匀量化A律13折线的全面解析

# 摘要 数字音频编码技术是现代音频处理和传输的基础,本文首先介绍数字音频编码的基础知识,然后深入探讨非均匀量化技术,特别是A律压缩技术的原理与实现。通过A律13折线模型的理论分析和实际应用,本文阐述了其在保证音频信号质量的同时,如何有效地降低数据传输和存储需求。此外,本文还对A律13折线的优化策略和未来发展趋势进行了展望,包括误差控制、算法健壮性的提升,以及与新兴音频技术融合的可能性。 # 关键字 数字音频编码;非均匀量化;A律压缩;13折线模型;编码与解码;音频信号质量优化 参考资源链接:[模拟信号数字化:A律13折线非均匀量化解析](https://wenku.csdn.net/do
recommend-type

arduino PAJ7620U2

### Arduino PAJ7620U2 手势传感器 教程 #### 示例代码与连接方法 对于Arduino开发PAJ7620U2手势识别传感器而言,在Arduino IDE中的项目—加载库—库管理里找到Paj7620并下载安装,完成后能在示例里找到“Gesture PAJ7620”,其中含有两个示例脚本分别用于9种和15种手势检测[^1]。 关于连线部分,仅需连接四根线至Arduino UNO开发板上的对应位置即可实现基本功能。具体来说,这四条线路分别为电源正极(VCC),接地(GND),串行时钟(SCL)以及串行数据(SDA)[^1]。 以下是基于上述描述的一个简单实例程序展示如
recommend-type

网站啄木鸟:深入分析SQL注入工具的效率与限制

网站啄木鸟是一个指的是一类可以自动扫描网站漏洞的软件工具。在这个文件提供的描述中,提到了网站啄木鸟在发现注入漏洞方面的功能,特别是在SQL注入方面。SQL注入是一种常见的攻击技术,攻击者通过在Web表单输入或直接在URL中输入恶意的SQL语句,来欺骗服务器执行非法的SQL命令。其主要目的是绕过认证,获取未授权的数据库访问权限,或者操纵数据库中的数据。 在这个文件中,所描述的网站啄木鸟工具在进行SQL注入攻击时,构造的攻击载荷是十分基础的,例如 "and 1=1--" 和 "and 1>1--" 等。这说明它的攻击能力可能相对有限。"and 1=1--" 是一个典型的SQL注入载荷示例,通过在查询语句的末尾添加这个表达式,如果服务器没有对SQL注入攻击进行适当的防护,这个表达式将导致查询返回真值,从而使得原本条件为假的查询条件变为真,攻击者便可以绕过安全检查。类似地,"and 1>1--" 则会检查其后的语句是否为假,如果查询条件为假,则后面的SQL代码执行时会被忽略,从而达到注入的目的。 描述中还提到网站啄木鸟在发现漏洞后,利用查询MS-sql和Oracle的user table来获取用户表名的能力不强。这表明该工具可能无法有效地探测数据库的结构信息或敏感数据,从而对数据库进行进一步的攻击。 关于实际测试结果的描述中,列出了8个不同的URL,它们是针对几个不同的Web应用漏洞扫描工具(Sqlmap、网站啄木鸟、SqliX)进行测试的结果。这些结果表明,针对提供的URL,Sqlmap和SqliX能够发现注入漏洞,而网站啄木鸟在多数情况下无法识别漏洞,这可能意味着它在漏洞检测的准确性和深度上不如其他工具。例如,Sqlmap在针对 "http://www.2cto.com/news.php?id=92" 和 "http://www.2cto.com/article.asp?ID=102&title=Fast food marketing for children is on the rise" 的URL上均能发现SQL注入漏洞,而网站啄木鸟则没有成功。这可能意味着网站啄木鸟的检测逻辑较为简单,对复杂或隐蔽的注入漏洞识别能力不足。 从这个描述中,我们也可以了解到,在Web安全测试中,工具的多样性选择是十分重要的。不同的安全工具可能对不同的漏洞和环境有不同的探测能力,因此在实际的漏洞扫描过程中,安全测试人员需要选择合适的工具组合,以尽可能地全面地检测出应用中存在的漏洞。 在标签中指明了这是关于“sql注入”的知识,这表明了文件主题的核心所在。SQL注入是一种常见的网络攻击方式,安全测试人员、开发人员和网络管理员都需要对此有所了解,以便进行有效的防御和检测。 最后,提到了压缩包子文件的文件名称列表,其中包含了三个文件:setup.exe、MD5.exe、说明_Readme.html。这里提供的信息有限,但可以推断setup.exe可能是一个安装程序,MD5.exe可能是一个计算文件MD5散列值的工具,而说明_Readme.html通常包含的是软件的使用说明或者版本信息等。这些文件名暗示了在进行网站安全测试时,可能涉及到安装相关的软件工具,以及进行文件的校验和阅读相应的使用说明。然而,这些内容与文件主要描述的web安全漏洞检测主题不是直接相关的。
recommend-type

【GPStoolbox使用技巧大全】:20个实用技巧助你精通GPS数据处理

# 摘要 GPStoolbox是一个广泛应用于GPS数据处理的软件工具箱,它提供了从数据导入、预处理、基本分析到高级应用和自动化脚本编写的全套功能。本文介绍了GPStoolbox的基本概况、安装流程以及核心功能,探讨了如何