
基础数学
qq-120
这个作者很懒,什么都没留下…
展开
-
矩阵形式的泰勒展开
如何求任意n维空间任意m个列向量构成的平行多面体的体积用楔形积求任意n维空间任意m个列向量构成的平行多面体的体积公式为原创 2022-07-10 18:43:30 · 2039 阅读 · 0 评论 -
平面旋转任意角度x
a,b是n维空间中的任意两个不相关的列向量,试构建一个在a,b所在的平面旋转任意角度x的矩阵?References:https://blog.youkuaiyun.com/qq_41585683/article/details/115838461原创 2022-07-10 18:39:04 · 1165 阅读 · 0 评论 -
微分方程组、无约束优化、线性规划、非线性规划、0-1规划、二次型规划
fprintf('解析解.\n');syms t x y[x,y] = dsolve('D2x=-2*x-3*Dx+exp(-5*t)',' D2y==2*x-3*y-4*Dx-4*Dy-sin(t)','x(0)=1','Dx(0)=2','y(0)=3','Dy(0)=4')x = simplify(x); y = simplify(y);fprintf('数值解.\n');[t1,x1] = ode45('myfun',[0,10],[1;2;3;4]);ezplot(x,[0,10]).原创 2022-05-05 19:02:51 · 558 阅读 · 0 评论 -
数据插值、数据拟合、最小二乘拟合、方程组解析解&迭代解法
(欠定方程)%解析解:A = [2 -9 3 -2 -1; 10 -1 10 5 0; 8 -2 -4 -6 3; -5 -6 -6 -8 -4];B = [-1 -4 0; -3 -8 -4; 0 3 3; 9 -5 3];[rank(A), rank([A B])]x0 = null(sym(A)); % 求AX=0的基础解x_analytical = sym(A)\B; syms a;x = a*[x0 x0 x0]+x_analyticalA*x - B输出:x = [ .原创 2022-05-05 18:48:19 · 803 阅读 · 0 评论 -
泰勒展开、级数求和、数值微分、数值积分
%inputsyms a x;f = exp(-5*x) * sin(3*x + pi/3);T = taylor(f,x,'Order', 5); pretty(T) % 此处展开到前 5 项T2 = taylor(f,x,a,'Order', 5); pretty(T2) % 此处展开到前 5 项%input%第一种方式syms n k; s = symsum(1/(1+k*pi/n^2),k,1,n) % 先求和limit(s/n,n,Inf) % 再求极限%第二种方式s.原创 2022-05-05 18:14:17 · 1077 阅读 · 0 评论 -
极限、隐函数导数、微分
%inputsyms x yf = (1-cos(x^2 + y^2))/((x^2 + y^2)*exp(x^2 + y^2));L = limit(limit(f,x,0),y,0)%outputL =0%inputsyms x y; f = atan(y/x)-log(x^2+y^2); pretty(-simplify(diff(f,x)/diff(f,y)))%output y + 2 x---------- 2 y + x%inputsyms t x.原创 2022-05-05 18:04:07 · 317 阅读 · 0 评论 -
解微分方程_matlab
(1)先计算右边式子的值;% Step 1syms t; u=exp(-5*t)*cos(2*t+1)+5;uu=5*diff(u,t,2)+4*diff(u,t)+2*u打印uu = 87*exp(-5*t)*cos(2*t + 1) + 92*exp(-5*t)*sin(2*t + 1) + 10(2)在计算左边式子的值;好像新版的matlab dsolve函数格式发生了一点点改变,使用dsolve函数时先help dsolve一下。% Step 2syms y(t)eqn .原创 2022-04-07 22:06:35 · 218 阅读 · 0 评论 -
数值分析与算法-喻文健-国科大
国科大数值分析资源列表注:公式太难打了,已传文件,内含亲调代码。数值分析代码第一二单元作业第三单元作业第四单元作业第五单元作业第六单元作业第七单元作业期末考试题原创 2021-11-17 11:17:44 · 2920 阅读 · 0 评论 -
《矩阵特征值、特征向量与对角化》内容小结、题型与典型题求解
转自 点这转载 2021-10-18 20:39:31 · 134 阅读 · 0 评论 -
Z逆变换(2020.10.21)
1. Z逆变换 求Z逆变换的常用方法有三种:围线积分法(留数法)、部分分式展开法和幂级数展开法(长除法)。 理论公式推导参考教材,下面就matlab内置函数进行讨论:2. Z逆变换函数 2.1 residues 2.2 residue 对拉普拉斯算子s的,适用于连续系统; 2.3 residuez 对Z变换算子z的,适用于离散系统; [r,p,k] = residuez(b,a)...原创 2020-10-21 15:19:39 · 14221 阅读 · 0 评论 -
方差分析的MATLAB实现(单因子)
单因子方差分析是比较俩组或多组数据的均值,它返回原假设——均值相等的概率,在matlab中我们可以用函数anova1进行单因子方差分析。函数anova1格式: p=anova1(X) %X的各列为彼此独立的样本观察值,其元素个数相同,p为各列均值相等的概率值,若p值接近于0,则原假设受到质疑,说明至少有一列均值与其余列均值有明显的不同p=anova1(X,group) %X和group为...原创 2019-04-23 20:29:06 · 5012 阅读 · 0 评论 -
matlab仿真随机数的产生
概率论和数理统计实验(matlab中实现)一.伯努利分布R=binornd(N,P); //N,P为二次分布的俩个参数,返回服从参数为N,P的二项分布的随机数,且N,P,R的形式相同。R=binornd(N,P,m); //m是一个1*2向量,它为指定的随机数的个数,其中N,P分别代表返回值R中行与列的维数;R=binornd(N,P,m,n); //m,n分别表示R的行数与列...原创 2019-03-01 17:27:56 · 4842 阅读 · 2 评论 -
十大概率密度函数
一.二项分布x=0:10;y=binopdf(x,10,0.4);plot(x,y,'*')二.卡方分布x=0:0.3:10;y=chi2pdf(x,4);plot(x,y)三.非中心卡方分布x=(0:0.2:10)';p1=ncx2pdf(x,3,2);p=chi2pdf(x,3);plot(x,p,'-',x,p1,'--')四.指数分布x=0:0.2:...原创 2019-03-01 21:07:08 · 11616 阅读 · 0 评论 -
mean,std,var
1、 均值数学定义: MATLAB中均值、方差、均方差的计算方法Matlab函数:meanx=[1 2 3];mean(x)ans = 2如果X是一个矩阵,则其均值是一个向量组。mean(X,1)为列向量的均值,mean(X,2)为行向量的均值。x=[1 2 3;4 5 6];mean(x,1)ans =2.5000 3.5000 4.5000若要求整个...原创 2019-04-23 23:49:46 · 3143 阅读 · 0 评论 -
二维矩阵与二维矩阵之间的卷积
最近在学习数字图像处理(Digital Image Processing,DIP)这门课,感觉有些吃力。由于在数字信号处理(Digital Singal Processing, DSP)这门课中只学了一维矩阵之间的卷积运算。第二步,需要将卷积核h的中心对准f(0,0),然后对应元素相乘后相加,没有元素的地方补0。第一步,需要将卷积核翻转180°,kernel矩阵 也就变成了。如果需要计算/了解多维的卷积,可看下面链接的博客。一定要主要conv2中A,B矩阵的顺序!假设我们的卷积核h为kernel矩阵(3。原创 2020-04-25 04:02:24 · 5905 阅读 · 2 评论