[scilab-Users] Matlab to Scilab conversion

[dmcomer at dmcmicro.com]

Lester,

If you do not have the plotlib toolbox installed, I recommend it. It will handle the semilogx command
without having to convert to the Scilab plot command. Not that the plot command is bad, just easier
to convert. Post if yo uhave problems getting plot lib working. I just did this myself. Plotlib can be found
on the Scilab contributions site.

Dave
-----Original Message-----
From: Lester Anderson [mailto:lester at arctica1.wanadoo.co.uk]
Sent: Wednesday, September 2, 2009 01:30 PM
To: users at lists.scilab.org
Subject: [scilab-Users] Matlab to Scilab conversion

Hello

I am fairly new to using Scilab (and Matlab) but need to convert some existing files. I have attached an example which I can't get to work and the image of the correct plot from Matlab.

Original Matlab file:

rho_m=3300; % in kg/m^3
rho_infill=2700; % in kg/m^3
E=1e11; % in kg/m/s^2
v=0.25;
g=9.81; % in m/s^2
l1=1e3:2e3:50e3; % in m
l2=75e3:25e3:1000e3; % in m
l=[l1 l2]; % in m
k=1./l; % in 1/m

for Te=0:10e3:100e3 % in m
 D=E*Te^3/(12*(1-v^2));
 phi=1./((D.*k.^4./((rho_m-rho_infill)*g))+1);
 semilogx(k,phi)
 hold on
 title('\Phi _e(k) for various T_e')
 xlabel('Wavenumber k')
 ylabel('\Phi _e(k)')
end
text(1.7e-6,0.5,'T_e=100km')
text(3.6e-5,0.5,'T_e=10km')
text(8e-5,0.95,'T_e=0km')

Output from the conversion below:

// Display mode
mode(0);

// Display warning for floating point exception
ieee(1);

rho_m = 3300;
// in kg/m^3
rho_infill = 2700;
// in kg/m^3
E = 100000000000;
// in kg/m/s^2
v = 0.25;
g = 9.81;
// in m/s^2
l1 = 1000:2000:50000;
// in m
l2 = 75000:25000:1000000;
// in m
l = [l1,l2];
// in m
k = 1 ./l;
// in 1/m

for Te = 0:10000:100000
 // in m
 D = (E*(Te^3))/(12*(1-v^2));
 phi = 1 ./((D .*(k .^4)) ./((rho_m-rho_infill)*g)+1);
 // !! L.14: Matlab function semilogx not yet converted, original calling sequence used.
 semilogx(k,phi)
 set(gca(),"auto_clear","off")
 title("\Phi _e(k) for various T_e")
 xlabel("Wavenumber k")
 ylabel("\Phi _e(k)")
end;
// !! L.20: Matlab function text not yet converted, original calling sequence used.
text(0.0000017,0.5,"T_e=100km")
// !! L.21: Matlab function text not yet converted, original calling sequence used.
text(0.000036,0.5,"T_e=10km")
// !! L.22: Matlab function text not yet converted, original calling sequence used.
text(0.00008,0.95,"T_e=0km")

Hopefully someone can see what is needed to get things working. Thanks in advance all.

Lester Anderson

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20090902/3c163c82/attachment.htm>