奇异值分解SVDC语言源代码

奇异值分解在某些方面与对称矩阵或Hermite矩阵基于特征向量的对角化类似。然而这两种矩阵分解尽管有其相关性，但还是有明显的不同。对称阵特征向量分解的基础是谱分析，而奇异值分解则是谱分析理论在任意矩阵上的推广。 资源提供的是奇异值分解的C语言实现。

#include “svd.h“
#include “stdio.h“
main（）
{
	int ij;
	/*double a[4][3]={{1.01.0-1.0}{2.01.00.0}{1.0-1.00.0}{-1.02.01.0}};
	double b[3][4]={{1.01.0-1.0-1.0}{2.01.00.02.0}{1.0-1.00.01.0}};
	double u[4][4]v[3][3]c[4][3]d[3][4];*/
	double a[12]={1.01.0-1.02.01.00.01.0-1.00.0-1.02.01.0};
	double b[12]={1.01.0-1.0-1.02.01.00.02.01.0-1.00.01.0};
	double u[16]v[9]c[12]d[12];

	for（i=0;i<16;i++）
		u[i]=0;
	double eps=0.000001;
	i=dluav（a43uveps5）;
	
	printf（“\n“）;
	printf（“i=%d\n“i）;
	printf（“\nMAT U Is:\n“）;
	for（i=0;i<=3;i++）
	{
		for（j=0;j<=3;j++）
			printf（“%e “u[i*4+j]）;
		printf（“\n“）;
	}
	printf（“\n“）;
	printf（“MAT V IS:\n“）;
	for（i=0;i<=2;i++）
	{
		for（j=0;j<=2;j++）
			printf（“%e “v[i*3+j]）;
		printf（“\n“）;
	}
	printf（“\n“）;
	printf（“MAT A Is:\n“）;
	for（i=0;i<=3;i++）
	{
		for（j=0;j<=2;j++）
			printf（“%e “a[i*3+j]）;
		printf（“\n“）;
	}	
	damul（ua443c）;
	damul（cv433a）;
	printf（“\nMAT UAV IS:\n“）;
	for（i=0;i<=3;i++）
	{
		for（j=0;j<=2;j++）
			printf（“%e “a[i*3+j]）;
		printf（“\n“）;
	}
	printf（“\n\n“）;
	printf（“第二个例子\n“）;
	i=dluav（b34vueps5）;
	printf（“i=%d\n“i）;
	printf（“\n“）;
	printf（“MATUS IS:\n“）;
	for（i=0;i<=2;i++）
	{
		for（j=0;j<=2;j++）
			printf（“%e “v[i*3+j]）;
		printf（“\n“）;
	}
	printf（“\n“）;
	printf（“MAT V IS:\n“）;
	for（i=0;i<=3;i++）
	{
		for（j=0;j<=3;j++）
			printf（“%e “u[i*3+j]）;
		printf（“\n“）;
	}
	printf（“\nMAT B IS:\n“）;
	for（i=0;i<=2;i++）
	{
		for（j=0;j<=3;j++）
			printf（“%e “b[i*4+j]）;
		printf（“\n“）;
	}
	printf（“\n\n“）;
	printf（“MAT UBV IS:\n“）;
	damul（vb334d）;
	damul（du344b）;
	for（i=0;i<=2;i++）
	{
		for（j=0;j<=3;j++）
			printf（“%e “b[i*4+j]）;
		printf（“\n“）;
	}
	printf（“\n“）;

}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1804  2004-06-07 09:37  svd\main.cpp

     文件      10381  2010-02-19 15:17  svd\svd.cpp

     文件       4367  2010-02-19 15:24  svd\svd.dsp

     文件        514  2010-02-19 15:08  svd\svd.dsw

     文件        573  2004-06-07 09:51  svd\svd.h

     文件      50176  2011-05-04 16:02  svd\svd.ncb

     文件      53760  2011-05-04 16:02  svd\svd.opt

     文件        240  2011-05-04 16:01  svd\svd.plg

     目录          0  2011-05-04 16:02  svd\Debug

     目录          0  2011-05-04 16:02  svd

----------- ---------  ---------- -----  ----

               121815                    10