Task

Choose two methods from Multidimensional data: direct methods
Make visualizations
Prepare presentation and upload.

library(pracma)

## Warning: paketas 'pracma' buvo sukurtas pagal R versiją 4.1.3

library(ggplot2)
library(plot3D)
library(rgl)
library(viscomplexr)
library(philentropy)

## Warning: paketas 'philentropy' buvo sukurtas pagal R versiją 4.1.3

library(treemap)

## Warning: paketas 'treemap' buvo sukurtas pagal R versiją 4.1.3

library(parcoords)

## Warning: paketas 'parcoords' buvo sukurtas pagal R versiją 4.1.3

x <- seq(-5, 5, len =  30)
y <- seq(-5, 5, len = 30)
k <- meshgrid(x, y)

s <- complex(real = k$X, imaginary = k$Y)
z <- zeta(k$X + k$Y*1i)
modulus <- abs(z)

min(abs(z))

## [1] 0.001432992

max(abs(z))

## [1] 4.187367

max(modulus[,])

## [1] 4.187367

surface3d(x = k$X, y = k$Y, z = modulus, 
        xlim = c(-5, 5), 
        ylim = c(-5, 5), 
        zlim = c(0, 5), 
        col = abs(z), 
        xlab = "kva kbcs", ylab = "", zlab = "",
        main = "Zeta function values and zero")
axes3d(edges = c("x--", "y--", "z-+"))

rglwidget()

x <- seq(-5, 5, len = 30)
x <- x[! x %in% c(0.79310345,  1.03448276 , 1.27586207)]
length(x)

## [1] 30

y <- seq(-40, 40, len = 30)
k <- meshgrid(x, y)
s <- complex(real = k$X, imaginary = k$Y)
z <- zeta(k$X + k$Y*1i)
min(abs(z))

## [1] 0.01467133

modulus <- (abs(z))^-1



surface3d(x = k$X, y = k$Y, z = modulus, scale = TRUE, plot = TRUE ,
        xlim = c(-5, 5),
        ylim = c(-5, 5),
        zlim = c(0, 5),
        col= (abs(z)),
        xlab = "Real part", 
        ylab = "Imeginary part",
        zlab = "Modulus values",
        main = "Zeta function values and zero")

axes3d(edges = c("x--", "y--", "z-+"))
rglwidget()

xx <- seq(-40, 10, len = 100)
yy <- seq(-2, 48, len = 100)
zz <- zeta(xx + yy*1i)

kkk <- meshgrid(xx, yy)

z <- zeta(kkk$X +kkk$Y*1i)


phasePortrait("exp(1/z)", xlim = c(-1,1), ylim = c(-1,1), main = "Zeta function phases")

## 
## Registering parallel workers ... 3 parallel workers registered ...
## Building z plane array ...
## .making 1 blocks ... parallel loop starting ... done.
## Evaluation loop starting ... 
## .processing block 1 ... parallel loop starting ... done.
## Transforming function values into colors ...

## 
## .transforming block 1 ... parallel loop starting ... done.
## Combine color rasters ... done.
## Plotting raster image ... done.
## Deleting temporary files ... done.
## 
## Parallel backend with 3 cores remains registered for convenience.
## Can be de-registered manually with 'foreach::registerDoSEQ()'.

Data Visualization Home work 6

Antanas Kaminskas

2022-04-20

Task