模糊数学求传递闭包(C语言)

一、问题描述

我们在模糊数学中利用传递闭包法进行聚类分析的时候,需要求模糊相似矩阵的传递闭包,对于阶数比较高的矩阵,手工计算太过于繁琐,我们可以借助程序实现。

我们先来看一下模糊矩阵的合成规则:

X=\{ x_1,x_2,\cdots,x_m\},Y = \{ y_1,y_2,\cdots,y_s\},Z = \{ z_1,z_2,\cdots,z_n\}为有限论域,

R_1 = (a_{ik})_{m\times s},R_2 = (b_{kj})_{s\times n}

R_1R_2的合成为

R_1\circ R_2 = C=(c_{ij})_{m\times n}

其中c_{ij} = \vee_{k=1}^s(a_{ik} \wedge b_{kj})

模糊矩阵的合成类似于矩阵的乘法,不过将 乘 换成了 取小,将 加 换成了 取大。

下面我们来看一个模糊矩阵合成的例子。

已知模糊矩阵R_1 = \begin{bmatrix} 0.2 & 0.9 &0.5 \\ 0.4& 0.1 &0.8 \\ 0.6&0.7 &0.3 \end{bmatrix}R_2 = \begin{bmatrix} 0.9 & 0.5 &0.2 \\ 0.4& 0.8 &0.7 \\ 0.6& 0.3 &0.1 \end{bmatrix}

则模糊矩阵R_1\circ R_2

R_1\circ R_2 = \begin{bmatrix} 0.2 &0.9 &0.5 \\ 0.4& 0.1 &0.8 \\ 0.6&0.7 &0.3 \end{bmatrix}\circ \begin{bmatrix} 0.9 &0.5 &0.2 \\ 0.4& 0.8 &0.7 \\ 0.6& 0.3 & 0.1 \end{bmatrix}=\begin{bmatrix} 0.5 &0.8 &0.7 \\ 0.6& 0.4 &0.2 \\ 0.6& 0.7 &0.7 \end{bmatrix} 

我们使用平方法求R的传递闭包.

二、代码实现

#include 

/// 
/// 取小
/// 
/// 
/// 
/// 
double getMin(double x, double y) {
	return x <= y ? x : y;
}

/// 
/// 取大
/// 
/// 
/// 
/// 
double getMax(double x, double y) {
	return x >= y ? x : y;
}

/// 
/// 矩阵合成
/// 
/// 
/// 
/// 
/// 
void Matrix_composition(double(*s)[10], double(*t)[10], double(*p)[10], int n) {
	int i, j, k;
	double sum = 0;
	for (i = 0; i < n; i++) {
		for (j = 0; j < n; j++) {
			sum = 0;
			for (k = 0; k < n; k++) {
				sum = getMax(getMin(s[i][k], t[k][j]), sum);
			}
			p[i][j] = sum;
		}
	}
}

/// 
/// 打印
/// 
/// 
/// 
void Print(double(*s)[10], int n) {
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			printf("%.2f  ", s[i][j]);
		}
		printf("\n");
	}
}

/// 
/// 判断矩阵是否相等
/// 
/// 
/// 
/// 
void equalMatrix(double(*s)[10], double(*p)[10], int n) {
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			if (s[i][j] != p[i][j]) {
				printf("不相等!\n");
				return;
			}
		}
	}
	printf("合成前后相等!\n");
}
int main() {

	double A[10][10];
	double B[10][10];
	int flag = 1;
	printf("请输入矩阵的阶数:");
	int n = 0;
	scanf("%d", &n);
	while (flag) {
		printf("------------------------------------------------------------------\n");
		printf("请输入矩阵:\n");
		for (int i = 0; i < n; i++) {
			for (int j = 0; j < n; j++) {
				scanf("%lf", &A[i][j]);
			}
		}

		Matrix_composition(A, A, B, n);

		printf("初始矩阵为:\n");
		Print(A, n);


		printf("合成后的矩阵为:\n");
		Print(B, n);

		equalMatrix(A, B, n);
		printf("------------------------------------------------------------------\n");

		printf("是否继续计算,如果是,请输入1,否则输入0\n");
		printf("请输入:");
		scanf("%d", &flag);
	}
	return 0;
}

三、效果演示

求传递闭包

R = 
1.0  0.9  0.7  0.5  0.9  1.0  0.4
0.9  1.0  0.5  0.4  0.5  0.4  0.1
0.7  0.5  1.0  0.5  0.7  0.1  0.4
0.5  0.4  0.5  1.0  0.5  0.7  0.9
0.9  0.5  0.7  0.5  1.0  0.7  0.1
1.0  0.4  0.1  0.7  0.7  1.0  0.5
0.4  0.1  0.4  0.9  0.1  0.5  1.0

模糊数学求传递闭包(C语言)_第1张图片

 模糊数学求传递闭包(C语言)_第2张图片

模糊数学求传递闭包(C语言)_第3张图片 

我们发现直到合成第三次才有R^8 = R^4,于是R的传递闭包t(R)  = R^4. 

你可能感兴趣的:(C语言系列,矩阵,学习,笔记)