劳斯稳定性判据程序代码

用于求解4阶及以下的劳斯判据，通过第一列数据的大小关系，判断系统的稳定性。

// judge_rosh.cpp : Defines the entry point for the console application.
//

#include “stdafx.h“
#include “stdio.h“
#include “stdlib.h“


int main（int argc char* argv[]）
{
	int ijm;
	float c[5]a[5][5]={{0}{0}{0}{0}{0}};
	float range（float a[5][5]float c[5]int k）;
	printf（“四阶及四阶以下特征方程劳斯判据判稳法\n\n输入系数时请把多余系数置零!\n“）;
	printf（“\nhow many coefficient do you want to input?\n“）;
	scanf（“%d“&m）;
	printf（“\nplease input 5 coefficient number!\n“）;
	float judge（float a[5][5]int k）;
	for（i=0;i<5;i++）
	{
		scanf（“%f“&c[i]）;
		/*rand（）;
		c[i]=rand（）%5;*/}
	/*for（i=0;i<5;i++）
	{
		printf（“%f\n“c[i]）;}
	printf（“\n“）;*/
	range（acm）;
	judge（am）;
	return 0;
}
float range（float a[5][5]float c[5]int k）
{
	int ijn;
	for（i=0;i<5;i++）
	{
		n=i;
		if（i		{
	    	for（j=0;j<3;j++）
			{
				if（i<2）
				{
            		a[i][j]=c[n];
		    		n=n+2; 
				}
   	        	else
				{
         	  	   a[i][j]=（a[i-1][0]*a[i-2][j+1]-a[i-2][0]*a[i-1][j+1]）/a[i-1][0];}
			}
		}
		else
		{
			a[i][0]=c[k-1];}
	}
	printf（“劳斯判据矩阵为:\n“）;
	for（i=0;i	{
		for（j=0;j		{
			printf（“%f  “a[i][j]）;
		    if（j==k-3）
			{printf（“\n“）;}
		}
	}
	return（a[5][5]）;
}
float judge（float a[5][5]int k）
{
	int ij;
	float sum=1;
	for（i=0;i	{
		if（a[i][0]<0）
		{
			sum=0;
			printf（“\nsystem is unstatable!\n\n“）;
			break;}
		else if（a[i][0]==0）
		{
			sum=0;
			for（j=0;j			{
				sum=sum+a[i][j];}
	    	if（sum==0）
			{
	    		printf（“\nsystem is unstatable!\n\n“）;
	   			break;}
		}
		
	}
	if（sum!=0）
	{	
		printf（“\nsystem is statable!\n\n“）;}
	return 0;
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-04-28 21:37  劳斯稳定性判据\
     目录           0  2013-04-28 21:36  劳斯稳定性判据\Debug\
     文件      208974  2013-04-28 21:36  劳斯稳定性判据\Debug\judge_rosh.exe
     文件      204136  2013-04-28 21:36  劳斯稳定性判据\Debug\judge_rosh.ilk
     文件        6853  2013-04-28 21:36  劳斯稳定性判据\Debug\judge_rosh.obj
     文件      220396  2013-04-28 15:36  劳斯稳定性判据\Debug\judge_rosh.pch
     文件      467968  2013-04-28 21:36  劳斯稳定性判据\Debug\judge_rosh.pdb
     文件        1952  2013-04-28 14:33  劳斯稳定性判据\Debug\StdAfx.obj
     文件       41984  2013-04-28 21:36  劳斯稳定性判据\Debug\vc60.idb
     文件       53248  2013-04-28 21:36  劳斯稳定性判据\Debug\vc60.pdb
     文件        1714  2013-04-28 21:36  劳斯稳定性判据\judge_rosh.cpp
     文件        3449  2013-04-28 15:36  劳斯稳定性判据\judge_rosh.dsp
     文件         528  2013-04-28 15:43  劳斯稳定性判据\judge_rosh.dsw
     文件       41984  2013-04-28 21:37  劳斯稳定性判据\judge_rosh.ncb
     文件       48640  2013-04-28 21:37  劳斯稳定性判据\judge_rosh.opt
     文件        1298  2013-04-28 21:36  劳斯稳定性判据\judge_rosh.plg
     文件        1232  2013-04-28 14:15  劳斯稳定性判据\ReadMe.txt
     文件         297  2013-04-28 14:15  劳斯稳定性判据\StdAfx.cpp
     文件         667  2013-04-28 14:15  劳斯稳定性判据\StdAfx.h