Reg :: Speed of this script

[Sumit Adhikari]

Hello All,
   Following is a script in which I am suffering from speed problem. 
Also if I do little further I get stacksize problem (the problem I am 
facing from
   last several years with scilab). Is it possible to do something with 
this script ? What experts has to say ?


printf("\n ---------------------------------- \n");
printf(" Calculating Frequency Response \n");
printf(" ---------------------------------- \n\n");
clear all ;

//stacksize(5000001);
stacksize('max');

a1 = 0.5 ;
a2 = 2.0 ;

z=poly(0,'z');

Hsdm_sig   = syslin('c',(a1*a2), (z^2 + (a2 - 2)*z + (a1*a2 - a2 + 1)));
Hsdm_noise = syslin('c',(z^2 + z -2.0), (z^2 + (a2 - 2)*z + (a1*a2 - a2 
+ 1)));

Hsdm = Hsdm_sig ;


M = 1024 ;
Hn = syslin('c',(z^M - 1.0), (z^M));
Hcic_num = Hn*Hn*Hn;

Hd = syslin('c',(z- 1.0), (z));
Hcic_den = Hd*Hd*Hd ;

Hcic = Hcic_num / Hcic_den ;



H = (Hsdm * Hcic)/2^14;

clear Hsdm_sig Hsdm_noise Hsdm Hd Hcic_den Hm Hcic_num Hcic ;

Npoints = 1024*1024;
[hzm,fr]=frmag(H,Npoints);

clear H ;

N = 512;
SBW = Npoints;
frm = fr*SBW;
//clf();
plot(log10(frm(2:N)),20*log10(hzm(2:N)),"red"); mtlb_grid("on");
//plot(frm(2:N),20*log10(hzm(2:N)),"red");
//bode(H);
//plzr(H);


Best Regards,

-- 

--------------------------------------------
Sumit Adhikari
System Design Engineer
austriamicrosystems AG
Business Unit : Automotive
Mob : 00-91-9885271710/00-91-9000161710