Octave Programming Examples

pdfs

octave

function [ N_turn , l_air_gap , B_max ] = Inductor_Winding_Calculation_Function( L , parameter , A_core , l_core , I_max , B_sat , method )
%% inputs
% L is the desired inductance
% A_core is the cross-section area of the core
% l_core is the length of the magnetic loop
% I_max is the maximum current
% B_sat is the saturation flux density
% parameter and method should correspond to each other:
% if the parameter is the mu_r of the core material,
% then method should be 'mu_r';
% if the parameter is the A_L of the core material,
% then method should be 'A_L'.
%% outputs
% N_turn is the result number of turns
% l_air_gap is the result air gap length
% B_max is the result maximum flux density
%%
mu_0 = 4 * pi * 1e-7;
switch method
    
    case 'mu_r'
        
        mu_r = parameter;
        
        % Magnetic Permeability
        mu = mu_0 * mu_r;
        % Reluctance
        R_m = l_core / mu / A_core;
                
    case 'A_L'
        
        A_L = parameter;
        R_m = 1 / A_L;
        mu = l_core / R_m / A_core;
        mu_r = mu / mu_0;
        
    otherwise
        
end
% Number of turns
N_turn = sqrt(L*R_m);
        
% Check Magnetic Saturation
B_max = N_turn * I_max / R_m / A_core;
if B_max > B_sat    % if saturated, then insert air gap
    N_turn = I_max * L / B_max / A_core;
    R_m = N_turn^2 / L;
    l_air_gap = mu_0 * A_core * R_m - l_core / mu_r;
else                % if not saturated, then no need to insert air gap
    l_air_gap = 0;
end
end

top3dnew.m

octave

x = [1 2 3 4 5 6 7 8 9 10];
y1 = [.16 .08 .04 .02 .013 .007 .004 .002 .001 .0008 ];
y2 = [.16 .07 .03 .01 .008 .003 .0008 .0003 .00007 .00002 ];

semilogy(x,y1,'-bo;y1;',x,y2,'-kx;y2;');
title('Plot title');
xlabel('X Axis');
ylabel('Y Axis');

print -dpng figure.png

f=exp(-0.691213)+((cos(0.691213)).^2)-3.*0.691213+1;

octave

f=exp(-0.691213)+((cos(0.691213)).^2)-3.*0.691213+1;

perceptron

octave

NoOfStudents = 6000;
TeachingStaff = 150;
NonTeachingStaff = 20;
Total = NoOfStudents + TeachingStaff ...
   + NonTeachingStaff;
disp(Total);

representation des fonction + some damn coooooool stuff_fractale de Sierpinski

octave

x = [1 2 3 4 5 6 7 8 9 10];
y1 = [.16 .08 .04 .02 .013 .007 .004 .002 .001 .0008 ];
y2 = [.16 .07 .03 .01 .008 .003 .0008 .0003 .00007 .00002 ];

% semilogy(x,y1,'-bo;y1;',x,y2,'-kx;y2;');
% title('Plot title');
% xlabel('X Axis');
% ylabel('Y Axis');
% print -dpng figure.png
%__________________________________


% plot([1 2 3],[2 4 6]);
% print -dpng figure(1).png

% scatter((1:20),rand(1,20));
% print -dpng figure(2).png

%__________________________________

% x = 0 : 0.1: 6*pi;
% y = sin(x);

% plot(x,y) 
% xlabel('x') ;
% ylabel('sin(x)');
% title('Sin(x) graph')
% grid on, axis equal
% print -dpng figure.png

%__________________________________

% x = [0 : 0.1 : 10];
% y = sin(x);
% z = cos(x);
% plot(x, y,'-bo;y;', x, z,'rx');
% legend('Sin(x)', 'Cos(x)');
% grid on

% print -dpng figure(1).png



%__________________________________

% x = 0:pi/100:2*pi;
% y = sin(x);
% plot(x,y,'bo;y;')

% hold on

% y2 = cos(x);
% plot(x,y2,'g--')
% legend('sin','cos')

% hold off
% print -dpng figure.png

%__________________________________

% x = [0 : 0.01 : 10];
% y = exp(-x).*sin(2*x + 3);
% plot(x,y,'r.-')
% axis([0 10 -1 1])        %axis ( [xmin xmax ymin ymax] )
% print -dpng figure.png

%__________________________________

% x = [0:0.01:5];
% y = exp(-1.5*x).*sin(10*x);

% subplot(1,2,1)
% plot(x,y)
% xlabel('x')
% ylabel('exp(-1.5x)*sin(10x)')
% axis([0 5 -1 1])
% print -dpng figure.png

% subplot(1,2,2)
% plot(x,y)
% xlabel('x')
% ylabel('exp(–2x)*sin(10x)')
% axis([0 5 -1 1])
% print -dpng figure.png
%__________________________________

figure(1);hold on;axis off;clf;
Nmax=100; % Nombre d iterations
X=[];Y=[];
A=[0;0];B=[1;0];C=[0;1];
% On tire un point au hasard
x=rand(1);y=rand(1);
if (x+y>1),x=x/2;y=y/2;end
  
  for bcl=1:Nmax
        
     % On choisit le point A, B ou C au hasard
     tmp=rand(1);
     if (tmp<1/3)        % alors point A
         x=x/2;y=y/2;
     elseif (tmp>2/3)    % alors point C
         x=x/2;y=(1+y)/2;
     else                % sinon point B
         x=(1+x)/2;y=y/2;
     end
     X=[X x]
     Y=[Y y]
  end

plot(X,Y,'b.')
print -dpng figure.png
% __________________________________

Execute MATLAB/Octave Online

octave

y=-3:.1:3;
x=-2:.1:2;
[ya,xa] = ndgrid(y,x);
a2=4; b2=9;
z=9*(xa.^2/a2-ya.^2/b2);
surfl(x,y,z);    
shading interp;                       
colormap(pink)

Execute MATLAB/Octave Online

octave

x = [1 2 3 4 5 6 7 8 9 10];
y1 = [.16 .08 .04 .02 .013 .007 .004 .002 .001 .0008 ];
y2 = [.16 .07 .03 .01 .008 .003 .0008 .0003 .00007 .00002 ];

semilogy(x,y1,'-bo;y1;',x,y2,'-kx;y2;');
title('Plot title');
xlabel('X Axis');
ylabel('Y Axis');

print -dpng figure.png

Execute MATLAB/Octave Online

octave

p = [
   1.8678107616286154e+003
   8.2582338306857423e+001
  -4.1881079073449552e+000
   9.5250880790185941e-002
  -9.3549730022735955e-004
   3.3431367262149853e-006
];
polyval(p,4)

Execute MATLAB/Octave Online

octave

p = [1 7 0  -5 9];
X = [1 2 -3 4; 2 -5 6 3; 3 1 0 2; 5 -7 3 8];
polyvalm(p, X)

flower170

octave

%in the name of God
turn = randi(120)
r1 = randi(20)
r2 = randi(1200)
speed = r1 + 1 / r2

tetha = 0:0.005:2*turn*pi;
r =10;

x = cos(tetha) * r + cos(speed*tetha) * 6;
y = sin(tetha) * r + sin(speed*tetha) * 6;
plot(x, y, 'r');
ylabel(['turn = ',num2str(turn),'         speed = ', num2str(speed)])
axis equal;
axis([-16.1 16.1 -16.1 16.1])
print -dpng figure.png

Follow

Newsletter