elapsetor

摘要： # MATLAB远程桌面打开失败解决办法 本文摘录自文章 MATLAB远程桌面打开失败解决办法 打开MATLAB安装目录，找到licenses文件夹，用记事本打开license_standalone.lic文件，在每一行加上 TS_OK，即可。 INCREMENT Aerospace_Blockse 阅读全文

摘要： 矩阵SVD奇异值分解用于图像压缩 clc; clear; I = imread('logo.png'); % 读入需要压缩的照片 Igray=rgb2gray(I);%转为灰度图像 [m,n] = size(Igray);%获取图像行列 k = 25;%设定压缩比率，一般在25-100 Igray 阅读全文

摘要： 保存Matlab图片的设置 用plot或者其他命令画图，之后保存图片 保存下来有白边，每次需要裁剪，可以用如下命令处理 f = figure; f.Position = [100,100,600,600]; %窗口位置大小 plot(x,y,'b'); grid off; set(gca,'loos 阅读全文

摘要： #Matlab无穷上限数值积分的近似方法 经常遇到从0积分到无穷的数值积分问题 首先要求被积函数有无穷限积分 数值积分可直接采用matlab 的 integral 有时候积分过程计算很复杂，可以取0到M（M是较大的数）的积分值代替无穷限积分 代码如下 function Isum = quadToIn 阅读全文

摘要： 导热问题有限元刚度矩阵二维平面四节点 二维平面4节点等参有限元的热传导刚度矩阵计算很简单，这里直接给了计算结果 function KE = elementMatVec2D(D0) KE1 = D0(1)*[2 -2 -1 1;-2 2 1 -1;-1 1 2 -2;1 -1 -2 2]; KE2 = 阅读全文

摘要： 导热问题有限元刚度矩阵三维立方体8节点 三维立方体8节点等参有限元的热传导刚度矩阵计算很复杂，这里直接给了计算结果 采用求导积分方法的matlab程序如下，DH是xyz三个方向的导热系数[kx,ky,kz] function KE = elementMatVec3DNew(DH) a=0.5; b= 阅读全文

摘要： matlab导出动画gif 获取当前绘图窗口，一帧一帧的追加写入Gif图片即可 % Write to the GIF File 追加Gif图片 frame = getframe(h); % 获得当前画图窗口h im = frame2im(frame); [imind,cm] = rgb2ind(im 阅读全文

摘要： 三张牌型概率模拟 模拟三张牌的牌形可能，不计（KA2）的顺子。 模拟1000万次的结果如下 单牌概率74.39%。豹子和顺金概率极低。 运行过程 clc;clear; cnt = 0; maxIter = 1e7; Tdata = zeros(maxIter,1); for i=1:maxIter 阅读全文

摘要： 极坐标三维插值问题 介绍 插值是在离散数据的基础上补插连续函数，使得这条连续曲线通过全部给定的离散数据点。 Matlab中有interp1,interp2函数通过给定一维或二维的点来插值。插值方法有近邻、线性、三次样条、多项式等等。 极坐标插值 ”三种笛卡尔坐标到极坐标转换插值算法比较“这篇博客说了 阅读全文

摘要： LU分解求线性方程组 解一维平板非稳态导热隐式格式时，需要求解线性方程组。LU分解适合线性方程组有唯一解的小规模求解。 function x=LUsolve(A,B,x0,eps,M) %LU分解法求方程组的解（矩阵公式求解） %A为方程组的系数矩阵；b为方程组的右端项 %x为线性方程组的解；x0为 阅读全文