傳遞函數(shù)的C語(yǔ)言實(shí)現(xiàn)

1、>>s=tf('s')Transferfunction:s>>sys=1/(3*s+1)Transferfunction:13s+1>>bode(sys)/Figure1匚舊反£ileEditView工IoolsIJesktcpWindowHelp口“幽|四曰|口國(guó)|國(guó)(mN)QJPW匚CTIraN30口45-r6sowgLIdBodeDiagram-2-1oi10W1010Frequencyfradfe')>>c2d(sys,0.0002,'tustin')Transferfunction:3.

2、333e-005z+3.333e-005z-0.9999Samplingtime(seconds):0.0002Y/X=3.333e-005+3.333e-005Z(-1)1-0.9999Z(-1)丫(1-0.9999Z(-1)=X(3.333e-005+3.333e-005Z(-1)Y=X*3.333e-005X+3.333e-005X(-1)+0.9999Y(-1)按照這個(gè)方程編寫(xiě)不對(duì)，因?yàn)橄禂?shù)精度太差了>>ab=tfdata(ans,'v')a=1.0e-004*0.3333222225925800.333322222592580b=1.00000000000

3、0000-0.999933335555481這才是差不多的系數(shù)"看"上門(mén)專(zhuān)門(mén)iihiCrCJiP:nrftfipsthfpaintfortriecnn'-inlr1.如府Mdu為ttU和'SIrtin*ihHNimtludB"math.h嚴(yán)dDubLeDiVransFunctdioubledlnput)<staricrioubipimpurp-staticdoublpnurput|2|-B.;IiquH0|ilTnpill；HurpuUJinpui«-(«."MU3222U)-iipw1*«Ti*Diu

4、ipun|1J;UutputipOulputiji；lnput1J-lnputtljintnin(intargcfchar*J)inti=o；FILEpFlle-MJLL；pFiZle=fopen("test-;for(i»0;i<2BBB(IO;i>><fprlntf(pFile/'Xn'CTransFuncd.0);>fclose(pFile);return0;>>step(sys)>>holdon>>step(ans)10.9StepResponse0.8sysansTime(secon

5、ds)>>plot(y)>>holdon>>plot(test(1:50001)x10/CTransFunc.cpp:Definestheentrypointfortheconsoleapplication./#include"stda僅.h"#include"stdio.h"#include"math.h"doubleCTransFunc(doubledInput)staticdoubleInput2=0.0,0.0;staticdoubleOutput2=0.0,0.0;Input0=dInput;Output0=Input0*(0.0000333322222592580)+(0.0000333322222592580)*Input1+0.999933335555481*Output1;Output1=Output0;Input1=Input0;returnOutput0;intmain(intargc,char*argv)inti=0;FILE*pFile=NULL;pFile=fopen("test.txt","w+");