【C语言】、杨氏矩阵

本文介绍了一种在递增有序的二维数组中查找指定元素的方法,该算法的时间复杂度优于O(N)。通过从数组右上角开始遍历,根据目标值与当前值的大小关系调整搜索方向,最终定位目标值或确定不存在。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

有一个二维数组.
数组的每行从左到右是递增的,每列从上到下是递增的.
在这样的数组中查找一个数字是否存在。

时间复杂度小于O(N)

代码如下:

#include<stdio.h>

void find(int arr[3][3], int *px, int *py, int key)
{
	int row = 0;
	int col = *py - 1;
	while ((col >= 0) && (col<*py))
	{
		if (arr[row][col] == key)
{
	*px = row;
	*py = col;
	return;
	}
		else if (arr[row][col]<key)
		{
	row++;
		}
		else
		{
		col--;
		}
	}
	*px = -1;
	*py = -1;
}
 int main()
{
	int arr[3][3] = { 1,2,3,4,5,6,7,8,9 };
	int key = 1;
	int x = 3;
	int y = 3;

	find(arr, &x, &y, key);
	printf("x = %d y = %d\n", x, y);
	system("pause");
	return 0;
}
结果:



### 关于C语言实现刚度矩阵求解 对于结构力学中的有限元方法,计算刚度矩阵是一个核心部分。下面展示了一个简单的例子来说明如何利用C语言构建并求解一个基本的一维杆单元的刚度矩阵。 #### 定义全局变量与函数原型声明 ```c #include <stdio.h> #include <stdlib.h> #define E 200e9 // 杨氏模量 (Pa) #define A 0.01 // 截面积 (m²) void assemble_global_stiffness(int n, double *K); double compute_local_stiffness(double length); int main() { int num_elements = 3; // 假设有三个元素 // 动态分配内存给全局刚度矩阵 K size_t matrix_size = sizeof(double)*(num_elements+1)*(num_elements+1); double* global_K = malloc(matrix_size); // 初始化全局刚度矩阵为零 for(size_t i=0;i<(num_elements+1)*(num_elements+1);i++) *(global_K+i)=0; // 组装全局刚度矩阵 assemble_global_stiffness(num_elements, global_K); printf("Global stiffness matrix:\n"); for(int row=0;row<=num_elements;row++){ for(int col=0;col<=num_elements;col++) printf("%lf ", *((global_K)+(row*(num_elements+1))+col)); putchar('\n'); } free(global_K); // 释放动态分配的空间 } ``` #### 计算局部刚度矩阵并将它们组装到全局刚度矩阵中 ```c // 函数用于将各个单元贡献加至整体刚度矩阵上 void assemble_global_stiffness(int n, double *K){ const float element_length = 1.0 / n; for(int e=0;e<n;++e){ // 对每一个单元执行操作 double local_k = compute_local_stiffness(element_length); // 更新对应的行和列位置处的值 (*((K)+e*n+(e))) += local_k; (*((K)+e*n+(e+1))) -= local_k; (*((K)+(e+1)*n+(e))) -= local_k; (*((K)+(e+1)*n+(e+1))) += local_k; } } // 返回单个线性弹簧模型下的局部刚度系数 k = EA/L double compute_local_stiffness(double length){ return E*A/length; } ``` 上述代码展示了怎样创建一个代表整个系统的大型稀疏方阵——即所谓的“全局”刚度矩阵,并通过迭代遍历所有独立的小型子系统(在这里指的就是那些组成更大结构体的一个个小梁),逐步填充这个大矩阵[^1]。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值