tic; % this method is transform from Ritz method %is used for solving two point BVP %this code was writen by HU.D.dong in February 11th 2017 %MATLAB 7.0 clear clc N=50; h=1/N; X=0:h:1; f=inline('pi^2/2*sin(pi/2*x)'); %以下是右端向量: for i=2:N fun1=@(x) f(X

tic; % this method is transform from Galerkin method %also call it as finit method %is used for solving two point BVP which is the first and second term. %this code was writen by HU.D.dong in February 11th 2017 %MATLAB 7.0 clear; clc; N=50; h=1/N; X=

tic; clear clc N=; M=*N; h1=/M; h2=/N; x=:h1:; y=:h2:; fun=inline('exp(x)*sin(pi*y)','x','y'); f=inline('(pi^2-1)*exp(x)*sin(pi*y)','x','y'); lamda1=inline(','y'); lamda2=inline('2*y','y'); lamda3=inline('2*x','x'); lamda4=inline('x^2','x'); kesai1=i

差分格式脚本文件: tic; clear clc M=32;%x的步数 N=16;%y的步数 h1=1/M;%x的步长 h2=1/N;%y的步长 x=0:h1:1; y=0:h2:1; u=zeros(M+1,N+1);%给数值解分配内存单元 U=u;%给精确解分配内存单元 u(1,:)=y.^3;%y边值 u(M+1,:)=1+y.^3;%y边值 u(:,1)=x.^3;%x边值.uo u(:,N+1)=1+x.^3;%x边值.un for i=1:M+1 for j=1:N+1 Accura

曾攀老师的<有限元分析基础教程>第三章有二维杆单元的推导,并结合一个例题进行了解析解和基于Matlab的程序求解.但是我感觉书中的MATLAB代码有点罗嗦,而且一些实现方法也比较麻烦,比如已经知道了杆单元的起点和终点坐标,仍然需要另外给出单元局部坐标与整体坐标的夹角,这完全没必要.于是我就用Python重构了这段程序,当然并不是把书中的MATLAB代码翻译成python(事实上完全可以这么干,而且很快!).比如我使用了面向对象的思想,把杆单元写成了一个类,这样思路比较清晰. #! /usr/b

之前转载了一篇博客http://blog.sina.com.cn/s/blog_6163bdeb0102dvay.html,讲Matlab网格划分程序Distmesh,看了看程序,感觉程序写得有很多值得学的地方,所以又自己重新看了一看,加了一些注释,最后再总结一下学到的东西吧. 源代码的地址已经改变,如下http://people.sc.fsu.edu/~jburkardt/m_src/distmesh/distmesh.html.程序给出了20个划分的例子,文件名为p1_demo.m~p20_

至于工具箱的安装说明参见:http://www.matlabsky.com/thread-120-1-1.html Maplesoft<Maple Toolbox for MATLAB> http://www.matlabsky.com/thread-236-1-1.html Sergiy Iglin<Graph Theory Toolbox>(图论工具箱)http://www.matlabsky.com/thread-295-1-1.html Koert Kuipers<B

p47.(实习题-李荣华)用线性元求下列边值问题的数值解 tic; % this method is transform from Galerkin method %also call it as finit method %is used for solving two point BVP which is the first and second term. % %MATLAB 7.0 clear; clc; N=; h=/N; X=:h:; f=inline('(0.5*pi^2)*sin

1.平面操作工具箱 http://cathy.ijs.si/~leon/planman.html 2.SimMechanics 工具箱 (这个好像不是免费的) http://www.mathworks.com/access/helpdesk/help/toolbox/physmod/mech/mech.shtml 3.gaot工具箱(遗传算法工具箱) http://www.dytrol.com/viewFile.asp?Boardid=11&ID=263 4.Optimization Toolb

思路 粘弹性边界因为能够考虑地基辐射阻尼而使得结构抗震的计算结果更趋于合理,所以在需要考虑结构地基相互作用的结构抗震计算时,是较为常用的地基边界处理和地震荷载施加方法.而ABAQUS软件是经常用来进行结构响应分析的有限元软件.下面介绍一种在ABAQUS中实现粘弹性边界及地震荷载施加的方法. 粘弹性边界是通过在有限元模型的地基边界节点上施加弹簧阻尼器实现的,在ABAQUS中的实现有以下几种方法:第一种,通过ABAQUS自有的弹簧单元spring单元和阻尼单元dashpot实现,具体的单元参数可以参

小白往往听到微分方程就觉得害怕,其实数学建模中的微分方程模型不仅没那么复杂,而且很容易写出高水平的数模论文. 本文介绍微分方程模型边值问题的建模与求解,不涉及算法推导和编程,只探讨如何使用 Python 的工具包,零基础求解微分方程模型边值问题. 通过 3个 BVP 案例层层深入,手把手教你搞定微分方程边值问题. 欢迎关注『Python小白的数学建模课 @ Youcans』系列,每周持续更新 1. 常微分方程的边值问题(BVP) 1.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