Matlab vs. Autolisp – Disegno di un grafico

Pubblicato il di

I dati analitici o sperimentali devono in alcuni casi essere presentati in modo accurato e unico. A volte non è sufficiente utilizzare tutte le opzioni grafiche dei programmi di matematica o dei fogli di calcolo. Per gestire i dati di un grafico nel modo più versatile e potente, può essere utile importare i dati numerici di un grafico all’interno di un software CAD. La rappresentazione grafica deve essere preparata in anticipo prima dell’esportazione, mentre la suddivisione di ogni grafico in un layer specifico permette di effettuare qualsivoglia modifica sul CAD.

Funzioni y=f(x) importate all’interno di un CAD
Il codice riportato è un esempio di creazione di un file Autolisp. Il file consente di disegnare due curve lineari e due istogrammi. I dati appartenenti a ciascun grafico sono caricati in quattro layer distinti all’interno del CAD. Il file lisp “grafici.lsp” viene caricato in autocad effettuando un load dalla riga di comando. Tutti i comandi in autolisp vanno racchiusi tra parentesi tonde. 
1

(load "d:\grafici.lsp")
Dopo aver caricato il file in autocad bisogna far partire la funzione presente, definita in precedenza. In questo caso la funzione creata si chiama: disegno_grafici (). Il codice è scritto in linguaggio matlab. 
1

(dis_grafici)
  
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53

function scrivi_lsp (fattore,larghezza,altezza)
 
% Autori
% Prof. Danilo Pasquini IPSIA Parodi Delfino, Colleferro (Roma)
% Prof. Paolo Sordi IPSIA A. Pacinotti, Pontedera (Pisa) 
% Creazione file LISP per importazione di grafici in Autocad
%
close all
 
x = -10:0.1:10;
y = 0.15* x.^2;
plot(x,y)
y1  = 10+0.01*x.^3; 
 
fid = fopen('d:\grafici.lsp', 'wt');
fprintf(fid, ['(defun disegno_grafici () \n']);
 
fprintf(fid, ['(command "._layer" "_M" "LAYER1" "")  \n']);
fprintf(fid, ['(command "._layer" "_M" "LAYER2" "")  \n']);
fprintf(fid, ['(command "._layer" "_M" "LAYER3" "")  \n']);
fprintf(fid, ['(command "._layer" "_M" "LAYER4" "")  \n']);
 
for i=2:length(y)
%       fprintf(fid, ['(command "_color" "T" "' num2str(100) ',' num2str(100) ',' num2str(100) '" )\n']);
        fprintf(fid, ['(setvar "CLAYER" "LAYER1")\n']);
		fprintf(fid, ['(command "_line" "' num2str(x(i-1),'% 10.2f') ',' num2str(y(i-1),'% 10.2f') '" "' num2str(x(i),'% 10.2f') ',' num2str(y(i),'% 10.2f') '" "") \n']);
 
        if (rem(i,3) == 0)
            fprintf(fid, ['(setvar "CLAYER" "LAYER3")\n']);
            min_y = min(y(i-1),y(i))
            fprintf(fid, ['(command "_line" "' num2str(x(i-1),'% 10.2f') ',' num2str(min_y,'% 10.2f') '" "' num2str(x(i-1),'% 10.2f') ',' num2str(0,'% 10.2f') '" "") \n']);
            fprintf(fid, ['(command "_line" "' num2str(x(i),'% 10.2f')   ',' num2str(min_y,'% 10.2f') '" "' num2str(x(i),'% 10.2f')   ',' num2str(0,'% 10.2f') '" "") \n']);
            fprintf(fid, ['(command "_line" "' num2str(x(i),'% 10.2f')   ',' num2str(min_y,'% 10.2f') '" "' num2str(x(i-1),'% 10.2f') ',' num2str(min_y,'% 10.2f') '" "") \n']);
            fprintf(fid, ['(command "_line" "' num2str(x(i),'% 10.2f')   ',' num2str(0,'% 10.2f')     '" "' num2str(x(i-1),'% 10.2f') ',' num2str(0,'% 10.2f') '" "") \n']);
        end
%       fprintf(fid, ['(command "_color" "T" "' num2str(150) ',' num2str(150) ',' num2str(150) '" )\n']);
        fprintf(fid, ['(setvar "CLAYER" "LAYER2")\n']);
		fprintf(fid, ['(command "_line" "' num2str(x(i-1),'% 10.2f') ',' num2str(y1(i-1),'% 10.2f') '" "' num2str(x(i),'% 10.2f') ',' num2str(y1(i),'% 10.2f') '"  "") \n']);
        if (rem(i,4) == 0)
            fprintf(fid, ['(setvar "CLAYER" "LAYER4")\n']);
            min_y = max(y1(i-1),y1(i))
            fprintf(fid, ['(command "_line" "' num2str(x(i-1),'% 10.2f') ',' num2str(min_y,'% 10.2f') '" "' num2str(x(i-1),'% 10.2f') ',' num2str(20,'% 10.2f') '" "") \n']);
            fprintf(fid, ['(command "_line" "' num2str(x(i),'% 10.2f')   ',' num2str(min_y,'% 10.2f') '" "' num2str(x(i),'% 10.2f')   ',' num2str(20,'% 10.2f') '" "") \n']);
            fprintf(fid, ['(command "_line" "' num2str(x(i),'% 10.2f')   ',' num2str(min_y,'% 10.2f') '" "' num2str(x(i-1),'% 10.2f') ',' num2str(min_y,'% 10.2f') '" "") \n']);
            fprintf(fid, ['(command "_line" "' num2str(x(i),'% 10.2f')   ',' num2str(20,'% 10.2f')     '" "' num2str(x(i-1),'% 10.2f') ',' num2str(20,'% 10.2f') '" "") \n']);
        end        
end
 
fprintf(fid, [') \n']);
plot (x,y)
hold on
plot (x,y1)
fclose(fid)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22

(defun dis_grafici () 
(command "._layer" "_M" "LAYER1" "")  
(command "._layer" "_M" "LAYER2" "")  
(command "._layer" "_M" "LAYER3" "")  
(command "._layer" "_M" "LAYER4" "")  
(setvar "CLAYER" "LAYER1")
(command "_line" "-10.00,15.00" "-9.90,14.70" "") 
(setvar "CLAYER" "LAYER2")
(command "_line" "-10.00,0.00" "-9.90,0.30"  "") 
(setvar "CLAYER" "LAYER1")
(command "_line" "-9.90,14.70" "-9.80,14.41" "") 
(setvar "CLAYER" "LAYER3")
(command "_line" "-9.90,14.41" "-9.90,0.00" "") 
(setvar "CLAYER" "LAYER2")
(command "_line" "-9.90,0.30" "-9.80,0.59"  "") 
(setvar "CLAYER" "LAYER1")
(command "_line" "-9.80,14.41" "-9.70,14.11" "") 
(setvar "CLAYER" "LAYER2")
(command "_line" "-9.80,0.59" "-9.70,0.87"  "") 
(setvar "CLAYER" "LAYER2")
(command "_line" "9.90,19.70" "10.00,20.00"  "") 
)