1. 分别用三点和四点Gauss-Chebyshev公式计算积分 并与准确积分值2arctan4比较误差.若用同样的三点和四点Gauss-Legendre公式计算,也给出误差比较结果. 2*atan(4) ans = 2.6516 Gauss-Chebyshev: function I = gausscheby(f,a,b,n) syms t; t= findsym(sym(f)); ta = (b-a)/2; tb = (a+b)/2; switch n case 3 I=pi/n*ta*(s

由两点坐标如何画出直线  方法1:利用直线方程 斜率加截距 方法2:数据拟合 %由两点坐标得数据拟合直线与画线 x = [,]; y = [,]; k = ((-)/(-));% 由两点坐标得到直线斜率 line = k*x+0.5;% 直线方程 xy = :;% 定义画线的 x 长度 line1 = k*xy+0.5; figure(),plot(xy,line1); % % x1=[,]; x2=[,]; X=[x1(),x2()];% 两点坐标的x值 Y=[x1(),x2()];% 两点坐

前言: 本贴写于2016年12与15日,UK.最近在学习EMD(Empirical Mode Decomposition)和HHT(Hilbert-Huang Transform)多分辨信号处理,FQ回百度(好奇怪,哪儿的墙)寻找资源,然而却发现大部分的资料对于基于Matlab的EMD的介绍和安装介绍,不够详细,而且在安装和使用过程中遇到的问题也未能总结一二,然而,还有一些C某DN的作者,Goolge免费的资源你要大家5分去下载?起码也得1分,本贴介绍的资源下载,放在C某DN,每个资源1分(为满

MATLAB地图工具箱学习总结(三)地图工具箱的基本知识 今天想要介绍的是一些比较基础的函数.了解了这些函数,地图投影的基本概念才能真正明白.而要想继续研究MATLAB中有关地图投影的函数,尤其是未来我要提到的投影文件源代码,知晓这些函数的功能必不可少.本篇文章将会罗列三个案例,并在后面一一进行讲解. 1                    作业案例:地图投影作业1 这次的案例从作业1开始.作业1是要求计算出地球椭球体的一些基本参数,包括子午圈曲率半径.卯酉圈曲率半径.平均曲率半径和纬圈半径等

MATLAB地图工具箱学习总结(二)大圆和恒向线 今天要和大家谈一谈大圆.恒向线航道的画法.还是先从案例开始说起,再分别介绍相关的函数. 1                    作业案例:地图投影作业6 在作业6中,李老师要求我们画出两点间的大圆航线.恒定线航线和两者混合的航线.我一开始没有找到gcwaypts函数,导致自己花了一点时间编写了一段代码实现,但效果并不好,存在误差.后来找到了神奇的gcwaypts函数,一切也就迎刃而解.在本文中,我将会把两种方法都简单介绍一下. 首先是gcway

FFT(Fast Fourier Transformation),即为快速傅氏变换,是离散傅氏变换(DFT)的快速算法. 采样得到的数字信号,做FFT变换,N个采样点,经过FFT之后,就可以得到N个点的FFT结果.为了方便进行FFT运算,通常N取2的整数次方. 假设信号: S=2+3*cos(2*pi*50*t-pi*30/180)+1.5*cos(2*pi*75*t+pi*90/180) 它含有:2V的直流分量 频率为50Hz.相位为-30度.幅度为3V的交流信号 频率为75Hz.相位为90度

2015-12-17晚,复习,甚是无聊,阅<复杂网络算法与应用>一书,得知最小生成树问题(Minimum spanning tree)问题.记之. 何为树:连通且不含圈的图称为树. 图T=(V,E),|V|=n,|E|=m,下列关于树的说法等价: T是一个树. T无圈,且m=n-1. T连通,且m=n-1. T无圈,但每加一新边记得到唯一一个圈. T连通,但任舍去一边就不连通. T中任意两点,有唯一道路相连. 何为生成树:若图G=(V,E)的生成子图是一棵树,则称该树为图G的生成树,也称支撑树

还是计算几何, 多边形的核可以这样理解:这个核为原多边形内部的一个多边形,站在这个叫核的多边形中,我们能看到原多边形的任何一个位置. 算法步骤如下: 1.根据原多边形最大和最小的x,y初始化核多边形,就是个矩形. 2.计算多边形当前处理的点的凹凸性. 3.用当前点与其后继点构成直线,判断当前点的前驱点在该直线的左边或右边. 4.用该直线将原核多边形分为两个部分,选择其中一个部分作为处理下一个点将用到的核,选择的依据有以下两点: 1)如果当前点为凸点,那么选择的核与3步中前驱点的所在方向相同. 2

又是计算几何,我感觉最近对计算几何上瘾了. 当然,工作上也会用一些,不过工作上一般直接调用boost的geometry库. 上次写过最小包围圆,这次是最小包围矩形,要比最小包围圆复杂些. 最小包围矩形可不一定是个直立的矩形,也可能像下图一样是倾斜的. 求法如下: 1.求多边形凸包,这里凸包直接调用系统函数了,细节可以参考这里,虽然当时写的不怎么样. 2.将凸包两个相邻的点连线作为矩形一条边. 3.寻找凸包上距离已得到的边最远的点,过该点做平行线,得到矩形第二条边. 4.将凸包上点向已求得的边投影

(一)线性方程组求解 包含n个未知数,由n个方程构成的线性方程组为: 其矩阵表示形式为: 其中 一.直接求解法 1.左除法 x=A\b; 如果A是奇异的,或者接近奇异的.MATLAB会发出警告信息的. 2.利用矩阵的分解来求解线性方程组(比单单进行左除速度快) (1)LU分解(只有方阵可以使用) LU分解就是分解成一个交换下三角矩阵(也就是说进行一定的操作后才是下三角矩阵)和一个上三角矩阵(不需要变换)的乘积形式.只要A是非奇异的,就可以进行LU分解. MATLAB提供的LU分解函数对于矩阵进行

disp函数直接将内容输出在Matlab命令窗口中, 关键是看disp函数怎么把字符和数字在一起进行显示. matlab中disp()就是屏幕输出函数,类似于c语言中的printf()函数 %%以下是一个通过给定两点显示直线方程的程序, %%该程序需要给出两个点的坐标,结果返回为y=kx+b的格式,且求得斜率 function [k,a1,b,type]=straight_line(A,B) % 输入,A,B两点坐标 V=B-A; a=inf; b=inf; type='undefined';

总结有两点: 1.如需添加jar包等附加库,可在待安装工具箱下,新建一个java文件夹,并将jar包等文件存放在里面,然后执行以下操作: 1)Create or open your preferences file:    edit(fullfile(prefdir,'javaclasspath.txt')) % Matlab 2012b and newer    edit(fullfile(prefdir,'classpath.txt'))    % Matlab 2012a and olde

虽然Modelsim的功能非常强大,仿真的波形可以以多种形式进行显示,但是当涉及到数字信号处理的算法的仿真验证的时候,则显得有点不足.而进行数字信号处理是Matlab的强项,不但有大量的关于数字信号处理的函数,而且图形显示功能也很强大,所以在做数字信号处理算法的FPGA验证的时候借助Matlab会大大加快算法验证的速度. 关于Matlab和Modelsim联合仿真,我从网上看到两种方法,一种是通过Link for Modelsim建立Matlab和Modelsim的联合仿真接口:另一种就是通过文

今天实现了<一类求解方程全部根的改进差分进化算法>(by 宁桂英,周永权),虽然最后的实现结果并没有文中分析的那么好,但是本文依然是给了一个求解多项式全部实根的基本思路.思路是对的,利用了代数原理. 求解全部根的理论还是很有必要说一下的.就是利用了多项式综合除法,在matlab中可以采用deconv(A,B)直接实现.同时为了确定多项式方程根的范围,还采用了代数方程根的分布理论,个人觉得这两点是值得借鉴的一种方法. % 首先定义常量,包括最大迭代次数.搜索范围.个体维度.缩放因子等.程序如下

多普勒失真信号采样Matlab模拟分析 方案 水声通信指的是使用声信号在水中数据传输. 相对而言.电磁信号在水中吸收严重衰减过快,光信号受水中悬浮颗粒的影响,也无法完毕远距离传输. 这两种信号的传播距离约为数百米,而声信号在水中传播距离达到数十千米.成为水下通信的首选方式. 水声通信的一个重大挑战就是运动带来的多普勒. 多普勒造成信号的频率漂移,影响系统性能,因此很有必要消除多普勒带来的影响.而多普勒在时域的表现为信号的长度压缩和扩展.即'y(t) = x((1+a) * t)',当中a为多普勒

这一篇博客主要学习了: 第一,自定义函数优点有很多,比如可重复利用:容易纠错,以后直接装入大工程里(更重要的是,实参是复制给形参的,在自定义函数中运行时形参独立显示,这一点和C不一样). 第二,通过help/lookfor 自定义函数名 快速查找到自己所写的函数(H1注释行功能) 第三,调试功能 自定义函数代码: function distance = dist(x1, y1, x2, y2) %写在首行 %自定义函数 %脚本文件:dist.m %函数功能: % 该函数用来计算笛卡尔坐标系中两点

这几天学习了matlab的逻辑数组功能,总的感觉就有两点: 第一,通过对原来的数组a进行逻辑判断后得到逻辑数组b 第二,逻辑数组进行某种运算符操作又回到原数组类型 第三,利用逻辑数组mask功能 测试如下: disp('逻辑数组学习及应用'); a = [ ; ; ]; b = a > ;%数组中各个元素与4作比较,逻辑为真则为1 whos 结果: 第二,当执行 c = b + 0后,c 数组为: 第三,它在算术运算中能提供一个屏蔽(mask).屏蔽(mask)是指一个数组, 它从另一个数组选择

Matlab图论工具箱的命令见表1 表1  matlab图论工具箱的相关命令 命令名 功能 graphallshortestpaths 求图中所有顶点对之间的最短距离 graphconncomp 找无向图的连通分支,或有向图的强弱连通分支 graphisdag 测试有向图是否含有圈,不含圈返回1,否则返回0 graphisomorphism 确定两个图是否同构,同构返回1,否则返回0 graphisspantree 确定一个图是否是生成树,是返回1,否则返回0 graphmaxflow 计算有向

转载于博客:各种距离 在做分类时常常需要估算不同样本之间的相似性度量(SimilarityMeasurement),这时通常采用的方法就是计算样本间的"距离"(Distance).采用什么样的方法计算距离是很讲究,甚至关系到分类的正确与否. 本文的目的就是对常用的相似性度量作一个总结. 本文目录: 1.欧氏距离 2.曼哈顿距离 3. 切比雪夫距离 4. 闵可夫斯基距离 5.标准化欧氏距离 6.马氏距离 7.夹角余弦 8.汉明距离 9.杰卡德距离& 杰卡德相似系数 10.相关系数

一.向量.矩阵的表示和使用 format long  %小数很多format short %默认4位小数format rat %显示最近的分数format short e %指数格式的数 尾数多少 exp(1) %自然对数的底exp(10) log(x) %x的自然对数 底是e log10(x) %以10为底 asin(pi) atan(pi/3) %反三角函数 向量(vector)一维数值数组.MATLAB 允许你创建列向量和行向量,列向量通过在方 括号内把数值用分号(;)隔开来创建,对元素的

<写在前面> 关于台风轨迹图,这次由于作业需要就画了一幅轨迹图,参考原图选自论文: LIU Zenghong, XU Jianping, SUN Chaohui, WU Xiaofen. An upper ocean response to Typhoon Bolaven analyzed with Argo profiling floats [J]. Acta Oceanol. Sin., 2014, 33(11), 90–101 大家有兴趣可以自行查找论文进行阅读,本文所绘制台风轨迹图选取