python - I have a table of values that I want to integrate with other functions -

-



i have table of lot of data points, , want use them function multiply other functions , integrate. i'm not sure how go this. figured out how integrate table using simpsons rule, , tried use numpy.polynomial.plynomial.polyfit best fit curve points values way off here have far: 

import numpy np scipy import stats, integrate import matplotlib.pyplot plt scipy.integrate import quad import math  s=5.670367e-8  a,b=np.array(np.loadtxt('cagn2.txt',  skiprows=1, unpack=true)) #change file  based on agn want. freq=np.arange(1e12,1e24,7e16)    a1=a*2.41798926215e14  def f(nu):     return nu**-1   #changed based on want  def total(nu):     return b * nu* f(nu)  """ here want place function describing data points instead of 'b'"""  e,err=quad(total,0,np.inf) #print(e) =integrate.simps(total,a1,even='first') #integrates table only, not want #print(i)  c1,c2,c3,c4,c5=np.polynomial.polynomial.polyfit(a1,b,4)   def y(x): return c1+ c2*x + c3*x**2+c4*x**3+ c5*x**4 n=y(freq)  ax=plt.gca() ax.set_xscale('log') ax.set_yscale('log') plt.plot(a1,b) plt.plot(freq,n) plt.show()

this graph got, orange polyfit. graph of polyfit

as can see way off. know how can closer fit or how use table directly 'b' in total function? thanks





wiki

Comments

Post a Comment

Popular posts from this blog

Asterisk AGI Python Script to Dialplan does not work -

-
i'm stuck project have check database of sip extension incoming caller id , dial sip extension. below receive when run script on asterisk cli. -- registered sip '9125' @ 192.168.1.102:50142 > saved useragent "zoiper rv2.8.40" peer 9125 [2017-08-22 13:53:33] notice[5095]: chan_sip.c:24558 handle_response_peerpoke: peer '9125' reachable. (17ms / 2000ms) == using sip rtp tos bits 184 == using sip rtp cos mark 5 -- executing [0112617769@from-trunk:1] answer("sip/obitrunk1-0000000d", "") in new stack > 0x1ca50d0 -- probation passed - setting rtp source address 192.168.1.98:16834 -- executing [0112617769@from-trunk:2] agi("sip/obitrunk1-0000000d", "test.py") in new stack -- launched agi script /var/lib/asterisk/agi-bin/test.py args: ['/var/lib/asterisk/agi-bin/test.py'] env line: agi_request: test.py env line: agi_channel: sip/obitrunk1-0000000d env line: agi_language: en env line: a
Read more

python - Read npy file directly from S3 StreamingBody -

-
i have many .npy (numpy) , .json files saved in s3, need download , load these files. for json files downloading files , loading dictionary using following code: obj = s3.object(bucket, key) response = obj.get() body = response['body'].read().decode('utf-8') data = json.loads(body) that way avoid writing .json file disk , reading it. now want same .npy files, equivalent? currently downloading file, saving temp file on disk , loading np.load avoid if possible. following https://stackoverflow.com/a/28196540/3509999 : obj = s3.object(bucket, key) io.bytesio(obj.get()["body"].read()) f: # rewind file f.seek(0) arr = np.load(f) wiki
Read more

kotlin - Out-projected type in generic interface prohibits the use of metod with generic parameter -

-
i start using kotlin , defined interface this: interface aadapter<vh : recyclerview.viewholder> { fun oncreateaviewholder(parent: viewgroup): vh fun onbindaviewholder(v: vh, position: int) } and when try use on code: class klasa ( private val adapter: aadapter<*> ) { fun dosth(){ //... val vh = this.adapter.oncreateaviewholder(parent) //on below line error adapter.onbindaviewholder(v, position) //... } } i error out-projected type 'aadapter<*>' prohibits use of 'public abstract fun onbindaviewholder(v: t, position: int): unit i tied add "in" or "out" definition i'm confused. how allow it. replace aadapter<*> aadapter<recyclerview.viewholder> . class klasa ( private val adapter: aadapter<recyclerview.viewholder> ) { // ... because have declared aadapter has type parameter vh type subclass of recyclerview.
Read more