Recent content by confused_engineer

  1. C

    Can I improve the sinc interpolation?

    Thanks. I will look into it.
  2. C

    Can I improve the sinc interpolation?

    Hello everyone. I am looking for a way to interpolate data and obtaining in return an expression which can be used to evaluate the expression at any point. An example of this would be the Shannon interpolation (see slides 15 to 18) which return a sum of cardinal sines and can be evaluated...
  3. C

    Can I improve the sinc interpolation?

    Hello everyone. I am working with mathematica, where I have developed a two-dimensional shannon interplation, just as can be seen in the slides 15 to 18 of this presentation. The code is as follows: savedX = Table[XposX = mat[[All, 1]]; YposX = mat[[All, 2]]; windXVal = mat[[All, i]]...
  4. C

    I Implementation of Correlated Gaussian Random Fields Model

    Hello everyone. I have been recently working in an optimization model in the presence of uncertainty. I have read https://www.researchgate.net/publication/310742108_Efficient_Simulation_of_Stationary_Multivariate_Gaussian_Random_Fields_with_Given_Cross-Covariance in which, a methodology for...
  5. C

    MATLAB Problem with Karhunen-Loève expansion in matlab

    I have copied the code of the accepted answer to this post in the official Matlab forums, since I am interested in performing the KL expansion myself. clc clear all y=[1,2,4;2,3,10]; y=y' %Reasons for transposing will become clear when you will read the second point given below...
  6. C

    KarhunenLoeveDecomposition function in Mathematica

    Hello everyone. I have a vector, stochasticData.mat, it contains a matrix of size 211302*50, being 211302 measurements of 50 realizations of a stochsatic process. I want to use the Karhunen-Loève expansion and the software Mathematica to calculate the uncorrelated random variables. For that, I...
  7. C

    I Question about weights using Chebyshev polynomials as quadrature

    First of all, thanks for answering my question. I am confused because the article tallks about ωs and ωN-s and I find extrange that following the paper one arrives to ω3=ω3=4/N*... for this particular example. Therefore, I understand from your answer that ω3=ω3=4/N*... is the weight of the...
  8. C

    I Question about weights using Chebyshev polynomials as quadrature

    Hello everyone. I am studying this article since I am interested in optimization. The article makes use of Clenshaw–Curtis quadrature scheme to discretize the integral part of the cost function to a finite sum using Chebyshev polynomials. The article differentiates between the case of odd...
  9. C

    I Question about the roots of Chebyshev polynomials

    First of all thanks for your answer. If I understand propperly, the node points are the roots of the polynomial of order N? but then, what is the use of the cosine? I have plotted the values if the cosine for N=60 and these values cluster arround the endpoints of the interval [-1, 1], I think...
  10. C

    I Question about the roots of Chebyshev polynomials

    Hello everyone. I am trying to construct an optimization problem using Chebyshev pseudospectral method as described in this article. For that, I need to calculate the zeros of the Chebyshev polynomial of any order. In the article is sugested to do it as tk=cos(πk/N) k=0, ..., N...
  11. C

    Python How can I evaluate a Chebishev polynomial in python?

    Thanks for the answer, but I am afraid that is exactly the problem. I need a function that returns the evaluation of a polynomial of order n evaluated at x. For that, I need to calculate the polynomial of order n based on the polynomials of order (n-1) and (n-2). Currently, I am using...
  12. C

    Python How can I evaluate a Chebishev polynomial in python?

    Hello everyone. I need to construct in python a function which returns the evaluation of a Chebishev polynomial of order k evaluated in x. I have tested the function chebval form these documents, but it doesn't provide what I look for, since I have tested the third one, 4t^3-3t and import numpy...
Back
Top