参考代码:(共三个文件)
主函数:main.m
clc
clear all
close all
tic
x = 1:3; % 自变量,请根据实际情况自行修改
r = zeros(3,3); % 组装后的整体刚度矩阵初始化,具体矩阵大小请自行修改
for i=1:3
y = m1(x,i); % 第一个函数,模拟计算单元刚度矩阵
r(i,:) = m2(y); % “另一个函数”输出组装入整体刚度矩阵
end
disp(r);
toc
第一个函数m1.m
function y = m1(x,i)
y = x.^i;
end
第二个函数m2.m
function r = m2(y) % 题干所述“另一个函数”
r = y*2; % 拟定函数,请自行修改
end
输出结果
2 4 6
2 8 18
2 16 54
------------------------------------
备注:以上方法适用于要拼装矩阵数已知的情况
还有一种矩阵拼装的方法为:
设有两个矩阵A=[1 2 3], B= [4 5 6]
则[A B]=[1 2 3 4 5 6],[A; B]=[1 2 3; 4 5 6]
(也可按照这个思路编程,适合于要拼装矩阵数未知的情况)