% matlab script to demonstrate use of Gauss quadrature clear all close all % first derive the 2-point Gauss quadrature rule eq1 = 'w1*1 + w2*1 = 2'; eq2 = 'x1*w1 + x2*w2 = 0'; eq3 = 'x1^2*w1 + x2^2*w2 = 2/3'; eq4 = 'x1^3*w1 + x2^3*w2 = 0'; [w1,w2,x1,…

% script to perform adaptive quadrature clear all, close all global pts % function to be integrated defined in routine f f = 'integrand'; a = 1; b = 3; pts = [a;b]; tol = input('Enter error tolerance: '); % this is just to plot the graph % it's usual…

以下例子均来自网络,只是稍作了编辑,方便今后查阅. 子目录 (一) Fortran调用C语言 (二) C语言调用Fortran (三) C++ 调用Fortran (四) Fortran 调用 C++ 需要说明的是,(一)和(二)对GCC编译器的版本要求并不高:而(三)和(四)对GCC编译器的要求比较高,需要GCC在4.7及以上才能编译通过,这是由于自Fortran 2003一代语言起,增加了一个名为“iso_c_binding”的模块,可以很方便地在Fortran和C++之间传递参数,简化了两…

Triangular DGM 1. Basis functions decomposing the domain \(\Omega\) into \(N_e\) conforming non-overlapping triangular elements \(\Omega_e\). \[\begin{equation} \Omega = \bigcup_{e = 1}^{N_e} \Omega_e \end{equation}\] nonsingular mapping \(x = \Psi(\…

MATLAB数值积分法 作者:凯鲁嘎吉 - 博客园http://www.cnblogs.com/kailugaji/ 一.实验目的 许多工程技术和数学研究中要用到定积分,如果无法直接算不出精确值(如含在积分方程中的积分)或计算困难但可用近似值近似时,就用数值积分法方法加以解决.常用的算法有:复化梯形.辛甫生(Simpson).柯特斯(Cotes)求积法; 龙贝格(Romberg)算法:高斯(Gauss)算法. 二.实验原理 三.实验程序 下面给出复化Simpson求积法程序(梯形及柯特斯复化求积…

OpenCASCADE Gauss Integration eryar@163.com Abstract. Numerical integration is the approximate computation of an integral using numerical techniques. The numerical computation of an integral is sometimes called quadrature. The most straightforward nu…

支持向量机是建立在统计学习理论基础之上的新一代机器学习算法,支持向量机的优势主要体现在解决线性不可分问题,它通过引入核函数,巧妙地解决了在高维空间中的内积运算,从而很好地解决了非线性分类问题. 构造出一个具有良好性能的SVM,核函数的选择是关键．核函数的选择包括两部分工作:一是核函数类型的选择,二是确定核函数类型后相关参数的选择．因此如何根据具体的数据选择恰当的核函数是SVM应用领域遇到的一个重大难题,也成为科研工作者所关注的焦点,即便如此,却依然没有得到具体的理论或方法来指导核函数的选取 1.…

本文系<数字图像处理原理与实践(MATLAB版)>一书之代码系列的Part6,辑录该书第281至第374页之代码,供有须要读者下载研究使用.代码运行结果请參见原书配图,建议下载代码前阅读下文: 关于<数字图像处理原理与实践(MATLAB版)>一书代码公布的说明 http://blog.csdn.net/baimafujinji/article/details/40987807 P338 i=double(imread('vase.tif'));[C,S]=wavedec2(i,2,…

基于MATLAB的中值滤波均值滤波以及高斯滤波的实现 作者:lee神 1.   背景知识 中值滤波法是一种非线性平滑技术,它将每一像素点的灰度值设置为该点某邻域窗口内的所有像素点灰度值的中值. 中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术,中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替,让周围的像素值接近的真实值,从而消除孤立的噪声点. 方法是用某种结构的二维滑动模板,将板内像素按照像素值的大小进行排序,生成单调上升(或下降)的为二维数…

在matlab中经常需要对数据进行曲线拟合,如最常见的多项式拟合,一般可以通过cftool调用曲线拟合工具(curve fit tool),通过图形界面可以很方便的进行曲线拟合,但是有些时候也会遇到不方便用图形工具.因此这里简单的记下两种常用的拟合方法. 1 多项式拟合(polyfit和polyval) polyfit可以对数据进行拟合(自定义用几次多项式),返回相应的参数,然后用polyval生成拟合后的数据点,下面的例子中我们对抛物线y=3x2+6x+5进行拟合. x = -:; y = *…

% matlab script to derive the 2-point Gauss-Laguerre quadrature rule % and use it on an example % inelegantly set up equations from method of undetermined coefficients % and solve clear all close all format long eq1 = 'w1*1 + w2*1 = 1'; eq2 = 'w1*x1…

1.随机高斯测量矩阵 function [ Phi ] = GaussMtx( M,N ) %GaussMtx Summary of this function goes here % Generate Bernoulli matrix % M -- RowNumber % N -- ColumnNumber % Phi -- The Gauss matrix %% Generate Gauss matrix Phi = randn(M,N); %Phi = Phi/sqrt(M); end 2…

苏金明．2005.电子工业     1 语句末尾加 : 可以不显示到屏． who 查看变量 whos 列出变量信息 exist t 判断变量是否在空间中． help 函数 doc 函数 : doc format ;   2 数据类型 常数 : ans, eps浮点相对精度, realmax  relmin ,  pi ,   ij  虚数单位, inf 无限值, NaN 不合法值,computer计算机类型, version 版本． 变量名长度 namelengthmax 保留字 iskeywo…

matlab 初级画图 1.plot() plot(x,y)   plots each vector pairs (x,y) 画图函数画出每个点   每组变量 plot (y)   plots each vector pairs(x,y),where x=[1...n],n=length(y) 仅有一个变量时  仅画出  当x=整数点时的对应y值得点 example: plot (cos(0:pi/20:2*pi)); 从0-2π   间隔为pi/20的  那些cos值 x=0:1/1000:2…

MATLAB实例:散点密度图 作者:凯鲁嘎吉 - 博客园http://www.cnblogs.com/kailugaji/ MATLAB绘制用颜色表示数据密度的散点图 数据来源:MATLAB中“fitgmdist”的用法及其GMM聚类算法,将数据保存为gauss.txt 1. demo.m % 用颜色表示数据密度的散点图 data_load=dlmread('E:\scanplot\gauss.txt'); X=data_load(:,1:2); scatplot(X(:,1),X(:,2),'…

转:https://blog.csdn.net/cqfdcw/article/details/84995904 小波与小波包.小波包分解与信号重构.小波包能量特征提取   (Matlab 程序详解) -----暨 小波包分解后解决频率大小分布重新排列问题 本人当前对小波理解不是很深入,通过翻阅网络他人博客,进行汇总总结,重新调试Matlab代码,实现对小波与小波包.小波包分解与信号重构.小波包能量特征提取,供大家参考,后续将继续更新! 本人在分析信号的过程中发现,按照网上所述的小波包分解方法理解…

原文链接:https://blog.csdn.net/humanking7/article/details/46826105 核心提示 在Matlab中高斯滤波非常方便,主要涉及到下面两个函数: 函数: fspecial函数: imfilter   代码实现   clear all;clc;%----------------------------------------------%对图像进行高斯滤波,并显示图像%----------------------------------------…

来源: 1.https://ww2.mathworks.cn/help/images/ref/fspecial.html?searchHighlight=fspecial&s_tid=doc_srchtitle#d117e81237 2.https://www.cnblogs.com/leegod/p/8202731.html 简单的原理: 基于MATLAB的中值滤波均值滤波以及高斯滤波的实现 作者:lee神 1.   背景知识 中值滤波法是一种非线性平滑技术,它将每一像素点的灰度值设置为该点某…

前言:在地球物理勘探,流体空间分布等多种场景中,定位空间点P(x,y,x)的物理属性值Q,并绘制三维空间分布图,对我们洞察空间场景有十分重要的意义. 1. 三维立体图的基本要件: 全空间网格化 网格节点的物理属性值 2.数据准备 数据不易贴,我放在了百度网盘:点击下载数据 大概如下形式: TIP: 这里的数据矩阵为v(5276),可以看成一本27页纸,每页绘制了5*6的网格,然后27页纸叠在一起.当你理解本图绘制后,数据可以随意制作. 3.主要函数:slice.isosurface.patch…

前言:在地球物理勘探,流体空间分布等多种场景中,定位空间点P(x,y,x)的物理属性值Q,并绘制三维空间分布图,对我们洞察空间场景有十分重要的意义. 1. 三维立体图的基本要件: 全空间网格化 网格节点的物理属性值 2.数据准备 数据不易贴,我放在了百度网盘:点击下载数据 大概如下形式: TIP: 这里的数据矩阵为v(5276),可以看成一本27页纸,每页绘制了5*6的网格,然后27页纸叠在一起.当你理解本图绘制后,数据可以随意制作. 3.主要函数:slice.isosurface.patch…

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

工作需要在MATLAB中绘制质点轨迹并保存成GIF以便展示. 绘制质点轨迹动图可用comet和comet3命令,使用例子如下: t = 0:.01:2*pi;x = cos(2*t).*(cos(t).^2);y = sin(2*t).*(sin(t).^2);z = t;comet(x,y,0.1); %绘制二维%comet3(x,y,z,0.1); %绘制三维 这有一个问题在于comet或comet3无法控制绘制的动画的速度,为了实现这一点,一个可行的方法是重写comet和comet3,在每…

在安装好的matlab下有MCR(MatlabCompilerRuntime)在matlab2011/toolbox/compiler/deploy/glnxa64下找到MCRInstaller.zip安装包解压后执行install命令,选择路径进行安装即可 环境变量设置:setenv LD_LIBRARY_PATH .:$(MCR_ROOT)/v716/bin/glnxa64:$(MCR_ROOT)/v716/runtime/glnxa64:$(MCR_ROOT)/v716/sys/os/gl…

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

Atitit MATLAB 图像处理 经典书籍attilax总结 1.1. MATLAB数字图像处理1 1.2. <MATLAB实用教程(第二版)>((美)穆尔 著)[简介_书评_在线阅读] - 当当图书.html1 1.3. 数字图像处理(MATLAB版)(第二版)(本科教学版)2 1.1. MATLAB数字图像处理 第1章 图像处理与MATLAB2007a简介 第2章 图像的编码和解码 第3章 图像复原 第4章 图像处理的相关操作 第5章 图像频域变换 第6章 图像处理中的代数运算及几何变…

Atitit MATLAB 图像处理attilax总结 1.1. 下载 Matlab7.0官方下载_Matlab2012 v7.0 官方简体中文版-办公软件-系统大全.html1 1.2. Matlab cli接口 命令行执行m文件1 1.3. Matlab编程语言2 1.4. Matlab helo world3 1.5. matlab如何调用外部函数_百度知道.html3 1.6. matlab中如何定义函数3 1.7. 获取matlab的输出console3 2. qa4 2.1. Err…

Atitit java c# php c++ js跨语言调用matlab实现边缘检测等功能attilax总结 1.1. 边缘检测的基本方法Canny最常用了1 1.2. 编写matlab边缘检测代码,保存为m文件1 1.3. 通过cli接口调用Matlab执行m文件1 1.4. 效果如图1 1.1. 边缘检测的基本方法Canny最常用了 1.2. 编写matlab边缘检测代码,保存为m文件 fprintf('Hello World'); imag = imread('C:\00edge\a.jp…

 试验报告 一.试验原理: 图像点处理是图像处理系列的基础,主要用于让我们熟悉Matlab图像处理的编程环境.灰度线性变换和灰度拉伸是对像素灰度值的变换操作,直方图是对像素灰度值的统计,直方图均衡是对灰度值分布的变换. 1.灰度线性变换 (1)线性变换函数 原图向灰度值为g,通过线性函数f(x)=kx+b转换为f(g)得到灰度的线性变换. (2)代码实现 Matlab中支持矩阵作为函数参数传入,定义一个线性转换函数,利用Matlab矩阵操作,用一行代码即可对整个二维图像矩阵中所有点的灰度进行线.…

实验一图像的滤波处理 一.实验目的 使用MATLAB处理图像,掌握均值滤波器和加权均值滤波器的使用,对比两种滤波器对图像处理结果及系统自带函数和自定义函数性能的比较,体会不同大小的掩模对图像细节的影响. 二.实验内容 使用n=3,5,9的正方形均值滤波器和加权均值滤波器对图像Fig1的滤波处理.观察处理前后图像效果,分析实验结果和算法特点. 三.实验原理 1.均值滤波器(平滑线性滤波器):其响应是包含模板内像素平均值,低通滤波器. R = (领域内所有点的值*对应点的掩模系数之和)/(掩模系数之…

当我第一次用matlab语言编写一个工程项目时,发现自己编写的脚本里循环特别多,导致编程效率很低,这让我特别苦恼.有一次导师让我阅读他编写的一个Matlab脚本,并按照新要求对其进行更改.我发现脚本里多次用到find()函数,猛然间豁然开朗,原来有大部分循环可以用find()函数解决!对我而言,find()函数是我从Matlab小白开始进阶的标志,仅以此文,献给find()函数! [1] find()函数基本功能 find()函数的基本功能是返回向量或者矩阵中不为0的元素的位置索引. >> X…