九灵猴-CSDN博客

原创 matlab向量中插入元素

c=[2 4 7 2 2 3];a=ones(1,9);b=[2 5 6];a(b)=0;index=find(a);a(index)=c(1,:)输出a=2 0 4 7 0 &nb...

2019-04-12 10:42:03 13528

转载在行向量或列向量的某个位置插入一个数

a = [2 4 6 8 12]; %已知行向量b=10; %插入的数n=5; %插入的位置aa(1,n)=b; %插入数aa(1,1:n-1)=a(1,1:n-1); %前面不变aa(1,n+1:length(a)+1)=a(1,n:length(a)) %后面的数后移一位aa =[2 4 6 8 10 12]...

2019-04-11 21:34:06 4319

转载平面桁架静力计算程序

format short %数据显示设置%从文件中读取数据fid = fopen(‘data.txt’,‘rt’);if fid == -1disp(‘数据文件打开错误！’);returnend%结构信息[N,count] = fscanf(fid,’%d’,1);%结点总数[NM,count] = fscanf(fid,’%d’,1);%单元总数[NR,count] = fs...

2019-04-10 14:50:43 2658

转载采用Newmark-β法求解振动方程——MATLAB

以四层刚架为例进行计算分析。某四层框架结构在顶部受一个简谐荷载的作用：力的作用时间t1=5.0s，F0=100N。试计算响应的时间为20s（计算分2000步完成）。阻尼矩阵由Rayleigh阻尼构造。MATLAB程序编写如下：m=[1,2,3,4];m=diag(m);k= [800 -800 0 0; -800 2400 -1600 0; 0 -1600 4800 -3200; 0 ...

2019-01-22 15:58:00 18742 17

转载 Matlab中生成向量的三种方法

方法一、直接输入法如我们要生成向量A=[7,8,9,4,5,6];这里A是个一维向量，其中的分量是7,8,9,4,5,6。方法二、“：”生成法这种方法适用于元素与元素之间存在等距步长（差值）的情况，即当元素间呈等差数列时可以使用。冒号生成法，其语法格式为：1.向量名=[第一个元素数值：步长：最后一个元素的数值]；2.向量名=第一个元素数值：步长：最后一个元素的数值如我们要生成向量A=[0,...

2019-01-22 15:08:00 62064

转载 matlab中triu（A，x）函数用法

triu–抽取上三角矩阵1、U=triu(X,k)抽取矩阵中第k条对角线及其以上的元素。k=0指主对角线，k>0指主对角线以上的第k条对角线，k<0指主对角线以下的第k条对角线。2、U=triu(X)返回矩阵的上三角矩阵，即默认k=0。...

2019-01-20 19:49:46 13559 1

转载 matlab中if语句使用方法和应用实例

1、我们首先需要了解if语句是一种选择判断语句，可以和for语句结合使用，也可以单独使用，这里创建了3个变量，作出了一个简单判断如果a小于b，b等于a，需要注意的是if语句后面要有end语句结束。运行脚本之后，在命令行窗口可以看到b已经等于a了。如下图所示：2、介绍一下if嵌套语句，这里作出的判断是如果a小于b,b等于a，如果继续b还大于c，b等于c，如下图所示：3、介绍if语句中的e...

2019-01-18 16:18:23 131556

转载 matlab中for语句使用方法和应用实例

for循环语句1、一般格式为：for x（循环变量） = array（数组）commands（执行的循环代码）end2、array可以是一个数字，也可以是数组，例如输入：for a=5for a=1:5for a=1:1:5（以1为步长到5）只不过在a=1:5和a=1:1:5时，会显示之间的结果，a=5时只显示最后结果。a的变动就是第一次循环a=1，第二次循环a=2，第三次循环a...

2019-01-18 15:39:59 72048 1

转载 matlab中size函数总结

matlab中size函数总结size(A)函数是用来求矩阵的大小的。比如说一个A是一个3×4的二维矩阵：1、size（A） %直接显示出A大小输出：ans=3 42、s=size（A）%返回一个行向量s，s的第一个元素是矩阵的行数，第二个元素是矩阵的列数输出：s=3 43、[r,c]=size（A）%将矩阵A的行数返回到第一个输出变量r，将矩阵的列数返回到第二个输出变量c输出...

2019-01-18 09:50:18 426

weixin_44563389的博客