c语言 传递函数,传递函数的C语言实现

传递函数的C语言实现

>> s=tf( s ) Transfer function: s >> sys=1/(3*s+1) Transfer function: 1 ------- 3 s + 1 >> bode(sys) >> c2d(sys,0.0002, tustin ) Transfer function: 3.333e-005 z + 3.333e-005 ------------------------- z - 0.9999 Sampling time (seconds): 0.0002 Y/X = 3.333e-005+3.333e-005Z(-1) ------------------------- 1 - 0.9999Z(-1) Y(1 - 0.9999Z(-1))=X(3.333e-005+3.333e-005Z(-1)) Y = X*3.333e-005X +3.333e-005X(-1)+ 0.9999Y(-1) 按照这个方程编写不对，因为系数精度太差了 >> [a b]=tfdata(ans, v ) a = 1.0e-004 * 0.2580 0.2580 b = 1.0000 -0.5481 这才是差不多的系数 >> step(sys) >> hold on >> step(ans) >> plot(y) >> hold on >> plot(test(1:50001)) // CTransFunc.cpp : Defines the entry point for the console application. // #include “stdafx.h“ #include “stdio.h“ #include “math.h“ double CTransFunc(double d) { static double [2] ={0.0,0.0}; static double Output[2]={0.0,0.0}; [0]=d; Output[0]=[0]*(0.)+(0.)*[1]+0.5481*Output[1]; Output[1]=Output[0]; [1]=[0]; return Output[0]; } int main(int argc, char* argv[]) { int i=0; FILE *pFile=NULL; pFile=(“test.txt“,“w+“); for (i=0;i