亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频

? 歡迎來到蟲蟲下載站！ | ?? 資源下載 ?? 資源專輯 ?? 關(guān)于我們
? 蟲蟲下載站

?? 征求求實矩陣特征值和相應(yīng)特征向量的c程序（可另開貼子加分）.htm

?? 數(shù)值計算方法
?? HTM
?? 第 1 頁 / 共 2 頁
字號: 
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<!-- saved from url=(0050)http://book.chinaz.com/CSDN/phppost10/php81380.htm -->
<HTML><HEAD><TITLE>征求求實矩陣特征值和相應(yīng)特征向量的c程序（可另開貼子加分）</TITLE>
<SCRIPT src="征求求實矩陣特征值和相應(yīng)特征向量的c程序（可另開貼子加分）.files/1.js"></SCRIPT>

<META http-equiv=Content-Type content="text/html; charset=GB2312">
<META content="MSHTML 6.00.2600.0" name=GENERATOR></HEAD>
<BODY bgColor=#ffffff leftMargin=0 topMargin=0 marginheight="0" marginwidth="0">
<CENTER>
<TABLE width="100%">
  <TBODY>
  <TR>
    <TD align=right>
      <TABLE cellSpacing=0 cellPadding=0 width="100%" align=right border=0>
        <TBODY>
        <TR>
          <TD colSpan=2><A 
            href="http://book.chinaz.com/CSDN/index.html">返回首頁</A> - <A 
            href="http://book.chinaz.com/CSDN/index.htm">當(dāng)前索引</A> - </TD></TR>
        <TR>
          <TD colSpan=2 height=15></TD></TR></TBODY></TABLE></TD></TR>
  <TR>
    <TD align=left>
      <TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>
        <TBODY>
        <TR>
          <TD vAlign=top align=right width=100>主　　題：</TD>
          <TD align=right width=5></TD>
          <TD><PRE>征求求實矩陣特征值和相應(yīng)特征向量的c程序（可另開貼子加分）</PRE></TD></TR>
        <TR bgColor=#666666>
          <TD colSpan=3 height=1></TD></TR>
        <TR height=10>
          <TD align=right></TD>
          <TD align=right></TD>
          <TD></TD></TR></TBODY></TABLE></TD></TR>
  <TR>
    <TD><PRE>
征求求任意實矩陣特征值和相應(yīng)特征向量的c程序。如果完成以下兩點中任意一點，相當(dāng)于完成本問題：
1、將任意實矩陣初等相似變換為實對稱陣；或
2、求出齊次線性方程組無數(shù)解中的任意一組。
求求各位大哥哥大姐姐，拉小弟一把。最好附有c,c++的源程序。
</PRE></TD></TR>
  <TR>
    <TD>
      <TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>
        <TBODY>
        <TR>
          <TD width=10 bgColor=#f1f5fa></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#3366cc>回復(fù)人：<B> 
            saint001(saint001) </B></FONT><FONT color=#000000>( ) 信譽(yù)：105 
          </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>2002-09-18 
            10:45:04Z </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>得分:</FONT><FONT 
            color=#000000>100</FONT> </TD></TR>
        <TR>
          <TD bgColor=#0066cc colSpan=4 height=1></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR>
        <TR>
          <TD><FONT color=#ffffff>?</FONT></TD>
          <TD colSpan=3><PRE>
我現(xiàn)在用Numerical Recipes in C程序庫做了求矩陣特征值的程序,不過特征向量還沒有求出來
這個c的程序庫感覺要比matlab用的fortran的那幾個好用,畢竟是c的,并且有pdf格式的Numerical Recipes in C這本書下載,非常好
但功能不是很強(qiáng),比如化為Hessenburg陣后它的子程序只能求出特征值,而沒有給出特征向量
對于對稱矩陣,它給出的求解特征值和特征向量的方法還是Jacobi法
</PRE></TD></TR>
        <TR>
          <TD align=right colSpan=4 height=10><FONT color=#ff6633><B><A 
            href="http://book.chinaz.com/CSDN/phppost10/php81380.htm#top"><FONT 
            color=#3366cc>Top</FONT></A></B></FONT></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR></TBODY></TABLE>
      <TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>
        <TBODY>
        <TR>
          <TD width=10 bgColor=#f1f5fa></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#3366cc>回復(fù)人：<B> 
            saint001(saint001) </B></FONT><FONT color=#000000>( ) 信譽(yù)：105 
          </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>2002-09-18 
            11:11:53Z </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>得分:</FONT><FONT 
            color=#000000>0</FONT> </TD></TR>
        <TR>
          <TD bgColor=#0066cc colSpan=4 height=1></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR>
        <TR>
          <TD><FONT color=#ffffff>?</FONT></TD>
          <TD colSpan=3><PRE>
先把求特征值的主程序發(fā)上來
#include &lt;math.h&gt;
#include &lt;alloc.h&gt;
#include "c\elmhes.c"
#include "c\hqr.c"
#define N 3

void main()
{
	double b[N+1][N+1]={{0,0,0,0},{0,3,1,-1},{0,2,2,-1},{0,2,2,0}};/*這里矩陣無需對稱,若對稱,則特征值均為實數(shù)*/
	double **a;
	double wr[N+1],wi[N+1];/*特征值的實部和虛部*/
	int i,j;
	a=(double **) malloc((unsigned) (N+1)*sizeof(double*));
	for(i=0;i&lt;=N;i++)
		a[i]=b[i]; /*賦值完畢*/

	printf("Original A:\n");
    for(i=1;i&lt;=N;i++)
	{
		for(j=1;j&lt;=N;j++)
			printf("%f  ",a[i][j]);
		printf("\n");
	}


	elmhes(a,N); /*變換至上Hessenburg型,特征值不變*/
	printf("After Elmhes:\n");
	for(i=1;i&lt;=N;i++)
		for(j=1;j&lt;(i-1);j++)
			a[i][j]=0;

	for(i=1;i&lt;=N;i++)
	{
		for(j=1;j&lt;=N;j++)
			printf("%f  ",a[i][j]);
		printf("\n");
	}

	hqr(a,N,wr,wi);/*用QR方法求上Hessenburg型矩陣的特征值*/
	printf("Eigenvalue:\n");
	for(i=1;i&lt;=N;i++)
		printf("%f + %f * i\n",wr[i],wi[i]);


</PRE></TD></TR>
        <TR>
          <TD align=right colSpan=4 height=10><FONT color=#ff6633><B><A 
            href="http://book.chinaz.com/CSDN/phppost10/php81380.htm#top"><FONT 
            color=#3366cc>Top</FONT></A></B></FONT></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR></TBODY></TABLE>
      <TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>
        <TBODY>
        <TR>
          <TD width=10 bgColor=#f1f5fa></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#3366cc>回復(fù)人：<B> 
            saint001(saint001) </B></FONT><FONT color=#000000>( ) 信譽(yù)：105 
          </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>2002-09-18 
            11:14:04Z </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>得分:</FONT><FONT 
            color=#000000>0</FONT> </TD></TR>
        <TR>
          <TD bgColor=#0066cc colSpan=4 height=1></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR>
        <TR>
          <TD><FONT color=#ffffff>?</FONT></TD>
          <TD colSpan=3><PRE>
其中(修改過的)c\elmhes.c為
#include &lt;math.h&gt;

#define SWAP(g,h) {y=(g);(g)=(h);(h)=y;}

void elmhes(a,n)
double **a;
int n;
{
	int m,j,i;
	double y,x;

	for (m=2;m&lt;n;m++) {
		x=0.0;
		i=m;
		for (j=m;j&lt;=n;j++) {
			if (fabs(a[j][m-1]) &gt; fabs(x)) {
				x=a[j][m-1];
				i=j;
			}
		}
		if (i != m) {
			for (j=m-1;j&lt;=n;j++) SWAP(a[i][j],a[m][j])
			for (j=1;j&lt;=n;j++) SWAP(a[j][i],a[j][m])
		}
		if (x) {
			for (i=m+1;i&lt;=n;i++) {
				if (y=a[i][m-1]) {
					y /= x;
					a[i][m-1]=y;
					for (j=m;j&lt;=n;j++)
						a[i][j] -= y*a[m][j];
					for (j=1;j&lt;=n;j++)
						a[j][m] += y*a[j][i];
				}
			}
		}
	}
}
</PRE></TD></TR>
        <TR>
          <TD align=right colSpan=4 height=10><FONT color=#ff6633><B><A 
            href="http://book.chinaz.com/CSDN/phppost10/php81380.htm#top"><FONT 
            color=#3366cc>Top</FONT></A></B></FONT></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR></TBODY></TABLE>
      <TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>
        <TBODY>
        <TR>
          <TD width=10 bgColor=#f1f5fa></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#3366cc>回復(fù)人：<B> 
            baizhao(云散) </B></FONT><FONT color=#000000>( ) 信譽(yù)：100 </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>2002-09-19 
            20:37:37Z </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>得分:</FONT><FONT 
            color=#000000>0</FONT> </TD></TR>
        <TR>
          <TD bgColor=#0066cc colSpan=4 height=1></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR>
        <TR>
          <TD><FONT color=#ffffff>?</FONT></TD>
          <TD colSpan=3><PRE>
To Saint001:
我的情況可你一樣，我現(xiàn)在已經(jīng)能夠很好地把任意實矩陣化為Hessenburg矩陣，然后求出特征值，但是，求不出特征向量，因為實際上這是求齊次線性方程無數(shù)解中的一個，我就差這一步。
</PRE></TD></TR>
        <TR>
          <TD align=right colSpan=4 height=10><FONT color=#ff6633><B><A 
            href="http://book.chinaz.com/CSDN/phppost10/php81380.htm#top"><FONT 
            color=#3366cc>Top</FONT></A></B></FONT></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR></TBODY></TABLE>
      <TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>
        <TBODY>
        <TR>
          <TD width=10 bgColor=#f1f5fa></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#3366cc>回復(fù)人：<B> 
            saint001(saint001) </B></FONT><FONT color=#000000>( ) 信譽(yù)：105 
          </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>2002-09-20 
            08:26:23Z </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>得分:</FONT><FONT 
            color=#000000>0</FONT> </TD></TR>
        <TR>
          <TD bgColor=#0066cc colSpan=4 height=1></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR>
        <TR>
          <TD><FONT color=#ffffff>?</FONT></TD>
          <TD colSpan=3><PRE>
全部做出來了
整個程序在
http://wwwbbs.chinaren.com/read_post.php3?pid=144473&amp;bid=2005167&amp;s=0&amp;bname=user_saint001
不知能不能看到
我這幾天要數(shù)末
</PRE></TD></TR>
        <TR>
          <TD align=right colSpan=4 height=10><FONT color=#ff6633><B><A 
            href="http://book.chinaz.com/CSDN/phppost10/php81380.htm#top"><FONT 
            color=#3366cc>Top</FONT></A></B></FONT></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR></TBODY></TABLE>
      <TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>
        <TBODY>
        <TR>
          <TD width=10 bgColor=#f1f5fa></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#3366cc>回復(fù)人：<B> 
            saint001(saint001) </B></FONT><FONT color=#000000>( ) 信譽(yù)：105 
          </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>2002-09-20 
            10:34:53Z </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>得分:</FONT><FONT 
            color=#000000>0</FONT> </TD></TR>
        <TR>
          <TD bgColor=#0066cc colSpan=4 height=1></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR>
        <TR>
          <TD><FONT color=#ffffff>?</FONT></TD>
          <TD colSpan=3><PRE>
vc++編譯的
不過發(fā)上去連個格式都沒有了
不太好看
另外可能用指針的地方可以優(yōu)化
我不再做了
</PRE></TD></TR>
        <TR>
          <TD align=right colSpan=4 height=10><FONT color=#ff6633><B><A 
            href="http://book.chinaz.com/CSDN/phppost10/php81380.htm#top"><FONT 
            color=#3366cc>Top</FONT></A></B></FONT></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR></TBODY></TABLE>
      <TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>
        <TBODY>
        <TR>
          <TD width=10 bgColor=#f1f5fa></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#3366cc>回復(fù)人：<B> 
            baizhao(云散) </B></FONT><FONT color=#000000>( ) 信譽(yù)：100 </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>2002-09-21 
            18:50:11Z </FONT></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#000000>得分:</FONT><FONT 
            color=#000000>0</FONT> </TD></TR>
        <TR>
          <TD bgColor=#0066cc colSpan=4 height=1></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR>
        <TR>
          <TD><FONT color=#ffffff>?</FONT></TD>
          <TD colSpan=3><PRE>
謝謝，我已經(jīng)看到了，不知能不能用，我試一下，不管如何，我過兩天就把分給你。對了，你是否需要建模的資料，我參加過大學(xué)、全國、美國的建模競賽，有一些積累，有興趣發(fā)郵件給我:baizhao@263.net，以后多交流。
再次感謝。
                                         baizhao
</PRE></TD></TR>
        <TR>
          <TD align=right colSpan=4 height=10><FONT color=#ff6633><B><A 
            href="http://book.chinaz.com/CSDN/phppost10/php81380.htm#top"><FONT 
            color=#3366cc>Top</FONT></A></B></FONT></TD></TR>
        <TR>
          <TD colSpan=4 height=5></TD></TR></TBODY></TABLE>
      <TABLE cellSpacing=0 cellPadding=0 width="100%" border=0>
        <TBODY>
        <TR>
          <TD width=10 bgColor=#f1f5fa></TD>
          <TD bgColor=#f1f5fa height=25><FONT color=#3366cc>回復(fù)人：<B> 
            saint001(saint001) </B></FONT><FONT color=#000000>( ) 信譽(yù)：105 
          </FONT></TD>

?? 快捷鍵說明

復(fù)制代碼 Ctrl + C
搜索代碼 Ctrl + F
全屏模式 F11
切換主題 Ctrl + Shift + D
顯示快捷鍵 ?
增大字號 Ctrl + =
減小字號 Ctrl + -
亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频
视频在线观看91| 国产成人精品免费视频网站| 久久夜色精品一区| 精品国产123| 91久久精品一区二区三区| 国产精品一线二线三线精华| 亚洲成av人片在线| 中文字幕一区二区三| 亚洲精品第1页| 欧美日韩精品欧美日韩精品一| 国模一区二区三区白浆| 玉米视频成人免费看| 久久精品夜夜夜夜久久| 91麻豆精品国产无毒不卡在线观看 | 日韩天堂在线观看| 欧洲一区在线观看| 色综合久久88色综合天天免费| 国产精品亚洲午夜一区二区三区| 视频一区欧美日韩| 亚洲亚洲精品在线观看| |精品福利一区二区三区| 欧美国产成人精品| 欧美精品一区二区精品网| 91麻豆精品国产91久久久更新时间| 色偷偷久久一区二区三区| av不卡一区二区三区| 国产99久久久国产精品潘金网站| 91影院在线观看| 五月激情六月综合| 亚洲va欧美va人人爽| 一区二区在线观看免费| 国产精品萝li| 国产精品久久久久国产精品日日| 久久久午夜精品理论片中文字幕| 亚洲精品一区二区三区福利 | 久久影院视频免费| 精品国产乱码久久久久久牛牛| 日韩精品最新网址| 欧美r级电影在线观看| 欧美精品一区二| 久久日一线二线三线suv| 26uuu国产一区二区三区| 精品免费一区二区三区| 久久日韩粉嫩一区二区三区| 国产亚洲自拍一区| 中文字幕亚洲精品在线观看| 亚洲欧美在线另类| 一区二区三区四区高清精品免费观看| 亚洲欧美aⅴ...| 亚洲444eee在线观看| 日本亚洲一区二区| 韩国在线一区二区| fc2成人免费人成在线观看播放 | 亚洲高清免费一级二级三级| 亚洲精品水蜜桃| 亚洲123区在线观看| 爽爽淫人综合网网站| 久久成人综合网| 国产精品自在欧美一区| 99re在线精品| 欧美日韩国产a| 亚洲精品一区二区三区在线观看| 中文字幕欧美三区| 夜夜爽夜夜爽精品视频| 青青国产91久久久久久| 国产精品主播直播| 色噜噜狠狠色综合中国| 日韩一区二区麻豆国产| 国产亚洲视频系列| 亚洲图片一区二区| 国产一区二区精品在线观看| 99久久久无码国产精品| 7777女厕盗摄久久久| 国产欧美中文在线| 亚洲国产精品久久不卡毛片| 精品在线一区二区三区| 99久久精品国产网站| 欧美久久一二区| 国产欧美日本一区二区三区| 亚洲午夜久久久久久久久电影网| 久久99精品久久久久| 99国产精品视频免费观看| 在线成人av影院| 国产精品欧美久久久久一区二区| 婷婷国产v国产偷v亚洲高清| 国产精品一二三四五| 欧美婷婷六月丁香综合色| 2023国产一二三区日本精品2022| 一区二区三区精品在线观看| 久久av中文字幕片| 欧美午夜精品久久久久久孕妇 | 国产精品久久久久久久久晋中 | 欧美日韩中文字幕一区二区| 国产欧美精品一区二区三区四区| 亚洲主播在线观看| 国产成人免费av在线| 91精品免费观看| 亚洲视频一区在线| 国内外成人在线| 欧美三级在线视频| 国产精品高潮久久久久无| 蜜桃av噜噜一区| 欧洲精品一区二区| 中文字幕五月欧美| 国内欧美视频一区二区| 91.麻豆视频| 伊人一区二区三区| 成人免费av资源| 精品欧美久久久| 日韩av中文字幕一区二区| 在线看日韩精品电影| 中文字幕欧美区| 国产乱码精品一区二区三| 777亚洲妇女| 亚洲成av人片在www色猫咪| 91农村精品一区二区在线| 国产欧美1区2区3区| 欧美在线综合视频| 亚洲视频你懂的| 国产a精品视频| 久久久www成人免费毛片麻豆| 美女一区二区在线观看| 欧美丰满嫩嫩电影| 亚洲成人综合在线| 欧美日韩精品欧美日韩精品| 一区二区三区四区视频精品免费| 99国产麻豆精品| 亚洲欧洲一区二区三区| 不卡免费追剧大全电视剧网站| 久久网站热最新地址| 国内精品免费在线观看| 精品电影一区二区| 激情图片小说一区| 2023国产精品自拍| 国产不卡一区视频| 国产精品国产三级国产专播品爱网 | 中文字幕日韩av资源站| av在线免费不卡| 亚洲视频免费观看| 日本道色综合久久| 亚洲国产视频一区| 51精品久久久久久久蜜臀| 日韩成人一区二区| 日韩精品一区在线| 国产一区二区三区在线观看免费视频 | 国产精品久久久久一区| 91亚洲精品久久久蜜桃| 一区二区三区鲁丝不卡| 欧美日韩精品二区第二页| 视频在线观看一区| 欧美成人a在线| 国产风韵犹存在线视精品| 国产精品嫩草99a| 在线视频你懂得一区二区三区| 午夜伊人狠狠久久| 精品88久久久久88久久久| 国产盗摄一区二区| 亚洲欧洲制服丝袜| 欧美日韩黄视频| 精品一区二区三区免费视频| 日本一区免费视频| 91国产成人在线| 日韩不卡一区二区三区 | xvideos.蜜桃一区二区| 国产99久久久久| 亚洲免费在线观看视频| 4438亚洲最大| 成人avav影音| 性感美女极品91精品| 2021久久国产精品不只是精品| av不卡在线播放| 日韩和欧美一区二区| 欧美国产成人精品| 欧美精三区欧美精三区| 国产精品1区二区.| 亚洲一二三四在线观看| 久久综合九色综合欧美就去吻| 波多野结衣欧美| 麻豆国产欧美一区二区三区| 亚洲欧美色图小说| 日韩午夜av电影| 91极品视觉盛宴| 91精品国产综合久久福利| 国产91富婆露脸刺激对白| 亚洲一区在线电影| 国产欧美日韩亚州综合| 欧美精品精品一区| av中文字幕一区| 久久精品国产99| 亚洲国产精品自拍| 国产精品女同一区二区三区| 欧美一级片在线看| 91成人在线观看喷潮| 国产成人综合网站| 免费高清在线一区| 亚洲一区二区在线免费观看视频| 久久久噜噜噜久噜久久综合| 欧美伦理电影网| 99精品视频在线播放观看| 国产毛片一区二区|