-
Notifications
You must be signed in to change notification settings - Fork 11
/
Wav4D.m
executable file
·51 lines (44 loc) · 1.33 KB
/
Wav4D.m
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
function [z] = Wav4D(w,x,y,z,q)
% WAVELETS EN 4 DIMENSIONES
% function [z] = Wav4D(w,x,y,z,q)
% w - wavelet base
% x - rango eje x
% y - rango eje y
% z - rango eje z
% q - rango eje q
switch w
case 1
% Primer eje
z1 = (2 / (sqrt(3*sqrt(pi))) )*(1 - (x).^2);
z2 = exp( (-(x).^2) / 2);
z3 = z1.*z2;
% Segundo eje
z4 = (2 / (sqrt(3*sqrt(pi))) )*(1 - (y).^2);
z5 = exp( (-(y).^2) / 2);
z6 = z4.*z5;
% Tercer eje
z7 = (2 / (sqrt(3*sqrt(pi))) )*(1 - (z).^2);
z8 = exp( (-(z).^2) / 2);
z9 = z7.*z8;
% Cuarto eje
z10 = (2 / (sqrt(3*sqrt(pi))) )*(1 - (q).^2);
z11 = exp( (-(q).^2) / 2);
z12 = z10.*z11;
z = z3.*z6.*z9.*z12;
%z = z./4;
otherwise
z = x;
end
% Funciona
% close, clear all
% x = [-1:0.1:1];
% y = [-1:0.1:1];
% z = [-1:0.1:1];
% q = [-1:0.1:1];
% a = [1]; % a = [2....]
% bx = [0];
% by = [0];
% bz = [0];
% bq = [0];
% (2^(a(1)/2)) * Wav4D(1, ( (2^(a(1))*x) - bx(1)), ( (2^(a(1))*y) - by(1)), ( (2^(a(1))*z) - bz(1)) , ( (2^(a(1))*q) - bq(1)) )
% plot((2^(a(1)/2)) * Wav4D(1, ( (2^(a(1))*x) - bx(1)), ( (2^(a(1))*y) - by(1)), ( (2^(a(1))*z) - bz(1)) , ( (2^(a(1))*q) - bq(1)) ))