svd分解的C语言实现

一般实矩阵的svd分解的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“）;

}

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

     文件       6748  2011-04-15 16:53  svd\svd\Debug\main.obj

     文件     230320  2011-04-15 16:53  svd\svd\Debug\main.pch

     文件      25600  2011-04-15 16:53  svd\svd\Debug\main.pdb

     文件      33792  2011-04-15 16:53  svd\svd\Debug\vc60.idb

     文件      45056  2011-04-15 16:53  svd\svd\Debug\vc60.pdb

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

     文件       3377  2011-04-15 16:52  svd\svd\main.dsp

    .......         0  2000-11-05 00:44  svd\svd\main.h

     文件      25600  2011-04-15 16:53  svd\svd\main.ncb

     文件       1146  2011-04-15 16:53  svd\svd\main.plg

     文件      10381  2004-06-07 09:50  svd\svd\svd.cpp

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

     文件      19456  2004-06-07 09:51  svd\svd\svd.ncb

     文件        895  2004-06-06 09:25  svd\svd\svd.sln

    ..A..H.     11776  2004-06-07 09:51  svd\svd\svd.suo

     文件       3443  2004-06-06 14:41  svd\svd\svd.vcproj

     目录          0  2011-04-15 16:53  svd\svd\Debug

     目录          0  2011-04-15 16:52  svd\svd

     目录          0  2011-04-15 16:52  svd

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

               420185                    20