I'm trying to figure out a way to algorithmically get the amplitude and phase of a function that has sinusoidal terms in the Maxima computer algebra system. This only applies to steady state (as t -> infinity and the transients decay). For example, a trivial case would be: f(t) = 1 / w * sin(w * t + theta) + exp(-a * t) + 8 In this ...
I have an image in MATLAB: y = rgb2gray(imread('some_image_file.jpg')); and I want to do some processing on it: pic = some_processing(y); and find the local maxima of the output. That is, all the points in y that are greater than all of their neighbors. I can't seem to find a MATLAB function to do that nicely. The best I can come ...
How to find solution to overdetermined systems in Macsyma, Scilab, Octave? ...
I have run across several posts and articles that suggests using things like simulated annealing to avoid the local minima/maxima problem. I don't understand why this would be necessary if you started out with a sufficiently large random population. Is it just another check to insure that the initial population was, in fact, sufficient...
First, I've tried the stuff you get using "load(functs)", no thank you. I have some pretty complex equations that I'd like to evaluate in Maxima. But at one step in the process I split a register in it's most-significant and least-significant components, and I need to be able to do that in Maxima -- if I am to evaluate the class of te...
Hello, I'm trying to use linsolve( ) from within a subroutine in order to code up a solution algorithm that is independent of the dimension of the problem. The desire is to dynamically allocate the equation and variable arrays within the subroutine itself, pass these to linsolve( ) and hopefully get back the results. /* SUBROUTINE IDE...
I've got a specific question and a related more general one... Why does imextendedmax() not give for example 9 in A(3,3) as a max? Generally... what is the best way for finding multiple maxes/peaks? The nice thing about imextended max is it allows a threshold where presumably everything under that threshold is not counted whereas imregio...
I'm looking for a possibility to use different 'higher' math operations in combination with R. A link or integration between R and a CAS would be the perfect solution. Which integration of R and other (math & statistic related) systems or vice verse are out there? How well do they work? What would you suggest? How expansive (in time, mo...
I have to solve the following convolution related numerical integration problem in R or perhaps computer algebra system like Maxima. Integral[({k(y)-l(y)}^2)dy] where k(.) is the pdf of a standard normal distribution l(y)=integral[k(z)*k(z+y)dz] (standard convolution) z and y are scalars The domain of y is -inf to +inf. The integral in ...
Hi All, I'm still a maxima newbie so bear with me. I am trying to write my own formula for calculating the adjoint of a matrix (I know maxima already has one built-in, but I was trying my own as a learning exercise). So far I have (for a 3x3 matrix): /* cofactor of some submatrix of the matrix, by deleting row i and column j */ cof(i, ...