MATLAB中AR模型功率谱估计中AR阶次估计的实现-psd_my.rar(最近看了几个关于功率谱的问题,有关AR模型的谱估计,在此分享一下,希望大家不吝指正)(声明:本文内容摘自我的毕业论文——心率变异信号的预处理及功率谱估计)(按:AR模型功率谱估计是对非平稳随机信号功率谱估计的常用方法,但是其模型阶次的估计,除了HOSA工具箱里的arorder函数外,没有现成的函数可用,arorder函数是基于矩阵SVD分解的阶次估计方法,为了比较各种阶次估计方法的区别,下面的函数使用了'FPE','AIC','MDL','CAT'集中准则一并估计,并采用试验方法确定那一个阶次更好。
)………………………………以上省略……………………………………………………………………假设原始数据序列为x,那么n阶参数使用最小二乘估计在MATLAB中实现如下:Y=x;Y(1:n)=[];m=N-n;X=[];%构造系数矩阵fori=1:m  forj=1:n      X(i,j)=xt(ni-j);  endendbeta=inv(X'*X)*X'*Y';复制代码beta即为用最小二乘法估计出的模型参数。
此外,还有估计AR模型参数的Yule-Walker方程法、基于线性预测理论的Burg算法和修正的协方差算法等[26]。
相应的参数估计方法在MATLAB中都有现成的函数,比如aryule、arburg以及arcov等。
4.3.3AR模型阶次的选择及实验设计文献[26]中介绍了五种不同的AR模型定阶准则,分别为矩阵奇异值分解(SingularValueDecomposition,SVD)定阶法、最小预测定误差阶准则(FinalPredictionErrorCriterion,FPE)、AIC定阶准则(Akaika’sInformationtheoreticCriterion,AIC)、MDL定阶准则以及CAT定阶准则。
文献[28]中还介绍了一种BIC定阶准则。
SVD方法是对Yule-Walker方程中的自相关矩阵进行SVD分解来实现的,在MATLAB工具箱中arorder函数就是使用的该算法。
其他五种算法的基本思想都是建立目标函数,阶次估计的标准是使目标函数最小化。
以上定阶准则在MATLAB中也可以方便的实现,下面是本文实现FPE、AIC、MDL、CAT定阶准则的程序(部分):form=1:N-1  ……    %判断是否达到所选定阶准则的要求  ifstrcmp(criterion,'FPE')    objectfun(m1)=(N(m1))/(N-(m1))*E(m1);  elseifstrcmp(criterion,'AIC')    objectfun(m1)=N*log(E(m1))2*(m1);  elseifstrcmp(criterion,'MDL')    objectfun(m1)=N*log(E(m1))(m1)*log(N);  elseifstrcmp(criterion,'CAT')    forindex=1:m1        temp=temp(N-index)/(N*E(index));    end    objectfun(m1)=1/N*temp-(N-(m1))/(N*E(m1));  end    ifobjectfun(m1)>=objectfun(m)    orderpredict=m;    break;  endend复制代码orderpredict变量即为使用相应准则预测的AR模型阶次。
(注:以上代码为结合MATLAB工具箱函数pburg,arburg两个功率谱估计函数增加而得,修改后的pburg等函数会在附件中示意,名为pburgwithcriterion)登录/注册后可看大图程序1.JPG(35.14KB,下载次数:20352)下载附件 保存到相册2009-8-2820:54上传登录/注册后可看大图程序2.JPG(51.78KB,下载次数:15377)下载附件 保存到相册2009-8-2820:54上传下面本文使用3.2.1实验设计的输出结果即20例经预处理的HRV信号序列作为实验对象,分别使用FPE、AIC、MAL和CAT定阶准则预测AR模型阶次,图4.1(见下页)为其中一例典型信号使用不同预测准则其目标函数随阶次的变化情况。
从图中可以看出,使用FPE、AIC以及MDL定阶准则所预测的AR模型阶次大概位于10附近,即阶次10左右会使相应的目标函数最小化,符合定阶准则的要求,使用CAT定阶准则预测的阶次较小,在5~10之间。
图4.2(见下页)为另一例信号的阶次估计情况,从中也可以得到同样的结论。
(注,实验信号为实验室所得,没有上传)登录/注册后可看大图图片1.JPG(28.68KB,下载次数:5674)下载附件 保存到相册2009-8-2820:54上传
2025/6/27 16:08:25 6KB matlab
1
使用方法:直接复制到dev-c++目录下的Templates文件夹,新建工程之后就可以看到有Qt的选项.该资源是模仿hangxin1940上传的资源写的,不过增加了几个选项,并且让其只能创建CPP文件.
2025/6/27 7:10:48 32KB Dev-cpp QT Templates
1
纯html+JS完成的抽奖系统,友通数码港实力派团队开发。
1、号码抽取等级、数量可以自由定义前端定义,不用手工修改任何代码3、背景、主题和样式可根据会场需要自由更换修改2、全屏幕显示1024X768像素,界面美观大方,适合于投影仪投射晚会晚宴豪华场所抽奖系统。
3、可自由选择键盘及鼠标双重操作方式,其它无用键自动屏蔽。
4、环境要求低,只需要IE浏览器即可完成所有执行和配置工作。
5、未到场人员可以双击数字单独重新抽取。
6、开完配置一次抽一组或者抽一个。
7、自定义配置可以根据需要保存为配置文件。
8、抽奖完成后可以选择保存结果以便日后统计分析。
V1.0增加了键鼠屏蔽参数-----------------------------------------------V1.12010-12-28增加了是否可以重复抽取配置选项-----------------------------------------------V1.22011-01-05增加了每组是否逐个抽取配置选项2011-01-05对逐个抽取的号码显示做了美化-----------------------------------------------V1.2.12011-01-08修复了window.obo=1;时出现重复数字的错误,V1.1之前版本无此错误~感谢石路街道的领导提醒:)2011-01-08增加了可选滚动姓名的功能配置m_name=[];即可2011-01-08增加双击数字单独重抽时的提示选项“此号码在以后的环节是否还有机会抽到!”2011-01-08修改了结果显示部分的错位样式2011-01-08修复了组抽号码重复错误问题-----------------------------------------------V1.2.22011-01-14兼容了火狐浏览器,但本人建议仍然用IE,要用火狐只有按F11手动全屏了2011-01-14部分键F5、退格、Ctrl+R、Ctrl+N、Shift+F10、Alt+F4被屏蔽,避免一些意外发生2011-01-14优化了一些繁琐的结构,效率有所提升!2011-01-16去掉了配置参数r、r_name,奖项配置变为window.ini,使配合样式表更加灵活-----------------------------------------------V1.2.22011-02-10在配置文件增加了速度控制变量,以便于在不同机器环境中做速度微调-----------------------------------------------V1.2.3*2011-08-08本版本为政府单位定制版本不对外开放-----------------------------------------------V2.02011-09-13突破性的解决了手工修改代码烦恼增加了图形化界面2011-09-13可以通过图形化设置界面保存配置文件(只在IE下有效)2011-09-13配置文件从HTML内提取出来,变为ini.js2011-09-13图形化配置界面做了一些兼容性优化2011-09-13奖项框除了可以填写样式名还可以直接填写样式代码,程序可以自动识别2011-09-13增加了手动和自动整理列表功能-----------------------------------------------V2.0.12011-12-09bug处理:列表获取处理错误,已修正!使用请看压缩包内说明,用过了别忘了给好评,欢迎光临苏州友通数码港!
1
webAPI示例,webAPI增加修改删除
2025/6/25 21:12:58 936KB webApi
1
MODBUS是OSI模型第7层上的应用层报文传输协议,它在连接至不同类型总线或网络的设备之间提供客户机/服务器通信。
自从1979年出现工业串行链路的事实标准以来,MODBUS使成千上万的自动化设备能够通信。
目前,继续增加对简单而雅观的MODBUS结构支持。
互联网组织能够使TCP/IP栈上的保留系统端口502访问MODBUS。
MODBUS是一个请求/应答协议,并且提供功能码规定的服务。
MODBUS功能码是MODBUS请求/应答PDU的元素。
本文件的作用是描述MODBUS事务处理框架内使用的功能码。
2025/6/25 17:32:53 1.05MB modbus RTU modbus TCP
1
定义一个存折类CAccount,存折类具有帐号、姓名、余额等数据成员,可以实现存款、取款和查询余额的操作,取款金额必须在余额范围内。
从存折类派生出信用卡类CCreditcard,信用卡类增加了透支限额一项数据成员,对存款操作进行修改,允许在限额范围内透支金额。
在main函数中建立这两个类的对象并测试之。
2025/6/25 13:15:50 995B 存折存款取款
1
在网上找了很多十六进制数据转ASCII码的工具,都不十分满意。
于是自己写了个小工具,分享给大家,欢迎指正。
(该工具不适用于unicode转换,后续可能会增加该功能)
2025/6/24 6:26:50 635KB ascii 十六进制 转换 工具
1
此固件为华为8145v和8145c专用使用tftp方式和WEB方式刷机使用WEB网页直接刷机方式(推荐),无需任何工具(本人所发华为固件全部支持WEB直刷,并且已经自带补全shell,包括以前的所发的固件)或tftp刷机方式命令:loadpackbytftpsvrip192.168.1.2remotefile8145V-R018CPC110-ff-v20190509.bin更新:(以下增加全部为电信或其它运营商界面中的功能,主要以电信为主)各更固件更新略有差异,具体以实际为准1,增加电信(或其它运营商)界面,WEB刷机方式,包括普通帐号登陆也有此功能管理--->设备管理--->固件升级WEB刷机轻松搞定以后再也不用去满世界找刷机工具了2,增加指示灯状态管理3,QOS增加上下行控制管理4,增加第三方插件管理5,增加全局无线高级设置6,等等一些其它小更新刷机有风险,风险自己担。
由于各运营商支持版本不相同,请选择适合自己的版本v20190509加强版列表:8145v017118加强v20190509版8145v018100加强v20190509版8145v018108加强v20190509版8145v018110加强v20190509版
2025/6/23 18:50:52 206.11MB 8145v
1
这个是很经典的问题实验题目:生产者与消费者(综合性实验)实验环境:C语言编译器实验内容:①由用户指定要产生的进程及其类别,存入进入就绪队列。
  ②调度程序从就绪队列中提取一个就绪进程运行。
如果申请的资源被阻塞则进入相应的等待队列,调度程序调度就绪队列中的下一个进程。
进程运行结束时,会检查对应的等待队列,激活队列中的进程进入就绪队列。
运行结束的进程进入over链表。
重复这一过程直至就绪队列为空。
  ③程序询问是否要继续?如果要转直①开始执行,否则退出程序。
实验目的:通过实验模拟生产者与消费者之间的关系,了解并掌握他们之间的关系及其原理。
由此增加对进程同步的问题的了解。
实验要求:每个进程有一个进程控制块(PCB)表示。
进程控制块可以包含如下信息:进程类型标号、进程系统号、进程状态、进程产品(字符)、进程链指针等等。
系统开辟了一个缓冲区,大小由buffersize指定。
程序中有三个链队列,一个链表。
一个就绪队列(ready),两个等待队列:生产者等待队列(producer);
消费者队列(consumer)。
一个链表(over),用于收集已经运行结束的进程本程序通过函数模拟信号量的操作。
参考书目:1)徐甲同等编,计算机操作系统教程,西安电子科技大学出版社2)AndrewS.Tanenbaum著,陈向群,马红兵译.现代操作系统(第2版).机械工业出版社3)AbranhamSilberschatz,PeterBaerGalvin,GregGagne著.郑扣根译.操作系统概念(第2版).高等教育出版社4)张尧学编著.计算机操作系统教程(第2版)习题解答与实验指导.清华大学出版社实验报告要求:(1)每位同学交一份电子版本的实验报告,上传到202.204.125.21服务器中。
(2)文件名格式为班级、学号加上个人姓名,例如:电子04-1-040824101**.doc  表示电子04-1班学号为040824101号的**同学的实验报告。
(3)实验报告内容的开始处要列出实验的目的,实验环境、实验内容等的说明,报告中要附上程序代码,并对实验过程进行说明。
基本数据结构:PCB*readyhead=NULL,*readytail=NULL;//就绪队列PCB*consumerhead=NULL,*consumertail=NULL;//消费者队列PCB*producerhead=NULL,*producertail=NULL;//生产者队列over=(PCB*)malloc(sizeof(PCB));//over链表intproductnum=0;//产品数量intfull=0,empty=buffersize;//semaphorecharbuffer[buffersize];//缓冲区intbufferpoint=0;//缓冲区指针structpcb{/*定义进程控制块PCB*/intflag;//flag=1denoteproducer;flag=2denoteconsumer;intnumlabel;charproduct;charstate;structpcb*processlink;……};processproc()---给PCB分配内存。
产生相应的的进程:输入1为生产者进程;
输入2为消费者进程,并把这些进程放入就绪队列中。
waitempty()---如果缓冲区满,该进程进入生产者等待队列;
linkqueue(exe,&producertail);//把就绪队列里的进程放入生产者队列的尾部voidsignalempty()boolwaitfull()voidsignalfull()voidproducerrun()voidcomsuerrun()voidmain(){processproc();element=hasElement(readyhead);while(element){exe=getq(readyhead,&readytail);printf("进程%d申请运行,它是一个",exe->numlabel);exe->flag==1?printf("生产者\n"):printf("消费者\n");if(exe->flag==1)producerrun();elsecomsuerrun();element=hasElement(readyhead);}printf("就绪队列没有进程\n");if(ha
1

微信小程序是一种轻量级的应用开发平台,主要针对移动端,尤其是微信用户。
"微票"是一个典型的微信小程序示例,主要用于票务管理或者活动报名等场景。
这个源码可以帮助开发者了解和学习微信小程序的开发流程、架构设计以及功能实现。
1. **微信小程序的基本概念** 微信小程序是腾讯公司推出的一种不需要下载安装即可使用的应用,它实现了“触手可及”的梦想,用户扫一扫或搜索即可打开应用。
同时,小程序也体现了“用完即走”的理念,不占手机内存,体验接近原生App。
2. **开发环境与工具** 开发微信小程序,你需要使用微信官方提供的开发者工具,该工具集成了代码编辑、预览、调试和发布等功能。
通过这个工具,你可以编写WXML(微信小程序标记语言)和WXSS(微信小程序样式语言),以及JavaScript来实现业务逻辑。
3. **WXML与WXSS** - **WXML**:类似于HTML,用于结构布局,但它独有的绑定机制可以将数据和视图紧密关联。
- **WXSS**:类似CSS,用于定义样式,但增加了微信小程序特有的选择器和单位。
4. **JavaScript的运用** 在微信小程序中,JavaScript主要负责数据管理和页面逻辑。
通过`Page`和`Component`进行页面和组件的生命周期管理,使用`wx.request`进行网络请求,`wx.setStorage`和`wx.getStorage`处理本地存储等。
5. **微票源码结构分析** "weapp-weipiao-master"可能包含以下结构: - `app.js`:全局配置,包括初始化数据、全局函数等。
- `app.json`:项目配置,定义小程序的页面路径、窗口表现、网络超时时间等。
- `app.wxss`:全局样式表。
- `pages/`:存放各个页面的文件夹,每个页面都有对应的`.wxml`, `.wxss`, `.js`, `.json`文件。
- `utils/`:可能包含一些辅助函数和工具模块。
- `components/`:自定义组件的目录,提高代码复用性。
6. **微票功能解析** "微票"小程序可能包括以下功能: - 用户登录与注册:通过微信授权快速登录。
- 票务展示:显示各类活动或演出的票务信息。
- 票务购买:用户选择座位,支付购票。
- 订单管理:查看、取消、确认订单。
- 活动详情:介绍活动背景、时间、地点等信息。
- 评论与分享:用户对活动进行评价,分享到社交平台。
7. **学习与实践** 通过研究"微票"源码,开发者可以了解如何在微信小程序中实现数据绑定、事件处理、网络请求、页面跳转、动画效果等常见功能,为自己的小程序项目积累经验。
8. **优化与扩展** 对于已有的"微票"源码,开发者可以考虑优化用户体验,如增加推送通知、添加优惠券功能、集成第三方支付等,以满足更多实际需求。
"微信小程序-微票源码"是一个很好的学习资源,对于想要涉足微信小程序开发的开发者来说,通过深入理解并实践这个源码,能够快速掌握小程序开发的核心技术,并能在此基础上创建自己的小程序应用。
2025/6/20 1:32:45 769KB
1
共 1000 条记录 首页 上一页 下一页 尾页
在日常工作中,钉钉打卡成了我生活中不可或缺的一部分。然而,有时候这个看似简单的任务却给我带来了不少烦恼。 每天早晚,我总是得牢记打开钉钉应用,点击"工作台",再找到"考勤打卡"进行签到。有时候因为工作忙碌,会忘记打卡,导致考勤异常,影响当月的工作评价。而且,由于我使用的是苹果手机,有时候系统更新后,钉钉的某些功能会出现异常,使得打卡变得更加麻烦。 另外,我的家人使用的是安卓手机,他们也经常抱怨钉钉打卡的繁琐。尤其是对于那些不太熟悉手机操作的长辈来说,每次打卡都是一次挑战。他们总是担心自己会操作失误,导致打卡失败。 为了解决这些烦恼,我开始思考是否可以通过编写一个全自动化脚本来实现钉钉打卡。经过一段时间的摸索和学习,我终于成功编写出了一个适用于苹果和安卓系统的钉钉打卡脚本。
2024-04-09 15:03 15KB 钉钉 钉钉打卡