Autor Tema: Graficando en Scilab la Molificación de la función valor absoluto.

0 Usuarios y 1 Visitante están viendo este tema.

11 Febrero, 2023, 02:44 am
Leído 551 veces

hikmath

  • $$\Large \color{#6a84c0}\pi$$
  • Mensajes: 24
  • País: cl
  • Karma: +0/-0
Hola. Estoy tratando de graficar las diferentes aproximaciones mediante convoluciones con la función estandar de Mollifier para la función \( f(x)=|x| \), sin embargo, mi código no logra funcionar. ¿Qué está mal? (No sé mucho de Scilab)

El código es:

Código: [Seleccionar]
x = -1:0.1:1;
f = (x > -1 & x < 1).*abs(x)+(x <= -1 | x >= 1).*0;
phi = (x > -1 & x < 1).*(exp(1/(x.^2-1))/sqrt(2*%pi))+(x <= -1 | x >= 1).*0;
y = conv(f, phi, 'same') * 0.1;
plot(x, f, 'b', x, y, 'r');
legend('f(x)', 'Convolution f * phi');
el cual, en la página de Scilab online, me entrega el siguiente error:
Código: [Seleccionar]
!--error 9999
Inconsistent element-wise operationat line       4 of function %b_x_s called by : 
phi = (x > -1 & x < 1).*(exp(1/(x.^2-1))/sqrt(2*%pi))+(x <= -1 | x >= 1).*0;
at line      10 of exec file called by :   
835eDLO7bMHkoFoEOTsnsEQVr4VFb8.sci", 2)

Las funciones a las cuales me refiero son \( f(x)=|x| \) para \( |x|<1 \) y \( f(x)=0 \) en otro caso, y la función \( \varphi_{\epsilon}=\frac{1}{\epsilon}\varphi(x/\epsilon) \) donde \( \varphi(x)=\frac{1}{\sqrt{2\pi}}\exp(1/(x^2-1)) \) para \( |x|<1 \) y 0 en otro caso.

pd: Mi objetivo es poder obtener gráficos como la siguiente imagen: