源代码就是源程序
源代码，是指未编译的文本代码。是一系列人类可读的计算机语言指令。
在现代程序语言中，源代码可以是以书籍或者磁带的形式出现，但最为常用的
格式是文本文件，这种典型格式的目的是为了编译出计算机程序。计算机源代
码的最终目的是将人类可读的文本翻译成为计算机可以执行的二进制指令，这
种过程叫做编译，通过编译器完成。
作用
源代码主要功用有如下2种作用：
生成目标代码，即计算机可以识别的代码。
对软件进行说明，即对软件的编写进行说明。为数不少的初学者，甚至少数有
经验的程序员都忽视软件说明的编写，因为这部分虽然不会在生成的程序中直
接显示，也不参与编译。但是说明对软件的学习、分享、维护和软件复用都有
巨大的好处。因此，书写软件说明在业界被认为是能创造优秀程序的良好习惯，一些公司也硬性规定必须书写。
需要指出的是，源代码的修改不能改变已经生成的目标代码。如果需要目标代
码做出相应的修改，必须重新编译。
代码组合
源代码作为软件的特殊部分，可能被包含在一个或多个文件中。一个程序不必
用同一种格式的源代码书写。例如，一个程序如果有C语言库的支持，那么就
可以用C语言；而另一部分为了达到比较高的运行效率，则可以用汇编语言编写。
较为复杂的软件，一般需要数十种甚至上百种的源代码的参与。为了降低种复
杂度，必须引入一种可以描述各个源代码之间联系，并且如何正确编译的系统。在这样的背景下，修订控制系统（RCS）诞生了，并成为研发者对代码修订的必备工具之一。
还有另外一种组合：源代码的编写和编译分别在不同的平台上实现，专业术语
叫做软件移植。
版权
如果按照源代码类型区分软件，通常被分为两类：自由软件和非自由软件。自
由软件一般是不仅可以免费得到，而且公开源代码；相对应地，非自由软件则
是不公开源代码。所有一切通过非正常手段获得非自由软件源代码的行为都将
被视为非法。
质量
对于计算机而言，并不存在真正意义上的“好”的源代码；然而作为一个人，
好的书写习惯将决定源代码的好坏。源代码是否具有可读性，成为好坏的重要
标准。软件文档则是表明可读性的关键。
效率
虽然我们可以通过不同的语言来实现计算机的同一功能，但在执行效率上则存
在不同。普遍规律是：越高级的语言，其执行效率越低。这也是为什么汇编语言生成的文件比用VB语言生成文件普遍要小的原因。
系统执行软件源程序代 码 Company Document number：WUUT-WUUY-WBBGB-BWYTT-1982GT
附录一系统执行软件源程序代码一、系统多机通信从机源程序代码 ** *Lu Bo* ------------------------------------------------------------------------------------------------------------------------------------- ORG 0000H SJMP MAIN ORG 0023H ；串行口中断入口 SJMP SJCS ------------------主程序---------------------------------------------------------------------------------------------------------- MAIN: SLAVE EQU 12 DK EQU MOV SCON,#0B0H MOV TMOD,#20H MOV TL1,#0FDH MOV TH1,#0FDH MOV IE,#90H SETB TR1 LCALL WDCJ SETB RS1 CLR RS0 LCALL BCDZH LCALL DISPLAY CLR RS1 CLR RS0 SJMP MAIN ------------------串行口中断服务子程序---------------------------------------------------------------------------------- SJCS: CLR RI PUSH A PUSH PSW MOV A,SBUF XRL A,#SLAVE JZ SJCS0 SJCS1: SETB SM2 POP PSW POP A RETI
1.硬币模拟试验 源代码： clear; clc; head_count=0; p1_hist= [0]; p2_hist= [0]; n = 1000; p1 = 0.3; p2=0.03; head = figure(1); rand('seed',sum(100*clock)); fori = 1:n tmp = rand(1); if(tmp<= p1) head_count = head_count + 1; end p1_hist (i) = head_count /i; end figure(head); subplot(2,1,1); plot(p1_hist); grid on; hold on; xlabel('重复试验次数'); ylabel('正面向上的比率'); title('p=0.3试验次数N与正面向上比率的函数图'); head_count=0; fori = 1:n tmp = rand(1); if(tmp<= p2) head_count = head_count + 1; end p2_hist (i) = head_count /i; end figure(head); subplot(2,1,2); plot(p2_hist); grid on; hold on; xlabel('重复试验次数'); ylabel('正面向上的比率'); title('p=0.03试验次数N与正面向上比率的函数图'); 实验结果：
2.不同次数的随机试验均值方差比较 源代码： clear ; clc; close; rand('seed',sum(100*clock)); Titles = ['n=5时' 'n=20时' 'n=25时' 'n=50时' 'n=100时']; Titlestr = cellstr(Titles); X_n_bar=[0]; %the samples of the X_n_bar X_n=[0]; %the samples of X_n N=[5,10,25,50,100]; j=1; num_X_n = 100; num_X_n_bar = 100; h_X_n_bar = figure(1);
C语言程序设计销售管理 系统 The latest revision on November 22, 2020
销售信息管理系统设计文档 组号：第二组 成员：汪品高齐瑞智欧阳鑫赵飞锋章森 目录 一．题目与要求 二．功能实现设计（模块化设计）三．运行结果（附带抓图）四．源程序
报告正文 1、题目与要求： 2、某公司有四个销售员（编号：1-4），负责销售五种产品（编号：1-5）。每个销售员汇总当天出售的每种产品的销售情况。 2 每天销售情况包括（销售员的代号、产品的代号、这种产品当天的销售额） 1、编写一个销售信息管理系统，读取上个月的销售情况 1、计算上个月没人每种产品的销售额 2、按销售额对销售员进行排序，输出排序结果（销售员代号） 3、统计每种产品的总销售额，对这些产品从高到低的顺序输出结果 （须输出产品的代号和销售额） 4、输出统计报表 2、需求分析 本程序将用到文件系统，其中单行内容代表一张纸条。程序运行后首先将文件系统的原始数据读取并保存到程序的结构体数组当中，所以应提供文件输入的操作；由于纸条数据零散，所以要对原始数据进行分类操作并提供一个数组用于保存分类信息；在程序中要进行统计工作，所以要提供显示、排序等操作；另外应提供键盘式选择菜单实现功能选择。 3、功能实现设计 总体设计
根据上面的分析，可以将这个系统分为如下五大模块：查看每人每种产品的销售额、按销售额对销售员排序、对每种产品总销售额从高到低排序、输出统计报表和退出。其中主要模块是前四个模块。 流程图 详细设计: #include<> #include<> #include<> #include<> #define M 【功能函数声明】 vvoid detail(); void foreach(); void menu(); void product(); void readdata(); void seller(); void total(); void choice(); 【数据结构】 原始数据有销售员代号、产品代号和销售额，所以定义一个全局结构体数组scrip用于保存原始数据。然后定义全局变量n保存纸条数目，定义数组s用于保存分类后的信息。 struct ss { int sid; id,&ss[n].pid,&ss[n].sn); } id)
完整代码： 1.客户类（class person）： package客户管理系统; public class Person { private int num; private String name; private String address; private String sex; public Person(){} public Person(int num,String name,String address,String sex){ this.num=num; https://www.doczj.com/doc/6c4633945.html,=name; this.address=address; this.sex=sex; } public void setNum(int num){ this.num=num; } public int getNum(){ return num; } public void setName(String name){ https://www.doczj.com/doc/6c4633945.html,=name; } public String getName(){ return name; } public void setAddress(String address){ this.address=address; } public String getAddress(){ return address;
} public void setSex(String sex){ this.sex=sex; } public String getSex(){ return sex; } public String toString() { return"\t"+num+"\t"+name+"\t"+address+"\t"+sex; } } 2.菜单类(class Menu)： package客户管理系统; import java.util.Scanner; public class Menu { Scanner input=new Scanner(System.in); public Menu() {} public int showMenu() { System.out.println("***********************************" ); System.out.println("客户管理系统"); System.out.println("===================================" ); System.out.println("【1】增加"); System.out.println("【2】删除"); System.out.println("【3】修改"); System.out.println("【4】查询"); System.out.println("【5】浏览"); System.out.println("【6】退出"); System.out.println("***********************************" ); System.out.print("请选择操作：");
matlab语音识别系统(源代码)最新版
目录 一、设计任务及要求 (1) 二、语音识别的简单介绍 2.1语者识别的概念 (2) 2.2特征参数的提取 (3) 2.3用矢量量化聚类法生成码本 (3) 2.4VQ的说话人识别 (4) 三、算法程序分析 3.1函数关系 (4) 3.2代码说明 (5) 3.2.1函数mfcc (5) 3.2.2函数disteu (5) 3.2.3函数vqlbg (6) 3.2.4函数test (6) 3.2.5函数testDB (7) 3.2.6 函数train (8) 3.2.7函数melfb (8) 四、演示分析 (9) 五、心得体会 (11) 附：GUI程序代码 (12)
一、设计任务及要求 用MATLAB实现简单的语音识别功能； 具体设计要求如下： 用MATLAB实现简单的数字1~9的语音识别功能。 二、语音识别的简单介绍 基于VQ的说话人识别系统，矢量量化起着双重作用。在训练阶段，把每一个说话者所提取的特征参数进行分类，产生不同码字所组成的码本。在识别(匹配)阶段，我们用VQ方法计算平均失真测度(本系统在计算距离d时，采用欧氏距离测度)，从而判断说话人是谁。 语音识别系统结构框图如图1所示。 图1 语音识别系统结构框图 2.1语者识别的概念 语者识别就是根据说话人的语音信号来判别说话人的身份。语音是人的自然属性之一，由于说话人发音器官的生理差异以及后天形成的行为差异，每个人的语音都带有强烈的个人色彩，这就使得通过分析语音信号来识别说话人成为可能。用语音来鉴别说话人的身份有着许多独特的优点，如语音是人的固有的特征，不会丢失或遗忘；语音信号的采集方便，系统设备成本低；利用电话网络还可实现远程客户服务等。因此，近几年来，说话人识别越来越多的受到人们的重视。与其他生物识别技术如指纹识别、手形识别等相比较，说话人识别不仅使用方便，而且属于非接触性，容易被用户接受，并且在已有的各种生物特征识别技术中，是唯一可以用作远程验证的识别技术。因此，说话人识别的应用前景非常广泛：今天，说话人识别技术已经关系到多学科的研究领域，不同领域中的进步都对说话人识别的发展做出了贡献。说话人识别技术是集声学、语言学、计算机、信息处理和人工智能等诸多领域的一项综合技术，应用需求将十分广阔。在吃力语音信号的时候如何提取信号中关键的成分尤为重要。语音信号的特征参数的好坏直接导致了辨别的准确性。
有一个小型超市，出售N（N>＝10）种商品，设计并实现一个 系统，完成下列功能： 1．保存及输出。超市中的各种商品信息保存在指定文件中，可以把它们输出显示。 2．计算并排序。计算每类商品的总价值（sum，单精度）及平均价（aver，单精度，输出一位小数），将每类商品按平均价从大到小的顺序排序打印出来。 3．统计。统计输出库存量低于100的货号及类别。统计输出有两种以上（含两种）商品库存量低于100的商品类别。 1.2总体结构 本程序主要分为八个模块：主模块、信息输出修改模块、新建信息模块、排序模块、计算模块、统计模 1、统计模 2、打印模块。 1）主模块：通过调用各分模块实现功能； 2）信息输出修改模块：输出显示文件中商品信息内容，添加商品信息，删除商品信息，修改商品信息； 3）新建商品信息模块：建立一个新结构体，为链表存信息用，并且将信息保存在指定的文件中； 4）排序模块：把文件中顺序零乱的商品信息按单价的大小从高到低进行排序，放到链表里存储； 5）计算模块：将所有商品的价格与库存量进行累加求和； 6）打印模块：将商品信息按每类平均价格排序（从高到低）按顺序打印出来； 7）统计模块1：统计库存量低于100的货名及类别；
8）统计模块2：统计商品库存量有2种以上（含2种）低于100的商品类别。 附录（程序清单） #include "stdio.h" /*输入，输出头文件*/ #include "stdlib.h" /*申请空间头文件*/ #include "string.h" /*对字符串加工头文件*/ #include "conio.h" /*清屏头文件*/ FILE *fp; int n=0; /*定义文件指针类型*/ int i,j,a[4],m; /*定义整数类型*/ float aver[4],sum[4],g[4],h; /*定义浮点类型*/ char c[5]="elec"; /*定义字符数组类型*/ char e[5]="food"; /*定义字符数组类型*/ char f[5]="offi"; /*定义字符数组类型*/ struct good /*定义结构体*/ { int num; /*商品编号*/ char name[20]; /*商品名称*/ char kind[40]; /*商品类型*/ float price; /*商品价格*/ char unit[10]; /*商品单位*/
#include
软件著作权-源代码范本 注意事项：常见的源代码包含：C语言，VB,C++,JAVA,.NET等。 提交的代码必须是源代码的开头载入程序，第30页必须断开，第60页是软 件的程序结尾，代码中不得出现与申请表内容不符合的日期，著作权人，软件名 字等，不能出现开源代码，不能出现任何版权纠纷。 格式要求：一、源代码应提交前、后各连续30页，不足60页的，应当全部提交。 、源代码页眉应标注软件的名称和版本号，应当与申请表中名称完全一致，页 眉右上应标注页码，源代码每页不少于50行。 范例如下: #i nclude
#in elude
%*************************电力系统直角坐标系下的牛顿拉夫逊法潮流计算********** clear clc load E:\data\IEEE014_Node.txt Node=IEEE014_Node; weishu=size(Node); nnum=weishu(1,1); %节点总数 load E:\data\IEEE014_Branch.txt branch=IEEE014_Branch; bwei=size(branch); bnum=bwei(1,1); %支路总数 Y=(zeros(nnum)); Sj=100; %********************************节点导纳矩阵******************************* for m=1:bnum; s=branch(m,1); %首节点 e=branch(m,2); %末节点 R=branch(m,3); %支路电阻 X=branch(m,4); %支路电抗 B=branch(m,5); %支路对地电纳 k=branch(m,6); if k==0 %无变压器支路情形 Y(s,e)=-1/(R+j*X); %互导纳 Y(e,s)=Y(s,e); end if k~=0 %有变压器支路情形 Y(s,e)=-(1/((R+j*X)*k)); Y(e,s)=Y(s,e); Y(s,s)=-(1-k)/((R+j*X)*k^2); Y(e,e)=-(k-1)/((R+j*X)*k); %对地导纳 end Y(s,s)=Y(s,s)-j*B/2; Y(e,e)=Y(e,e)-j*B/2; %自导纳的计算情形 end for t=1:nnum; Y(t,t)=-sum(Y(t,:))+Node(t,12)+j*Node(t,13); %求支路自导纳 end G=real(Y); %电导 B=imag(Y); %电纳 %******************节点分类************************************* * pq=0; pv=0; blancenode=0; pqnode=zeros(1,nnum); pvnode=zeros(1,nnum); for m=1:nnum; if Node(m,2)==3 blancenode=m; %平衡节点编号 else if Node(m,2)==0 pq=pq+1; pqnode(1,pq)=m; %PQ 节点编号 else if Node(m,2)==2 pv=pv+1; pvnode(1,pv)=m; %PV 节点编号 end end end end %*****************************设置电压初值********************************** Uoriginal=zeros(1,nnum); %对各节点电压矩阵初始化 for n=1:nnum Uoriginal(1,n)=Node(n,9); %对各点电压赋初值 if Node(n,9)==0;
源代码范例参考 (一套完套的设计至少包括可以运行源程序,论文,开题,答辩,外文,以下列出的源代码仅为参考,仅是源程序中的部分代码) ------------------------------------------------------ 5 新闻系统的主要开发过程 5.1 数据库设计 该系统使用Access作为应用程序的数据库，建立了有6张数据表：Lb_Admin表用于存储帐户的信息，Lb_Class表用于存储新闻分类信息，Lb_Links 表用于存储友情链接信息, Lb_News表用于存储新闻内容, Lb_Setup表用于存储系统设置信息,Lb_Skins表用于存储页面风格。Lb_News表与Lb_Class表相关联。 在6个数据表中，主要介绍一下新闻表和账户表如下： Lb_News表如表1所示： 表1 新闻表 n_id为该表主键，n_cid为所属新闻分类ID，n_hit为点击数，n_commend 用于判断是否为推荐新闻，n_index用于判断是否为首页新闻，n_title为新闻标题，n_date为发布新闻的日期，n_author为新闻发布者名，n_content为新闻内容，n_ispic用于判断是否为首页滚动图片，n_spic为首页所滚动图片的存储地址，,n_check用于判断是否为审核新闻(值为0表示未通过审核，值1表示通过
审核)，n_intro为新闻简介，n_publisher为账户类型，如为管理员、审核员或新闻输入员。 Lb_Admin表如表2所示： 表2 账户表 A_Id为该表主键，A_Name为用户名，A_Password为密码，A_Grade为用户级别。 5.2各主要窗体页的图形界面介绍 ◆首页：主要功能是显示最热门、最新和推荐的前8条新闻，显示新闻搜索功能、页面风格设置、新闻分类、滚动新闻图片和网站公告等； ◆新闻详细内容查看页：在首页里点击新闻标题所进入的页面，显示当前所点击新闻的详细内容，新闻标题正中加粗显示，在添加新闻时，根据所编辑的新闻内容对应显示；
第二章 任务一：闪烁广告灯的设计 利用89C51单片机的端口控制两个LED ( DO和D1 ),编写程序，实现两个LED互闪。 #include
旗开得胜CRM系统：开创源码购买先河 信息化时代的到来，让互联网开始迅速融入我们的生活工作中，我们不仅需要依靠网络来了解社会动态，而且还需要依靠网络来统筹企业的发展。尤其是一些成熟运作的办公软件，不仅能更好节约工作中的人力物力，甚至还能助力企业的稳定发展。 旗开得胜CRM系统即是在网络信息化发展浪潮中应运而生的一款客户管理系统，它维系着企业与客户之间的良好关系，亦承载着企业发展过程中的重要商业战略。它不仅有着同行业CRM 系统的优势，同时也领先行业的发展，开创了源码购买的先河。依托源码开发构建自主品牌 就目前互联网行业中的各种付费系统及软件而言，基本上还处于一个“拿来主义”阶段，虽然使用方式快捷方便，但是也要依托系统本身的套路来走，不能根据自身的需求来制定“规则”。而旗开得胜的源码开发恰恰弥补了客户管理系统中的“套路式”发展，更有利于成就企业发展中的“独一无二”。 CRM源码开发的重要性就相当于买房子产权的重要性，没有源码开发权的CRM系统就好比只有使用权却没有产权的房子，我们就算住着也不踏实。 目前，很多从事互联网软件行业的人士面临着发展的困境，自建团队成本太高，开发周期太长，而代理别人的品牌成本又太高，旗开得胜的CRM源码招商方案则有效解决了这两大难题。只要拥有源码你就可以创建属于自己的品牌，并且可以在系统标
准版的基础上，做出各行业的版本，满足自己现有客户资源的需求，快速帮助客户实现增值服务。 移动系统助力员工旗开得胜 除却利用源码开发对自主品牌的搭建外，旗开得胜CRM系统本身还有着无法比拟的优势。它结合移动端的创新功能应用，将移动、社交、云智慧融合，实现了PC端和移动端的完美结合，整合了碎片化时间，让员工可以随时随地掌握工作进展。 此外，为了保证数据的安全、系统的稳定，旗开得胜CRM 系统还采用了银行级的数据运算方案，不仅实现了客户资料的集中化记录，杜绝了因为人员变动而导致客户流失状况的出现，而且让销售管理变得更精细化，更自动化。 简约、便捷、易用的经营理念使旗开得胜CRM系统在短短十年内成为了数以万计企业选择的标准，而此次源码招商方案的启动，更是软件系统发展史上又一具有里程碑意义的尝试。 不管你是软件开发人员想要创业，还是一些软件代理商苦恼没有自己的品牌，甚至一些销售人员空有客户资源，只需要一点资金，您就可以选择旗开得胜CRM源码构建属于自己的系统平台，不仅省去了招聘开发团队的成本、时间与精力，同时还可以享受政府的补贴，实现真正意义上的高收益，高回报！
最常用的一些图像处理Matlab源代 码 #1：数字图像矩阵数据的显示及其傅立叶变换 #2：二维离散余弦变换的图像压缩 #3：采用灰度变换的方法增强图像的对比度 #4：直方图均匀化 #5：模拟图像受高斯白噪声和椒盐噪声的影响 #6：采用二维中值滤波函数medfilt2对受椒盐噪声干扰的图像滤波 #7：采用MATLAB中的函数filter2对受噪声干扰的图像进行均值滤波 #8：图像的自适应魏纳滤波 #9：运用5种不同的梯度增强法进行图像锐化 #10：图像的高通滤波和掩模处理 #11：利用巴特沃斯（Butterworth）低通滤波器对受噪声干扰的图像进行平滑处理 #12：利用巴特沃斯（Butterworth）高通滤波器对受噪声干扰的图像进行平滑处理 1.数字图像矩阵数据的显示及其傅立叶变换 f=zeros(30,30); f(5:24,13:17)=1; imshow(f, 'notruesize'); F=fft2(f,256,256); % 快速傅立叶变换算法只能处矩阵维数为2的幂次，f矩阵不 % 是，通过对f矩阵进行零填充来调整 F2=fftshift(F); % 一般在计算图形函数的傅立叶变换时，坐标原点在 % 函数图形的中心位置处，而计算机在对图像执行傅立叶变换 % 时是以图像的左上角为坐标原点。所以使用函数fftshift进 %行修正，使变换后的直流分量位于图形的中心； figure,imshow(log(abs(F2)),[-1 5],'notruesize');
2 二维离散余弦变换的图像压缩I=imread('cameraman.tif'); % MATLAB自带的图像imshow(I); clear;close all I=imread('cameraman.tif'); imshow(I); I=im2double(I); T=dctmtx(8); B=blkproc(I,[8 8], 'P1*x*P2',T,T'); Mask=[1 1 1 1 0 0 0 0 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; B2=blkproc(B,[8 8],'P1.*x',Mask); % 此处为点乘(.*) I2=blkproc(B2,[8 8], 'P1*x*P2',T',T); figure,imshow(I2); % 重建后的图像 3.采用灰度变换的方法增强图像的对比度I=imread('rice.tif'); imshow(I); figure,imhist(I); J=imadjust(I,[0.15 0.9], [0 1]); figure,imshow(J); figure,imhist(J);
超市商品管理系统源代码 import java.awt.*; import java.awt.event.*; import javax.swing.*; import java.sql.*; public class sale2 extends JFrame implements ActionListener //销售界面 { private JPanel sM=new JPanel(); //supermarket面板 JLabel label1,label2; JTextField field1,field2; JButton button1,button2; String num=null,name=null,time=null,add=null,t1=null,t2=null,t3=null; int price=0,countBG,n1;// 价格/商品的卖出数量 private JTextArea goodsShow=new JTextArea(); //文本区 private JScrollPane jsp=new JScrollPane(goodsShow); //滚动区域设定内容为商品显示的文本区 private JButton[] pS={new JButton("销售界面"), new JButton("销售管理"), new JButton("商品管理"), new JButton("用户管理")}; //personty 按钮 private JLabel[] CD={new JLabel("清单号",JLabel.CENTER), new JLabel("编号",JLabel.CENTER), new JLabel("名称",JLabel.CENTER), new JLabel("价格",JLabel.CENTER), new JLabel("生产日期",JLabel.CENTER), new JLabel("产地",JLabel.CENTER), new JLabel("数量",JLabel.CENTER)}; //帐单查询时的商品显示标签 public sale2() { label1=new JLabel("按编号查询",JLabel.CENTER); label2=new JLabel("按帐单查询",JLabel.CENTER); field1=new JTextField(30); field2=new JTextField(30); button1=new JButton("查询1"); button2=new JButton("查询2"); sM.setLayout(null);
正弦波的源程序： (一)，用到的函数 1,f2t函数 function x=f2t(X) global dt df t f T N %x=f2t(X) %x为时域的取样值矢量 %X为x的傅氏变换 %X与x长度相同并为2的整幂 %本函数需要一个全局变量dt(时域取样间隔) X=[X(N/2+1:N),X(1:N/2)]; x=ifft(X)/dt; end 2,t2f函数。 function X=t2f(x) global dt df N t f T %X=t2f(x) %x为时域的取样值矢量 %X为x的傅氏变换 %X与x长度相同,并为2的整幂。 %本函数需要一个全局变量dt(时域取样间隔) H=fft(x); X=[H(N/2+1:N),H(1:N/2)]*dt; end (二),主程序。 1，%（1）绘出正弦信号波形及频谱 global dt df t f N close all k=input('取样点数＝2^k, k取10左右'); if isempty(k), k=10; end f0=input('f0=取1(kz)左右'); if isempty(f0), f0=1; end N=2^k; dt=0.01; %ms df=1/(N*dt); %KHz T=N*dt; %截短时间
Bs=N*df/2; %系统带宽 f=[-Bs+df/2:df:Bs]; %频域横坐标 t=[-T/2+dt/2:dt:T/2]; %时域横坐标 s=sin(2*pi*f0*t); %输入的正弦信号 S=t2f(s); %S是s的傅氏变换 a=f2t(S); %a是S的傅氏反变换 a=real(a); as=abs(S); subplot(2,1,1) %输出的频谱 plot(f,as,'b'); grid axis([-2*f0,+2*f0,min(as),max(as)]) xlabel('f (KHz)') ylabel('|S(f)| (V/KHz)') %figure(2) subplot(2,1,2) plot(t,a,'black') %输出信号波形画图grid axis([-2/f0,+2/f0,-1.5,1.5]) xlabel('t(ms)') ylabel('a(t)(V)') gtext('频谱图') 最佳基带系统的源程序： (一)，用到的函数 f2t函数和t2f函数。代码>> (二),主程序 globaldt t f df N T close all clear Eb_N0 Pe k=input('取样点数＝2^k, k取13左右'); if isempty(k), k=13; end z=input('每个信号取样点数＝2^z, z close all clear echo on clc % NEWFF——生成一个新的前向神经网络 % TRAIN——对 BP 神经网络进行训练 % SIM——对 BP 神经网络进行仿真 % 定义训练样本 % P为输入矢量 P=[0.7317 0.6790 0.5710 0.5673 0.5948;0.6790 0.5710 0.5673 0.5948 0.6292; ... 0.5710 0.5673 0.5948 0.6292 0.6488;0.5673 0.5948 0.6292 0.6488 0.6130; ... 0.5948 0.6292 0.6488 0.6130 0.5654; 0.6292 0.6488 0.6130 0.5654 0.5567; ... 0.6488 0.6130 0.5654 0.5567 0.5673;0.6130 0.5654 0.5567 0.5673 0.5976; ... 0.5654 0.5567 0.5673 0.5976 0.6269;0.5567 0.5673 0.5976 0.6269 0.6274; ... 0.5673 0.5976 0.6269 0.6274 0.6301;0.5976 0.6269 0.6274 0.6301 0.5803; ... 0.6269 0.6274 0.6301 0.5803 0.6668;0.6274 0.6301 0.5803 0.6668 0.6896; ... 0.6301 0.5803 0.6668 0.6896 0.7497]; % T为目标矢量 T=[0.6292 0.6488 0.6130 0.5654 0.5567 0.5673 0.5976 ... 0.6269 0.6274 0.6301 0.5803 0.6668 0.6896 0.7497 0.8094]; % Ptest为测试输入矢量 Ptest=[0.5803 0.6668 0.6896 0.7497 0.8094;0.6668 0.6896 0.7497 0.8094 0.8722; ... 0.6896 0.7497 0.8094 0.8722 0.9096]; % Ttest为测试目标矢量 Ttest=[0.8722 0.9096 1.0000]; % 创建一个新的前向神经网络 net=newff(minmax(P'),[12,1],{'logsig','purelin'},'traingdm'); % 设置训练参数 net.trainParam.show = 50; net.trainParam.lr = 0.05; net.trainParam.mc = 0.9; net.trainParam.epochs = 5000; net.trainParam.goal = 0.001; % 调用TRAINGDM算法训练 BP 网络 [net,tr]=train(net,P',T); % 对BP网络进行仿真 A=sim(net,P'); figure; plot((1993:2007),T,'-*',(1993:2007),A,'-o'); title('网络的实际输出和仿真输出结果，*为真实值，o为预测值'); xlabel('年份'); ylabel('客运量'); % 对BP网络进行测试 A1=sim(net,Ptest');BP神经网络matlab源程序代码
相关主题
文本预览