打印杨氏矩阵

这段代码实现了一个杨辉三角的生成算法,根据输入的行数n,打印出n行的杨辉三角形。每一行的元素由上一行的元素通过特定规则计算得出:当前元素等于上一行的前一个元素加上上一行的当前元素。例如,输入5将打印出前5行的杨辉三角形。

1
1 1
1 2 1
1 3 3 1
1 4 6 4 1

思路:
杨氏矩阵第一列或者当行和列相等时,元素都为1;
其余元素等于它上一行元素与其前一个元素之和;

#define  _CRT_SECURE_NO_WARNINGS
#include<stdio.h>
#include<stdlib.h>
int main(){
	int i = 0;
	int j = 0;
	int n = 0;
	scanf("%d", &n);
	int arr[100][100] = { 0 };
	for (i = 0; i <n; i++){
		for (j = 0; j <= i; j++){
			if (j == 0 || i == j){
				arr[i][j] = 1;
			}
			else{
				arr[i][j] = arr[i - 1][j - 1] + arr[i - 1][j];
			}printf("%d ", arr[i][j]);

		}printf("\n");
	}
	system("pause");
	return 0;
}

输入要打印的行数:
在这里插入图片描述
在这里插入图片描述

在三维梁单元的有限元计算中,我们需要进行坐标转换,将局部坐标系下的应力、应变、刚度矩阵等转换到全局坐标系下。坐标转换矩阵(Coordinate Transformation Matrix)在这个过程中起到了关键作用。以下是 MATLAB 中三维梁单元的坐标转换矩阵示例代码: ```matlab %% 定义局部坐标系下的梁单元的节点坐标(三维空间) node = [ 0, 0, 0; 0, 0, 1; 0, 1, 0; 1, 0, 0; 1, 1, 0; 1, 0, 1; 0, 1, 1; 1, 1, 1 ]; %% 定义局部坐标系下的梁单元的单元拓扑(即节点编号) elem = [ 1, 2; 1, 3; 1, 4; 2, 7; 2, 6; 3, 7; 3, 5; 4, 6; 4, 5; 5, 8; 6, 8; 7, 8 ]; %% 定义梁的长度和杨氏模量 L = 1; % 梁的长度 E = 2e11; % 杨氏模量 %% 计算局部坐标系下的梁单元的刚度矩阵 k_local = zeros(6,6); for i = 1:size(elem,1) x1 = node(elem(i,1),1); y1 = node(elem(i,1),2); z1 = node(elem(i,1),3); x2 = node(elem(i,2),1); y2 = node(elem(i,2),2); z2 = node(elem(i,2),3); L_e = sqrt((x2-x1)^2 + (y2-y1)^2 + (z2-z1)^2); % 计算单元长度 k_e = E * [1,-1;-1,1] * [A/L_e, -A/L_e; -A/L_e, A/L_e]; % 计算单元刚度矩阵 k_local([2*i-1,2*i],[2*i-1,2*i]) = k_local([2*i-1,2*i],[2*i-1,2*i]) + k_e; % 将单元刚度矩阵组装到总刚度矩阵中 end %% 计算局部坐标系到全局坐标系的转换矩阵 theta = pi/4; % 定义旋转角度(弧度) Rz = [cos(theta), -sin(theta), 0; sin(theta), cos(theta), 0; 0, 0, 1]; % 绕 z 轴旋转的旋转矩阵 Ry = [cos(theta), 0, sin(theta); 0, 1, 0; -sin(theta), 0, cos(theta)]; % 绕 y 轴旋转的旋转矩阵 Rx = [1, 0, 0; 0, cos(theta), -sin(theta); 0, sin(theta), cos(theta)]; % 绕 x 轴旋转的旋转矩阵 R = Rz * Ry * Rx; % 组合三个旋转矩阵得到总的旋转矩阵 T = [R, zeros(3,3); zeros(3,3), R]; % 构造局部坐标系到全局坐标系的转换矩阵 %% 计算全局坐标系下的梁单元的刚度矩阵 k_global = T' * k_local * T; %% 打印全局坐标系下的梁单元的刚度矩阵 disp('Global stiffness matrix: '); disp(k_global); ``` 上述代码中,我们首先定义局部坐标系下的梁单元的节点坐标和单元拓扑,然后计算局部坐标系下的梁单元的刚度矩阵。接着,我们定义旋转角度并计算局部坐标系到全局坐标系的转换矩阵 T。最后,将局部坐标系下的刚度矩阵通过转换矩阵 T 转换到全局坐标系下,得到全局坐标系下的梁单元的刚度矩阵 k_global。
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值