Help with iteration problem

2 次查看(过去 30 天)
mohsen
mohsen 2012-11-19
>> x=0.7; % x is fixed value
>> syms s % s is symbols
>> v=char(x./(s^2+1)); % v is function of (x,s)
% numerical method to compute inversion Laplace transformation
% this function is available through remark at below
>> [t1,ft1]=INVLAP(v,0.0001,5,1000,6,39,89);
>> plot(t1,ft1)
please, help me the previous small code run easily, but i need (x not equal fixed), i.e., i need x- run over the interval (0:0.1:1), and also make all the following steps, then make 3-D figures between mesh(x,t1,ft1)
% remark
% INVLAP Numerical Inversion of Laplace Transforms
  1 个评论
mohsen
mohsen 2012-11-20
thank u, i do the step which u refer to it in your answer but their exist error and the following message is appeared
??? Error using ==> mesh at 80
Z must be a matrix, not a scalar or vector.
code according to your vision
for x=0:0.1:1;
syms s % s is symbols
v=char(x./(s^2+1)); % v is function of (x,s)
% numerical method to compute inversion Laplace transformation
% this function is available through remark at below
[t1,ft1]=INVLAP(v,0.0001,5,1000,6,39,89);
end
mesh(x,t1,ft1)

请先登录,再进行评论。

回答(1 个)

bym
bym 2012-11-19
seems like a simple 'for loop' would do the trick
x = 0:.1:1;
for x
v=char(x./(s^2+1)); % v is function of (x,s)
% modify call to INVLAP to store ft1 in a matrix
...
end
mesh(...)
  2 个评论
mohsen
mohsen 2012-11-20
thank u, i do the step which u refer to it in your answer but their exist error and the following message is appeared
??? Error using ==> mesh at 80
Z must be a matrix, not a scalar or vector.
code according to your vision
for x=0:0.1:1;
syms s % s is symbols
v=char(x./(s^2+1)); % v is function of (x,s)
% numerical method to compute inversion Laplace transformation
% this function is available through remark at below
[t1,ft1]=INVLAP(v,0.0001,5,1000,6,39,89);
end
mesh(x,t1,ft1)
bym
bym 2012-11-20
yes, that is why I put in my answer the comment
% modify call to INVLAP to store ft1 in a matrix

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Number Theory 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by