Hi all
I'm working on a project which is related to programming in ipython. now i need to draw a diagram according to following equation:
f(s)=(sigma from n=1 to 18)(c(n)/(s-a(n))
how i can define such equation?
the amounts of c(n) and a(n) are complex numbers and 0
# -*- coding: utf-8 -*-
from __future__ import print_function, absolute_import
import numpy as np
length=np.arange(101)
Data = np.random.random((18,2))
c = np.empty(Data.shape[:-1], dtype=complex)
c.real = Data[...,0]
c.imag = Data[...,1]
Data = np.random.random((18,2))
a = np.empty(Data.shape[:-1], dtype=complex)
a.real = Data[...,0]
a.imag = Data[...,1]
f=np.zeros(101, dtype=complex)
for s in length:
for n in np.arange(18):
f[s] += ( c[n] / (s-a[n]) )
I think you want something like that?
thanks for your answer dear Maxime
i just need a reference to understand how to fill blanks in upper program and a little definition about these codes.
can you suggest me any reference?!
I wrote upper code with my data as follow. but it says (for a and c), "too many indices for array."
# -*- coding: utf-8 -*-
from __future__ import print_function, absolute_import
import numpy as np
length=np.arange(101)
Data = np.random.random((18,2))
c = np.empty(Data.shape[:-1], dtype=complex)
c.real = Data[-3000,-83000,-5,-5,-20,-20,6000,6000,40,40,90,90,5000,5000,1000,1000,-5000,-5000]
c.imag = Data[0,0,7000,-7000,18000,-18000,45000,-45000,60000,-60000,10000,-10000,80000,-80000,45000,-45000,92000,-92000]
Data = np.random.random((18,2))
a = np.empty(Data.shape[:-1], dtype=complex)
a.real = Data[-4500,-41000,-100,-100,-120,-120,-3000,-3000,-200,-200,-1500,-1500,-500,-500,-1000,-1000,-2000,-2000]
a.imag = Data[0,0,5000,-5000,15000,-15000,35000,-35000,45000,-45000,45000,-45000,70000,-70000,73000,-73000,9000,-9000]
f=np.zeros(101, dtype=complex)
for s in length:
for n in np.arange(18):
f[s] += ( c[n] / (s-a[n]) )
what i have to do?
I think you have to learn python a bit more before attempt to solve equation.
the issue here is that you are indexing the Data list.
I wrote Data as a random number just because I didn't know your values.
for c you just have to set the real and imaginary values. For instance
c = np.empty(18, dtype=complex)
c.real = [-3000,-83000,-5,-5,-20,-20,6000,6000,40,40,90,90,5000,5000,1000,1000,-5000,-5000]
c.imag =[0,0,7000,-7000,18000,-18000,45000,-45000,60000,-60000,10000,-10000,80000,-80000,45000,-45000,92000,-92000]
- Badges
- Science topic
- Similar topics
- Computer Science
- Program