Interpolation in matlab

  • Thread starter sara_87
  • Start date
  • #1
763
0

Homework Statement



If i have an interpolation formula, say:

[tex]y(t_i)=3+h\sum^i_{j=1}jy(t_j)[/tex]

where i=0 to n
h=(a+b)/n
ti=i*h
and
when i=0, we have: y(t0)=0.

How would i write a code to find y(ti)?

Homework Equations





The Attempt at a Solution



function y=examp(a,b,n)
h=(a+b)/n
for i=1:n+1
t(i)=(i-1)*h
for j=1:i+1

and now im stuck because i dont know how write a code that uses the previous value of y to compute the next value of y.
 

Answers and Replies

  • #2
Before the first loop, I would set up the y's(actually, it would be faster for your code to init the t's to zeros to, but unless n is very large it shouldn't matter):

y = zeros(n+1,1)

Then in the loop it would be something like
y(j) = 3 + sum( y(1:j-1).*(1:j-1) )

the .* does an element-wise product
 

Related Threads on Interpolation in matlab

Replies
0
Views
1K
Replies
13
Views
2K
  • Last Post
Replies
1
Views
1K
  • Last Post
Replies
1
Views
28K
  • Last Post
Replies
1
Views
10K
  • Last Post
Replies
5
Views
3K
Replies
0
Views
1K
  • Last Post
Replies
2
Views
954
  • Last Post
Replies
1
Views
705
Top