clear;
Nx=35; // taille de la grille
Ny=27;
// Ny=5; // test
Nxaff = 30; // taille affichee
Nyaff = 26;
// Nyaff=4; // test
L=2; // largeur (loi empirique)
points_initiaux_x = zeros(Nx,Ny); // initialisation
points_initiaux_y = zeros(Nx,Ny);
for i=1:Nx
for j=1:Ny
points_initiaux_x(i,j) = i;
points_initiaux_y(i,j) = j;
end
end
points_finaux_x = points_initiaux_x;
// fonction
deff('[e]=profil(r)','e=exp(-(r/L)^2)'); // profil empirique
// programme principal
xmin = (Nx-Nxaff)/2+1;
xmax = Nx-xmin;
ymin = (Ny-Nyaff)/2+1;
ymax = Ny-ymin;
unite=0;dizaine=0; // pour le nom de fichier
for r=-1:0.5:Nx+1 // position du front d'onde
intensite = 1.5; // variation avec la distance
for i=1:Nx
for j=1:Ny
// e = deplacement du point
e = intensite*profil(i-r);
points_finaux_x(i,j) = points_initiaux_x(i,j) + e;
end
end
// dessin
clf;
plot2d(0,0,-1,"010"," ",[xmin,ymin,xmax,ymax])
for i=1:Nx
for j=1:Ny
xpoly(points_finaux_x(i,:),points_initiaux_y(i,:));
xpoly(points_finaux_x(:,j),points_initiaux_y(:,j));
end
end
nom='test\test'+code2str(dizaine)+code2str(unite)+'.gif';
xs2gif(0,nom,1);
unite = unite+1;
if unite == 10 then
dizaine = dizaine + 1;
unite = 0 ;
end
end
Лиценцованє
Ја, носилац ауторских права над овим делом, објављујем исто под следећим лиценцама:
Дате допущенє же би ше копировало, дистрибуовало и/або меняло тот документ под условиями ҐНУ-овей лиценци за шлєбодну документацию, верзия 1.2 або було хтора нова верзия хтору обяви Фондация за шлєбодни софтвер;без нєпременююцих оддзелькох и без текста на насловним и остатнїм боку. Текст лиценци можеце пречитац ту.http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue
шириц другим – копировац, розподзельовац и преношиц дїло
модификовац – преробйовац дїло
Под шлїдуюцима условиями:
авторство – Мушице дац одвитуюци заслуги, обезпечиц вязу ґу лиценци и назначиц чи вименки зробени. Можеце то зробиц на було хтори розумни способ, лєм нє на способ хтори предклада же лиценцор одобрує вас або вашо хаснованє.
Ова лиценца је додата на ову датотеку као део ажурирања GFDL лиценце.http://creativecommons.org/licenses/by-sa/3.0/CC BY-SA 3.0Creative Commons Attribution-Share Alike 3.0truetrue
{{Information |Description={{en|Propagation of a plane compression wave (impulse); made with [http://www.scilab.org/ Scilab] and Jasc Animation Shop 2.02}} {{fr|Propagation d'une onde de compression plane (impulsion) ; créé avec [http://www.scilab.org/