摘要: 前段时间改写老师的一个计算纳米细颗粒物的C语言程序,不过运行的效率实在是太低了,工作站跑都很吃劲,晚上看了冈萨雷斯的数字图像处理意外发现了问题所在。在将C语言的程序改写成Matlab程序的时候我还是保持的写C语言的想法,上来全是for循环而忽视了Matlab的矩阵操作。 这里举的例子是冈萨雷斯的书里面的一个例子,比较的是f(x,y)=A sin(u0x+v0y),用for循环和向量化代码实现,具体如下:function [rt,f,g] = twodsin(A,u0,v0,M,N)tic %开始时间for r = 1:M u0x = u0*(r-1); for c... 阅读全文
posted @ 2013-07-01 22:34 k-frank 阅读(615) 评论(0) 推荐(0) 编辑
摘要: 这几天给学弟的讲一些关于服务器的东西,这里总结了一些远程连接的工具的使用方法。 一般远程连接服务器的工具分为两大类:文字接口服务器和图形接口服务器,前者的代表有SSH、Putty,而后者的代表有XDMCP、VNC。我使用的是Putty和VNC,服务器操作系统是RedHat。1、文字接口服务器---Putty 选用Putty的一个主要考虑的是它是开放源代码软件。随着Linux在服务器端应用的普及,Linux系统管理越来越依赖于远程。Putty是一个免费的、Windows 32平台下的telnet、rlogin和ssh客户端,但是功能丝毫不逊色于商业的telnet类工具。 另外,因为我在使用的这. 阅读全文
posted @ 2013-05-14 12:51 k-frank 阅读(491) 评论(0) 推荐(0) 编辑
摘要: 好久没发博文了,这两天正好在做通原的实验,正好学学怎么用Live Writer。 通常我们会把模拟信号抽样、量化,直到最后变换成二进制符号的集成过程称为脉冲编码调制(Pulse Code Modulation),实现的主要包括三个步骤:抽样(sampling)、量化(quantization)和编码(coding)。 1、Sampling Matlab仿真过程中选定模拟信号sin(4*t),按抽样定理,抽样频率只需高于模拟信号最高频率的两倍即可,但是实际情况下会大于两倍的模拟最高频率。 2、Quantization 量化方法包括:均匀量化和非均匀量化。这里主要讨论非均匀量化,并模拟A率压缩,. 阅读全文
posted @ 2013-05-10 16:14 k-frank 阅读(4134) 评论(0) 推荐(0) 编辑
摘要: MCM/ICM结束有段日子了,过年回家几天天天在外面走亲戚吃喝的。比赛那会儿天天盼着放假回家过年,可以休息啊,结果,回家过个年,走亲戚吃饭比比赛还累人、、、 按一个基友的观点就是我从数学建模这儿得到的荣誉太多了,不过我更多的是觉得,在这里我学到了很多,成长了很多,这是我想说的。现在回首我的数学建模之路走的其实并不是很好,主要还是位置没认清楚,没站好位置。 参加数学建模竞赛第一个成长的就是团队协作能力。我想工科生男一般都有这毛病,沉默寡言,一天到晚只顾做自己的事儿。其实现在的我还是的,相比以前稍微好点。我所在的这只队伍从国赛,一直走到了美赛,肯定不是最优秀的,但肯定也蛮奇葩的。另外两个人... 阅读全文
posted @ 2013-02-22 21:57 k-frank 阅读(446) 评论(1) 推荐(0) 编辑
摘要: 下面四段代码主要实现了单极性归零、单极性不归零、双极性归零、双极性不归零四个信号的变化function y = dnrz(x)%本函数实现将输入的一段二进制代码编为相应的双极性非归零码输出%输入x为二进制码,输出y为编好的码%给出计算每一个码元的点数,因为我们只是用离散的点来得出连续的函数表达grid = 300;t = 0:1/grid:length(x);for i = 1:length(x) %计算码元的值 if(x(i) == 1) %如果信息为1 for j = 1:grid %该码元对应的点值为1 ... 阅读全文
posted @ 2012-11-08 14:55 k-frank 阅读(876) 评论(0) 推荐(0) 编辑
摘要: 我们知道Fourier分析是信号处理里很重要的技术,matlab提供了强大的信号处理能力,但是有一些细节部分需要我们注意。记信号f(t)的起始时间为t_start, 终止时间为t_end, 采样周期为t_s, 可以计算信号的持续时间Duration为 t_end – t_start, 信号离散化造成的采样点数 N = Duration/t_s + 1;根据Fourier分析的相关结论,我们知道时域的采样将会造成频域的周期化,该周期为采样频率f_s(著名的香农采样定理基于此).于是, 经过matlab的fft函数处理后,得到数据的横坐标为0:f_s/(N-1):f_s。相关代码如下所示:%mat 阅读全文
posted @ 2012-11-08 14:51 k-frank 阅读(3056) 评论(0) 推荐(0) 编辑
摘要: matlab神经网络工具箱主要用于BP网络分析与设计的函数,这面这些函数最好都弄懂:1、前向网络创建函数:newcf、newff和newfftd2、激励函数:logsig、dlogsig、(S型对数式)tansig、dtansig、(正切函数式)purelin、dpurelin(线性函数)3、学习函数:learngd(梯度下降权值/阀值学习函数)、learndm(梯度下降动量函数)4、训练函数:trainbfg(BFGS准牛顿BP算法)、traingd(梯度下降BP算法训练函数)、traingdm(梯度下降动量BP算法)等等5、性能函数:mse(均方误差性能函数)、msereg函数学习这一块在 阅读全文
posted @ 2012-11-02 15:09 k-frank 阅读(1955) 评论(0) 推荐(3) 编辑
摘要: clear all,clcxn = [1 0 1 1 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0]; %输入单极性码yn = xn; %输出yn初始化num = 0; %计数器初始化for k = 1:length(xn) if xn(k)==1 num = num+1; %'1'计数器 if num/2 == fix(num/2) %奇数输出-1,偶数输出1,进行极性交替 yn(k) = 1; else yn(k) ... 阅读全文
posted @ 2012-10-31 21:39 k-frank 阅读(1824) 评论(0) 推荐(0) 编辑
摘要: huffman.mfunction [h,l]=huffman(p)if (length(find(p<0))~=0) error('Not a prob,negative component'); end if (abs(sum(p)-1)>10e-10) error('Not a prob.vector,component do not add to 1') end n=length(p); q=p; m=zeros(n-1,n); for i=1:n-1 [q,l]=sort(q); m(i,:)=[l(1:n-i+1),zeros(1,i-1 阅读全文
posted @ 2012-10-27 17:26 k-frank 阅读(425) 评论(0) 推荐(0) 编辑
摘要: 微软那边下的directx11,刚开始学习语音视频开发,顺便做点笔记吧,。下载安装完后,打开vs2008--->工具--->选项--->---项目与解决方案--->C++目录,显示以下内容里选择“包含文件”,然后“新行”或者“ctrl-insert”,插入C:\Program Files\Microsoft DirectX SDK (June 2010)\Include和C:\Program Files\Microsoft DirectX SDK (June 2010)\Lib\x86,另外需要注意的是C:\Program Files\Microsoft DirectX 阅读全文
posted @ 2011-11-02 19:00 k-frank 阅读(495) 评论(0) 推荐(0) 编辑