Todd's Matlab讲义第6讲：割线法

割线法割线法求解方程$f(x)=0$的根需要两个接近真实根$x^\*$的初值$x_0$和$x_1$,于是得到函数$f(x)$上两个点$(x_0,y_0=f(x_0))$和$(x_1,y_1=f(x_1))$,连接这两点得到一条直线(割线): \begin{equation*}y-y_1=\frac{y_1-y_0}{x_1-x_0}(x-x_1)\end{equation*} 由于我们要求解$f(x)=0$,因此设$y=0$,由上式解出$x$,作为下次迭…

Todd's Matlab讲义第5讲：二分法和找根

二分法和if ... else ... end 语句先回顾一下二分法.要求方程$f(x)=0$的根.假设$c = f(a) < 0$和$d = f(b) > 0$,如果$f(x)$是连续函数,那么方程的根$x^*$一定位于$a$和$b$之间.然后,我们看一下$a$和$b$中点$x=(a+b)/2$,计算函数值$y=f(x)$,如果函数不为0,比较$c$.$d$和$y$的符号,确定新的二分区间.具体来说,如果$c$和$y$同号,新…

Todd's Matlab讲义第4讲：控制误差和条件语句

误差和残量数值求解方程$f(x)=0$的根,有多种方法测算结果的近似程度.最直接的方法是计算误差.第$n$步迭代结果与真值$x^\*$的差即为第$n$步迭代的误差: \begin{equation*}e_n=x_n-x^*\end{equation*} 但是,我们一般是不知道真实值$x^\*$的,否则,我们也不会费劲去算了.所以,直接计算误差是不可能的,需要我们另辟蹊径. 一个可能的方法是,程序一直运行,直到结果不再变化.这个方法通常还是很管用的.有时候,程序结果不再变化并…

Todd's Matlab讲义第3讲：牛顿法和for循环

方程数值求解下面几讲,我们将聚集如下方程的解法: \begin{equation} f(x)=0 \tag{3.1}\label{3.1} \end{equation} 在微积分课程中,我们知道,许多优化问题最终归结为求解上述形式的方程,其中$f$为你要求极值的函数$F$的导数.在工程问题中,函数$F$来源多种多样,有公式.微分方程的解.实验和模拟等. 牛顿迭代我们把方程\eqref{3.1}的解记为$x^\*$.方程的解法有三种:对分法.割线法和牛顿法.这三种方法都需要猜测…

Todd's Matlab讲义第2讲：Matlab 编程

Matlab也可以编程,可存为以.m为后缀的文件,称为M文件.M文件有两种:函数和脚本. 函数程序点击新建图标,在打开的窗口里输入如下内容: function y = myfunc (x) y = 2*x.^2 - 3*x + 1; 将文件保存为myfunc.m,保存在当前目录下.这个文件就可以直接在命令窗口使用了,用法如Matlab内置函数,如在命令窗口输入如下内容: >> x = -2:.1:2; >> y = myfunc(x); >> plot(x,y) 这里…

Todd's Matlab讲义第1讲：向量，函数和作图

向量 Matlab 中最基本的对象是矩阵,向量是特殊的矩阵.行向量是$1\times n$矩阵,列向量是$m\times 1$矩阵.输入如下行向量: >> v=[0 1 2 3] 回车,Matlab 将输出此行向量: v = 0 1 2 3 输入如下列向量: >> u = [9; 10; 11; 12; 13] 按如下方式访问向量的元素: >> u(2) 按如下方式改变元素的值: >> u(2)=47 按如下方式提取多个连续向量: >>…

基于MATLAB的GUI（Graphical User Interface）音频实时显示设计

摘要:本文章的设计主要讲基于matlab的gui音频实时显示设计,此次设计的gui相当于一个简洁的音乐播放器,界面只有”录音“和”播放“两个控件,哈哈,够简洁吧.通过”录音“按钮可以实现声音从电脑的声卡录入,并且实时显示录入声音的时域图形和频域图形:待录音结束,通过”播放“按钮可以播放刚录入的声音,并且一边播放一遍实时显示时域和频域图形.本设计的编码在matlab2013a上亲测,可以实现... 一.首先matlab的gui界面设计打开matlab—>在命令行执行guide—>出现gui编辑…

Bag of Words/Bag of Features的Matlab源码发布

2010年11月19日 ⁄ 技术, 科研 ⁄ 共 1296字 ⁄ 评论数 26 ⁄ 被围观 4,150 阅读+ 由于自己以前发过一篇文章讲bow特征的matlab代码的优化的<Bag-Of-Words中K-Means聚类的效率优化>,其中的代码也用过Spatial Pyramid Code的代码里面的几个函数,不过大部分还是根据本地的需要,根据大数据量计算的需要自己整合修改的,经过不少同学的改错,现在已经基本没有错误了,注释没怎么写,以后慢慢补上,如果有什么问题可以交流. 下载地址是我的实验室…

Coursera台大机器学习技法课程笔记03-Kernel Support Vector Machine

这一节讲的是核化的SVM,Andrew Ng的那篇讲义也讲过,讲的也不错. 首先讲的是kernel trick,为了简化将低维特征映射高维特征后的计算,使用了核技巧.讲义中还讲了核函数的判定,即什么样的函数K能使用kernel trick. 此外,核函数还可以衡量两个特征的相似度,值越大,越相似(可用于跟踪哦 ). 接着讲的是Polynomial Kernel,需要注意的是核函数的系数和常量,这会影响到最终的margin. 再接着讲了高斯Kernel,它能将原始数据映射到无限维!但是如果参数选的…

转自 z55250825 的几篇关于FFT的博文（一）

关于FFT,咱们都会迫不及待地 @ .....(大雾)(貌似被玩坏了...) .....0.0学习FFT前先orz FFT君. 首先先是更详细的链接(手写版题解点赞0v0) FFT的资料其实众所周知的最详细的算法解释在<算法导论>上...然后咱就是边看着那个边码理解的... 首先来看看多项式乘法和快速FFT的关系,然后咱们再来看能否聊到卷积什么的东西... 其实觉得还是去看算法导论最好. [一.多项式及其表达方式.] 首先什么是…

Deep Learning 学习随记（五）Deep network 深度网络

这一个多周忙别的事去了,忙完了,接着看讲义~ 这章讲的是深度网络(Deep Network).前面讲了自学习网络,通过稀疏自编码和一个logistic回归或者softmax回归连接,显然是3层的.而这章则要讲深度(多层)网络的优势. Deep Network: 为什么要使用深度网络呢?使用深度网络最主要的优势在于,它能以简洁的方式来表达比浅层网络大得多的函数集合.正式点说,可以找到一些函数,它们能够用k层网络简洁的表达出来(这里的简洁指的是使用隐层单元的数目与输入单元数目是多项式关系),但是对一…

WC2015流水账

THU那四场考试没考好,只有20+名.这也许是我OI生涯中最后一场吧(已确认是最后一场),真是感慨万千. day0 搬进浙大宿舍404房间(神房间号),四个人一间.中午发现学军伙食相当良心,是我参加的OI活动中最良心的 day1 徐寅展讲了一些动态维护图的连通性以及最小生成树的算法,太深了,基本上只听懂第一个维护图的连通性的算法.下面就简单讲一下这个算法吧. 首先,这个算法是在线算法,其时间复杂度是m lg^2.离线算法可以做到m lg,用LCT维护最大生成树就行了(边权为边被删除的时间).…

BaseAdapter自定义适配器——思路详解

BaseAdapter自定义适配器——思路详解引言: Adapter用来把数据绑定到扩展了AdapterView类的视图组.系统自带了几个原生的Adapter. 由于原生的Adapter视图功能太少,有时需要有自己的视图格式.并且在开发中经常用到. 实现效果: 实现过程: 重写视图区域: 即重写BaseAdapter 有四方要重写: public int getCount() 数据域大小方法 public Object getItem(int arg0) 返回每个Item的数据据 public…

GDKOI2018发烧记

偏远小渔村NOIP螺旋升天选手又一次来到了广州参加GDKOI...金实的初三爷们也来啦?要被碾啦T T Day 0 跟HR Lao爷拼(biao)车到了高铁站,上了高铁居然没有颓颓颓吃吃吃(雾),安心地打了打板子(姚老坐在旁边T T 到广州,愣是盯了半天地铁网,从广州南站坐到苏元站,woc怎么离市中心这么远啊...人生中第一次坐这么久地铁... 出了地铁站,面前就是一个van♂da广场,说好的郊区呢T T 怎么比偏远小渔村还繁华的样子啊喂... 到了酒店,房间距离其他大爷的房间隔了好几栋楼几百米…

hiho一下第115周：网络流一•Ford-Fulkerson算法（Edmond-Karp，Dinic，SAP）

来看一道最大流模板水题,借这道题来学习一下最大流的几个算法. 分别用Edmond-Karp,Dinic ,SAP来实现最大流算法. 从运行结过来看明显SAP+当前弧优化+gap优化速度最快. hiho一下第115周:网络流一•Ford-Fulkerson算法原题网址:http://hihocoder.com/contest/hiho115/problem/1 网络流一·Ford-Fulkerson算法时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述小Hi和…

【CSAPP】Attack Lab实验笔记

attacklab这节玩的是利用一个字符串进行缓冲区溢出漏洞攻击,就小时候想象中黑客干的事儿. 做题的时候好几次感叹这些人的脑洞,"这都可以攻击?还能这么注入?这还可能借力打力?"等自己注入的时候却是"啊?怎么又段错误了?怎么又算错地址了?"也是一次有趣的经历了. 小插曲:我拿到文件的时候直接去读得readme,看完了还迷惑这readme咋就这么点信息.后来知道了实验都要配合着writeup讲义看,不禁感叹我前两个实验没看讲义还能做出来真是个奇迹! level1 第…

matlab实现不动点迭代、牛顿法、割线法

不动点迭代 function xc = fpi( g, x0, tol ) x(1) = x0; i = 1; while 1 x(i + 1) = g(x(i)); if(abs(x(i+1) - x(i)) < tol) break end i = i + 1; end xc = x(i+1); end 牛顿法: function xk = funNewton(f, x0, max_steps, tol) syms x symbol_f = f(x); dif_f = matlabFunct…

Matlab高级教程_第二篇：MATLAB和C#对应数据类型的讲解（多讲一点儿C#的矩阵运算）

1. MATLAB对应C#的数据类型主要在引入的父类库MWArray当中.有如下对应规则 .NET TYPE MWArrayTYPE MATLAB Type System.Double MWNumericArray double System.Number MWNumericArray double System.Float MWNumericArray single System.Byte MwNumbericArray int8 System.Short MWNumbericArray …

来来来，有讲一个吐血的故事（matlab）之脚本运行路径是什么

脚本运行路径是什么,这真是太重要!! 重要1:你默认保存的路径重要2:你访问的相对路径先放图: 再看一幅图: 我的操作,点击左侧的文件夹,使上框的显示栏路径不一样,再点击运行,发现pwd指示的路径不一样.当然身亡.脚本的默认路径令人堪忧.难怪我总是错路径.addpath 总是出错. 解决方法: (1)按照上述的,点左边的文件夹,使上框显示的路径和你心仪的路径(你要运行的脚本路径)一致. (2)法2,run(‘你脚本的路径’) (3)附上永久修改默认路径的方法(不推荐) http://jing…

MATLAB 编程风格指南及注意事项

MATLAB编程风格指南Richard Johnson 著Genial 译MATLAB 编程风格指南Richard JohnsonVersion 1.5,Oct. 2002版权: Datatool 所有翻译:Genial @ USTC“Language is like a cracked kettle on which we beat tuned to dance to, while allthe time we long to move the stars to pity.”- Gustave…

非线性方程(组)：MATLAB内置函数 solve, vpasolve, fsolve, fzero, roots [MATLAB]

MATLAB函数 solve, vpasolve, fsolve, fzero, roots 功能和信息概览求解函数多项式型非多项式型一维高维符号数值算法 solve 支持,得到全部符号解若可符号解则得到根支持支持支持当无符号解时符号解方法:利用等式性质得到标准可解函数的方法基本即模拟人工运算 vpasolve 支持,得到全部数值解 (随机初值)得到一个实根支持支持 $\times$ 支持未知 fsolve 由初值得到一个实根由初值得到一个实根支持支持…

PYTHON替代MATLAB在线性代数学习中的应用(使用Python辅助MIT 18.06 Linear Algebra学习)

前言 MATLAB一向是理工科学生的必备神器,但随着中美贸易冲突的一再升级,禁售与禁用的阴云也持续笼罩在高等学院的头顶.也许我们都应当考虑更多的途径,来辅助我们的学习和研究工作. 虽然PYTHON和众多模块也属于美国技术的范围,但开源软件的自由度毕竟不是商业软件可比拟的. 本文是一篇入门性文章,以麻省理工学院(MIT) 18.06版本线性代数课程为例,按照学习顺序介绍PYTHON在代数运算中的基本应用. 介绍PYTHON代数计算的文章非常多,但通常都是按照模块作为划分顺序,在实际应用中仍然有较多…

Matlab 高斯_拉普拉斯滤波器处理医学图像

前言:本程序是我去年实现论文算法时所做.主要功能为标记切割肝脏区域.时间有点久,很多细节已经模糊加上代码做了很多注释,因此在博客中不再详述. NOTE: 程序分几大段功能模块,仔细阅读,对解决医学图像还是有一定的借鉴意义想借鉴本文的一定要仔细阅读代码和注释,中间有人机交互部分,空跑会抛异常 .dcm数据,我放到了我的百度云盘,有兴趣的可以下载,实测一下代码.dcm数据连接 clc,clear img_1=dicomread('10011.dcm');%读取dcm文件 (所谓的灰度值) meta…

我的公司培训讲义（1）：.NET开发规范教程

这是1年多以前我在公司所做讲座的讲义,现在与园友们分享,欢迎拿去使用.一起讨论.文中有若干思考题,对园友们是小菜一碟.另有设计模式讲义一篇,随后发布.博文上了首页,感谢博客园团队推荐,也感谢所有园友的支持.经园友提醒,我们推荐新手使用ReSharper插件,它可以自动提示不符合代码规范的地方. 1 概述 1.1 意义 “.NET平台开发规范”包括两方面内容:代码规范.最佳实践. (1) 规范和优美的代码风格有助于思路清晰. (2) 规范和优美的代码有助于团队沟通交流. (3) 规范和优美的代码有…