Help with iteration problem

조회 수: 1 (최근 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에서 Mathematics에 대해 자세히 알아보기

태그

Community Treasure Hunt

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

Start Hunting!

Translated by