problem with 'solve' function in matlab 2008 and upper

I have a function in matlab, I want to find its maximum.. ( my function is guassian signal, it has a 2 pulse, I want to find extermum point of each pulse and calculate distance between their extermum).

base on this data, I used this commands:

   syms t
      Z=[0.000001116863592   0.000004972731239  -0.547213084010968   1.015860491921958]
   alpha = 1.749262265091693e+013
   Fc = 4.334354827118762e+006
   phase2 = 0.991124769614895
   phase1 = -0.549482979269295
solve(diff((Z(4).*exp(-alpha.*(t-Z(2)).^2).*cos(2*pi.*Fc.*(t-Z(2))+phase2))))-solve(diff((Z(3).*exp(-alpha.*(t-Z(1)).^2).*cos(2*pi.*Fc.*(t-Z(1))+phase1))));

this commands work in Matlab 2007a, but dont work in upper version. What should I do?

[Walter Roberson]

Please clarify what you mean by "don't work".

The values you give for z and alpha and so on: are those representing outputs or commands? If they represent commands, you should probably be using subs() instead.

[fsalehi]

all of them are commands,

I meant if in Matlab 2007 the answer of these commands is 1.898976e-6, in upper version the answer is [matrix 1x1] ! I dont know why!?

would please tell me how should I use subs()?

[Walter Roberson]

You do not tell solve() what to solve for, and it is not going to be able to solve what you give if you give it in that order. You have to move the solve() command after the assignments. You should probably also use

tsol = solve(diff((Z(4).*exp(-alpha.*(t-Z(2)).^2).*cos(2*pi.*Fc.*(t-Z(2))+phase2))))-solve(diff((Z(3).*exp(-alpha.*(t-Z(1)).^2).*cos(2*pi.*Fc.*(t-Z(1))+phase1))), t);

to be sure that solve will solve for the variable t and return just that variable.

I am not certain at the moment why solve() is returning a 1x1 matrix. Try asking what the class() of the result is, and try asking what size() of it is. It might be a symbolic array of answers for which it happened to find only a single answer.

[fsalehi]

would you please see my answer?

Answer 1

[Kai Gehrs]

Hi,

I agree with Walter's recommendation. I just would like to add that the results returned by solve are SYM objects and scalar SYM objects (including numbers) are displayed as 1x1 SYM objects in the workspace/variable browser. Maybe this is the root cause for your question.

I run your example in R2011a and without checking the details it seems to work fine.

Best regarrds,

-- Kai

Answer 2

[fsalehi]

I try Walter's recommendation but it didnt work in Matlan 2009 and 2010a. I'm not sure about 2011a.

I should note that, If I change values of Z, my solve commands works in all version of matlab..

I used these codes in my gui program and it should work in all version of matlab (at least 2009 and upper), but these codes didnt answer in 2009 and 2010a.