一、“近似”的两种分类 一个复杂的函数,可以通过一系列的“基底函数”(base function)的组合来近似,也就是函数逼近,有两种典型的方法: 基于全域的逼近,如傅立叶级数展开; 基于子域的分段函数组合,如有限元方法。 第一种函数逼近方式,就是力学分析中经 ...
传递矩阵法求简支梁固有频率的近似解 clc clear syms p q Sp sym x 点传递矩阵 Sf sym 场传递矩阵 n input 输入划分单元数: S Sf Sp n Sf 两端支座之间的传递矩阵 求固有频率 xs solve S , S , S , S , xs sort double xs xt xs n xt sqrt xt xe :n pi :n . 精确解 xe xe ...
2016-08-03 18:45 0 1539 推荐指数:
一、“近似”的两种分类 一个复杂的函数,可以通过一系列的“基底函数”(base function)的组合来近似,也就是函数逼近,有两种典型的方法: 基于全域的逼近,如傅立叶级数展开; 基于子域的分段函数组合,如有限元方法。 第一种函数逼近方式,就是力学分析中经 ...
可参考文涛,基于Matlab语言的多自由度振动系统的固有频率及主振型计算分析,2007 对于无阻尼系统 [VEC,VAL]=eig(inv(A)*K) 对于有阻尼系统,参考振动论坛计算程序 输入M,D,K function [v,w,zeta]=vbr_sf(m,d,k ...
上图是《有限元分析基础教程》中的图。 这是《材料力学》(孙训方)里面给出的图。 之所以给出这两幅图,是因为在推导公式的时候,第一幅图让我误解了:红箭头标注的微端中的外荷载$\bar ...
matlab矩阵求逆矩阵 因为 所以该矩阵可逆,根据 ,其中 得到 计算矩阵A每个元素的代数余子式 ...
求矩阵的模: function count = juZhenDeMo(a,b) [r,c] = size(a);%求a的行列 [r1,c1] = size(b);%求b的行列 count = 0; for j=1:r-r1+1%所求的行数中取 for i=1:c-c1+1%所有的列数中 ...
行和 sum(a, 2) 列和 sum(a) 所有元素之和 sum(sum(a)) 某一列元素之和 sum(a(:,1)) %a矩阵的第一列之和 某一行元素之和 sum(a(1,:)) %a矩阵的第一行之和 ...
在写MATLAB代码时经常会用到求某个矩阵长度,我之前通常用size,后来发现别人喜欢用length。这篇文章旨在统计分析各大求长度指令的异同点。 size(a) -- 表矩阵每个维度的长度 结果是[3,4] length(a) -- 表矩阵的最大长度,相当于max(size ...
原始矩阵 \[\begin{bmatrix} 1 & 3 & 4 \\ 5 & 6 & 9 \\ 2 & 7 & 8 \\ \end{bmatrix} \] 求一阶导数 \[\begin{bmatrix} 4 & ...