คำถามติดแท็ก matlab

ฉันกำลังเปรียบเทียบรหัสของฉันบางส่วนกับรหัส "หุ้น" MATLAB ฉันประหลาดใจที่ผลลัพธ์ ฉันรันโค้ดตัวอย่าง (Sparse Matrix) n = 5000; a = diag(rand(n,1)); b = rand(n,1); disp('For a\b'); tic;a\b;toc; disp('For LU'); tic;LULU;toc; disp('For Conj Grad'); tic;conjgrad(a,b,1e-8);toc; disp('Inv(A)*B'); tic;inv(a)*b;toc; ผล : For a\b Elapsed time is 0.052838 seconds. For LU Elapsed time is 7.441331 seconds. For Conj Grad Elapsed time is …

36 linear-algebra  performance  matlab 

ฉันจำเป็นต้องประเมินตัวเลขอินทิกรัลด้านล่าง: ∫∞0sinc′(xr)rE(r)−−−−√dr∫0∞sinc′(xr)rE(r)dr\int_0^\infty \mathrm{sinc}'(xr) r \sqrt{E(r)} dr ที่,x∈R+และλ,κ,ν>0 นี่Kคือฟังก์ชัน Bessel ที่แก้ไขของชนิดที่สอง ในกรณีของฉันโดยเฉพาะฉันมีλ=0.00313,κ=0.00825และν=0.33E(r)=r4(λκ2+r2−−−−−−√)−ν−5/2K−ν−5/2(λκ2+r2−−−−−−√)E(r)=r4(λκ2+r2)−ν−5/2K−ν−5/2(λκ2+r2)E(r) = r^4 (\lambda\sqrt{\kappa^2+r^2})^{-\nu-5/2} K_{-\nu-5/2}(\lambda\sqrt{\kappa^2+r^2})x∈R+x∈R+x \in \mathbb{R}_+λ,κ,ν>0λ,κ,ν>0\lambda, \kappa, \nu >0KKKλ=0.00313λ=0.00313\lambda = 0.00313κ=0.00825κ=0.00825\kappa = 0.00825ν=0.33ν=0.33\nu = 0.33 0.33 ฉันใช้ MATLAB และผมก็มีความพยายามในตัวฟังก์ชั่นintegralและquadgkซึ่งจะช่วยให้ฉันมากของข้อผิดพลาด (ดูด้านล่าง) ฉันได้พยายามธรรมชาติสิ่งอื่น ๆ อีกมากมายเช่นกันเช่นการบูรณาการโดยส่วนและข้อสรุปปริพันธ์จากไป( k + 1 ) x πkxπkxπkx\pi(k+1)xπ(k+1)xπ(k+1)x\pi ดังนั้นคุณมีข้อเสนอแนะเกี่ยวกับวิธีที่ฉันควรลองต่อไปหรือไม่ UPDATE (คำถามเพิ่มเติม) ฉันอ่านกระดาษ @Pedro ที่เชื่อมโยงกับและฉันไม่คิดว่ามันยากเกินไปที่จะเข้าใจ อย่างไรก็ตามฉันมีคำถามสองสามข้อ: มันจะโอเคไหมที่จะใช้เป็นองค์ประกอบพื้นฐานψ kในวิธีของ univariate …

25 matlab  quadrature  special-functions 

คำถามนี้ถูกโยกย้ายจากการตรวจสอบข้ามเพราะมันสามารถตอบได้ในการแลกเปลี่ยนวิทยาศาสตร์กองซ้อน อพยพ 7 ปีที่ผ่านมา ฉันต้องการเรียนรู้การเขียนโปรแกรม MATLAB เพื่อให้ฉันสามารถทำการวิจัย / วิเคราะห์ได้ด้วยตัวเองและเพื่อให้ฉันสามารถศึกษา / แก้ไขสคริปต์ MATLAB ที่ฉันพบออนไลน์เป็นต้น อย่างไรก็ตามปัญหาคือฉันไม่สามารถจ่าย MATLAB ได้ GNU Octave จากสิ่งที่ฉันได้ยินมานั้นค่อนข้างเข้ากันได้กับ MATLAB ความท้าทายของ Octave คือว่าเอกสารนั้นกระจัดกระจายมาก ดังนั้น MATLAB ในด้านหนึ่งจึงมีราคาแพง (-ve) แต่มีเอกสารประกอบแบบฝึกหัดและสคริปต์ที่ดีจำนวนมากทางออนไลน์ (+) ในขณะที่ Octave ไม่มีค่าใช้จ่าย (เช่นเดียวกับเบียร์) แม้ว่าจะไม่มีเอกสารที่ครอบคลุม ฯลฯ . ฉันกำลังคิดที่จะดาวน์โหลดและติดตั้ง Octave แต่ใช้เอกสาร MATLAB (และสคริปต์) เพื่อช่วยให้ฉันดำเนินการต่อ ฉันทราบว่ามีความแตกต่างบางอย่างระหว่างอ็อกเทฟและ MATLAB - สิ่งที่ฉันพยายามทำคือความแตกต่างนั้นมีขนาดใหญ่พอที่จะขัดขวางความพยายามของฉันในการ "เรียนรู้ MATLAB โดยใช้อ็อกเทฟ" ได้อย่างมีประสิทธิภาพ …

23 software  matlab  languages  octave 

ฉันอยากรู้ว่ามีวิธีที่รวดเร็วในการคำนวณระยะทางแบบยุคลิดของเวกเตอร์สองตัวใน Octave หรือไม่ ดูเหมือนว่าไม่มีฟังก์ชั่นพิเศษสำหรับสิ่งนั้นดังนั้นฉันควรใช้สูตรด้วยsqrtหรือไม่

18 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

ฉันต้องการจำลองพฤติกรรมของระบบลูกตุ้มแบบสองเท่า ระบบนี้เป็นหุ่นยนต์หุ่นยนต์ 2 องศาอิสระที่ไม่ได้ถูกกระตุ้นและจะทำตัวเป็นเหมือนลูกตุ้มสองเท่าที่ได้รับผลกระทบจากแรงโน้มถ่วง ข้อแตกต่างที่สำคัญเพียงอย่างเดียวกับลูกตุ้มคู่คือประกอบด้วยวัตถุแข็งสองตัวที่มีคุณสมบัติมวลและความเฉื่อยที่จุดศูนย์กลางมวล โดยทั่วไปฉันตั้งโปรแกรมode45ภายใต้ Matlab เพื่อแก้ปัญหาระบบ ODE ของประเภทต่อไปนี้: ⎡⎣⎢⎢⎢10000M110M1200100M120M22⎤⎦⎥⎥⎥⎡⎣⎢⎢⎢x˙1x˙2x˙3x˙4⎤⎦⎥⎥⎥=⎡⎣⎢⎢⎢x2−V1−G1x4−V2−G2⎤⎦⎥⎥⎥[10000M110M1200100M120M22][x˙1x˙2x˙3x˙4]=[x2−V1−G1x4−V2−G2] \left[ \begin{array}{cccc} 1 & 0 & 0 & 0\\ 0 & M_{11} & 0 & M_{12}\\ 0 & 0 & 1 & 0\\ 0 & M_{12} & 0 & M_{22} \end{array} \right] \left[ \begin{array}{c} \dot{x}_1\\ \dot{x}_2\\ \dot{x}_3\\ \dot{x}_4 \end{array} …

15 matlab  ode  numerical-modelling 

คนวิทยาศาสตร์การคำนวณ: ฉันโพสต์คำถามนี้ที่ Math Stack Exchangeและมีคนแสดงความคิดเห็นว่าฉันอาจได้รับคำตอบ "ดีกว่า" ที่นี่: ฉันเป็นสามเณรด้วยวิธีการเชิงตัวเลขและ Matlab ฉันพยายามประเมินผลรวมของสามอินทิกรัลสองอันดังต่อไปนี้ (สามารถเขียนได้ง่ายกว่า แต่คุณยังไม่สามารถประเมินได้ว่าเป็นสัญลักษณ์ (?)) ฉันมีปัญหาในการทำให้ทำงานที่นี่ดังนั้นฉันลังเลที่จะแยกมันออกเป็นชิ้น ๆ ที่นี่: ฉันต้องการค้นหาผลรวมของLATEXLATEX\LaTeX 2((1/0.3)−1)2(∫1/0.31∫r11∫r1−r00F1(r0,r1,t)exp(−(0.3)2t24)dtdr0dr1),2((1/0.3)−1)2(∫11/0.3∫1r1∫0r1−r0F1(r0,r1,t)exp⁡(−(0.3)2t24)dtdr0dr1),\frac{2}{((1/0.3) - 1)^2}\left(\int_1^{1/0.3}\int_1^{r_1}\int_0^{r_1-r_0}F_1(r_0,r_1,t)\exp(-\frac{(0.3)^2 t^2}{4})\,dt\,dr_0\,dr_1 \right), และ 2((1/0.3)−1)2(∫1/0.31∫r11∫r1+r0r1−r0F2(r0,r1,t)exp(−(0.3)2t24)dtdr0dr1),2((1/0.3)−1)2(∫11/0.3∫1r1∫r1−r0r1+r0F2(r0,r1,t)exp⁡(−(0.3)2t24)dtdr0dr1),\frac{2}{((1/0.3) - 1)^2}\left(\int_1^{1/0.3}\int_1^{r_1}\int_{r_1-r_0}^{r_1+r_0} F_2(r_0,r_1,t)\exp(-\frac{(0.3)^2 t^2}{4})\,dt\,dr_0\,dr_1 \right), ที่ไหน F1(r0,r1,t)=t2r30∗(0.3)32r31π−−√F1(r0,r1,t)=t2r03∗(0.3)32r13πF_1(r_0,r_1,t)=\frac{t^2 r_0^3*(0.3)^3}{2r_1^3\sqrt{\pi}} และ F2(r0,r1,t)=(0.3)3π3/2(r0+r1−t)4(t2+2t(r0+r1)−3(r1−r0)2)2288(43πr30)(43πr31).F2(r0,r1,t)=(0.3)3π3/2(r0+r1−t)4(t2+2t(r0+r1)−3(r1−r0)2)2288(43πr03)(43πr13).F_2(r_0,r_1,t)=\frac{(0.3)^3\pi^{3/2}(r_0+r_1-t)^4 (t^2+2t(r_0+r_1)-3(r_1-r_0)^2)^2}{288(\frac{4}{3}\pi r_0^3)(\frac{4}{3}\pi r_1^3)}. แก้ไข (2 มีนาคม 2013): มีคนตอบว่าได้รับ Mathematica เพื่อทำอินทิกรัลเป็นสัญลักษณ์ ฉันพยายามทำสิ่งนี้ (ด้วยอินทิกรัลเวอร์ชันที่ง่าย) และ Mathematica สามารถทำได้เพียงสองส่วนนอกของอันแรกและหยุดชะงักในอันที่สอง ฉันขอขอบคุณความช่วยเหลือ …

15 matlab 

ฉันกำลังรู้สึกหงุดหงิดกับวิธีที่ matlab จัดการกับการรวมเชิงตัวเลขกับ Scipy ฉันสังเกตเห็นความแตกต่างต่อไปนี้ในรหัสทดสอบด้านล่าง: รุ่นของ Matlab ทำงานเร็วกว่างูหลามของฉันโดยเฉลี่ย24 เท่า ! เวอร์ชันของ Matlab สามารถคำนวณอินทิกรัลโดยไม่มีการเตือนขณะที่ python คืนค่า nan+nanj ฉันจะทำอะไรได้บ้างเพื่อให้แน่ใจว่าฉันจะได้รับประสิทธิภาพการทำงานเหมือนกันในงูใหญ่เทียบกับสองประเด็นที่กล่าวถึง? ตามเอกสารทั้งสองวิธีควรใช้ "การปรับพื้นที่สี่เหลี่ยมจัตุรัสส่วนกลาง" เพื่อประมาณค่าอินทิกรัล ด้านล่างเป็นรหัสในสองเวอร์ชัน (ค่อนข้างคล้ายกันแม้ว่างูใหญ่ต้องการให้มีการสร้างฟังก์ชั่นครบถ้วนเพื่อให้สามารถจัดการกับการรวมที่ซับซ้อนได้) หลาม import numpy as np from scipy import integrate import time def integral(integrand, a, b, arg): def real_func(x,arg): return np.real(integrand(x,arg)) def imag_func(x,arg): return np.imag(integrand(x,arg)) real_integral = integrate.quad(real_func, a, …

13 matlab  python  quadrature  numba 

ฉันได้อ่านแล้วว่าถ้าเช่นฉันมี double forloop ที่วิ่งอยู่เหนือ index ของเมทริกซ์จากนั้นการวางดัชนีการรันคอลัมน์ในลูปด้านนอกนั้นมีประสิทธิภาพมากกว่า ตัวอย่างเช่น: a=zeros(1000); for j=1:1000 for i=1:1000 a(i,j)=1; end end วิธีที่มีประสิทธิภาพมากที่สุดในการเขียนโค้ดคืออะไรถ้าฉันมีforลูปสามลูปขึ้นไป ตัวอย่างเช่น: a=zeros(100,100,100); for j=1:100 for i=1:100 for k=1:100 a(i,j,k)=1; end end end

12 matlab  efficiency 

เฮ้มีและขอบคุณที่ให้เวลาเพื่อดูคำถามของฉัน นี่เป็นคำถามที่ฉันได้โพสต์ก่อนหน้านี้ใน physics.stackexchange.com ขณะนี้ฉันกำลังศึกษาสปินทริปปินเตอร์ 2 มิติ Bose-Einstein Condensate และอยากรู้เกี่ยวกับสถานะพื้นฐานของระบบนี้ วิธีการทางคณิตศาสตร์ของการเดินทางไปสภาพพื้นดินที่เรียกว่าวิธีการที่เวลาจินตนาการ วิธีนี้ง่ายมากที่เวลาในกลศาสตร์ควอนตัมถูกแทนที่ด้วยจินตภาพหนึ่ง ทดแทนนี้ทำให้อนุภาคพลังงานสูงในระบบของฉันสลายตัวเร็วกว่าพลังงานต่ำ การทำให้จำนวนอนุภาคกลับมาเป็นปกติในทุกขั้นตอนของการคำนวณเราจะจบลงด้วยระบบของอนุภาคพลังงานต่ำสุดหรือที่รู้จัก สภาพพื้นดินt = - i τt=−iτ t = -i \tau สมการ (s) ในคำถามเป็นเชิงเรียกว่าไม่เชิงเส้นSchrödingerสมบางครั้งสม Gross-Pitaevskii เพื่อแก้ปัญหาที่ฉันใช้ Matlabs ode45 ซึ่งวิวัฒนาการระบบไปข้างหน้าในเวลาและในที่สุดก็ถึงสภาพพื้นดิน บันทึก! สมการไม่เชิงเส้นSchrödingerรวมถึง laplacian และคำที่แตกต่างอื่น ๆ ในอวกาศ ทั้งหมดนี้ได้รับการแก้ไขโดยใช้การแปลงฟูริเยร์ที่รวดเร็ว ในที่สุดเรามีเวลาเพียง ODE * * * * ปัญหาและคำถามของฉัน:การคำนวณไปจากเพื่อt_fode45 ถูกใส่ไว้ในfor for loop ดังนั้นมันจึงไม่คำนวณเวกเตอร์ยักษ์ในเวลาเดียวกัน รอบแรกจะเริ่มต้นด้วย ode45 …

12 matlab  ode  quantum-mechanics 

ฉันกำลังใช้ MATLAB เพื่อแก้ปัญหาที่เกี่ยวข้องกับการแก้A x = bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}ทุก ๆ เวลาที่ขb\mathbf{b}เปลี่ยนแปลงไปตามกาลเวลา ตอนนี้ฉันกำลังทำสิ่งนี้โดยใช้ MATLAB mldivide: x = A\b ฉันมีความยืดหยุ่นในการทำให้เป็น precomputations มากที่สุดเท่าที่จำเป็นเพื่อให้ผมสงสัยว่าถ้ามีวิธีที่เร็วขึ้นและ / mldivideหรือความถูกต้องมากกว่า โดยทั่วไปแล้วจะทำอะไรที่นี่ ขอบคุณทุกคน!

12 linear-algebra  matlab  linear-solver  matrix  linear-system 

ฉันกำลังพยายามคำนวณการจัดอันดับที่ดีสำหรับเมทริกซ์อย่างถูก ดังนั้นฉันจึงคำนวณการถอดรหัสคอลัมน์โดยใช้การถอดรหัส QRAAA [Q,R,E]=qr(A) ใน Matlab ฉันประเมินระดับAAAใช้ tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) มันใช้งานได้ดีและพล็อตในแนวทแยงทั้งหมดของ R มีลักษณะดังนี้: หากพอร์ตอัลกอริธึมทั้งหมดเป็น C / Fortran ฉันจะแทนที่ [Q, R, E] = qr (A) โดยใช้ DGEQP3 จาก LAPACK ซึ่งคำนวณคอลัมน์ที่มีการแยกย่อย QR แต่ถ้าฉันใช้การประมาณแบบเดียวกันกับอันดับฉันส่วนใหญ่เข้าใจผิด พล็อตเดียวกันสำหรับผลิตจาก DGEQP3 ดูเหมือนว่า RRR อินพุตเมทริกซ์เหมือนกันทั้งสองการทดลอง คำถามของฉันตอนนี้ฟังก์ชั่น LAPACK คือคอลัมน์ใดที่หมุนแกน QR การสลายตัวจาก Matlab พึ่งพา? ขอบคุณสำหรับความช่วยเหลือใด ๆ Grisu แก้ไข: …

12 matlab  lapack 

ฉันกำลังเรียนการคำนวณทางวิทยาศาสตร์และเราเพิ่งไปประมาณกำลังสองน้อยที่สุด คำถามของฉันเกี่ยวกับการประมาณโดยใช้พหุนาม ฉันเข้าใจว่าถ้าคุณมีจุดข้อมูล n + 1 คุณสามารถค้นหาพหุนามเฉพาะขององศา n ที่อธิบายจุดเหล่านี้ทั้งหมด แต่ฉันก็สามารถเห็นได้ว่าทำไมสิ่งนี้จึงไม่เหมาะเสมอไป คุณสามารถได้ยินเสียงรบกวนมากมายระหว่างจุดข้อมูลด้วยวิธีการดังกล่าว ฉันคิดว่าเป็นเรื่องดีที่จะได้พหุนามระดับต่ำกว่าซึ่งประมาณว่าข้อมูลของคุณดีพอ คำถามของฉันคือคุณจะตัดสินใจอย่างไรในการใช้พหุนามในระดับใด มีกฎง่ายๆหรือขึ้นอยู่กับปัญหาที่เกิดขึ้นเพียงลำพัง? เราต้องคำนึงถึงการแลกเปลี่ยนที่หลากหลายเมื่อตัดสินใจระหว่างองศาที่มากหรือน้อย? หรือฉันเข้าใจผิดบางอย่างที่นี่? ขอบคุณล่วงหน้า.

11 matlab  least-squares 

ตอนนี้ฉันกำลังผลิตภาพ 1d และ 2d โดยใช้ Matlab เพื่อเปรียบเทียบความแม่นยำกับรุ่นที่กำหนด ฉันต้องเปรียบเทียบวิธีการของฉันกับแบบจำลอง. wfn แบบ Gaussian มาตรฐานและฉันจะทำเช่นนั้นโดยการวิเคราะห์ความหนาแน่นของโมเลกุลเช่นเดียวกับ laplacian ฉันสนใจที่จะสร้างภาพสองมิติที่มีความแตกต่างในการประมาณของฉันไปยังการประมาณ. wfn (อาจผ่านทางแผนที่เส้นโค้ง) และฉันก็สนใจในคุณสมบัติ (ความหนาแน่น, laplacian ฯลฯ ) ตามเส้นทางพันธบัตรและเส้นตรงระหว่างสองอะตอม ตามที่ฉันได้กล่าวไปก่อนหน้านี้ฉันกำลังใช้ Matlab เพื่อสร้างภาพของฉันแม้ว่าภาพเหล่านี้ส่วนใหญ่จะเป็นของฉันและผู้ร่วมงานของฉัน พวกเขาดูเหมือนจะไม่เป็นประเภทของภาพที่ฉันมักจะเห็นในเอกสารและสิ่งพิมพ์อื่น ๆ ซอฟต์แวร์ประเภทใดที่ฉันควรใช้ (หรือผู้อื่นใช้) เพื่อสร้างกราฟและภาพสิ่งพิมพ์คุณภาพ?

11 matlab  computational-chemistry  publications 

ได้รับการตั้งค่าโดยพลการของ (ตัวเลข) การฝึกอบรมที่ซับซ้อนตาราง= { 1 , 2 , ⋯ , เมตร } , ฉันสนใจในการคำนวณพีชคณิตจริงเมทริกซ์ที่สร้างขึ้นโดยเรียกว่าL นั่นคือฉันต้องการพื้นฐานสำหรับ L = s พีn R { B : B ∈ ∪ ∞ k = 1 C k } ที่C kถูกกำหนดให้เป็นซ้ำA={A1,A2,⋯,Am}A={A1,A2,⋯,Am}\mathcal{A}=\{A_1,A_2,\cdots,A_m\}AA\mathcal{A}LALA\mathcal{L_\mathcal{A}}LA=spanR{B:B∈∪∞k=1Ck}LA=spanR{B:B∈∪k=1∞Ck} \mathcal{L_\mathcal{A}} = \mathbb{span_R}\{B:B\in\cup_{k=1}^{\infty}\mathcal{C}_k\} CkCk\mathcal{C}_k C k + 1 ={[X,Y]:X,Y∈ ∪ k J = 1 C …

11 linear-algebra  matlab  basis-set 

ฉันกำลังใช้กระดาษ " การขนส่งมวลชนที่เหมาะสมที่สุดสำหรับการลงทะเบียนและแปรปรวน " เป้าหมายของฉันคือการทำให้มันออนไลน์เพราะฉันไม่สามารถหารหัสขนส่งมวลชนแบบ eulerian ออนไลน์ได้และสิ่งนี้น่าสนใจอย่างน้อยสำหรับชุมชนการวิจัยในการประมวลผลภาพ กระดาษที่สามารถสรุปได้ดังนี้: - หาแผนที่เริ่มต้นuuuใช้ 1D จ้อ histogram พร้อม x และ y พิกัด - แก้ปัญหาสำหรับจุดคงที่ของut=1μ0Du∇⊥△−1div(u⊥)ut=1μ0Du∇⊥△−1div(u⊥)u_t = \frac{1}{\mu_0} Du \nabla^\perp\triangle^{-1}div(u^\perp)u⊥u⊥u^\perp△−1△−1\triangle^{-1}DuDuDudt&lt;min|1μ0∇⊥△−1div(u⊥)|dt&lt;min|1μ0∇⊥△−1div(u⊥)|dt<\min|\frac{1}{\mu_0}\nabla^\perp\triangle^{-1}div(u^\perp)| สำหรับการจำลองเชิงตัวเลข (ดำเนินการในตารางปกติ) พวกเขาระบุว่าใช้poicalcของmatlabสำหรับการแก้สมการปัวซองพวกเขาใช้ความแตกต่างที่มีศูนย์กลาง จำกัด สำหรับอนุพันธ์เชิงพื้นที่ยกเว้นDuDuDuซึ่งคำนวณโดยใช้แบบ upwind การใช้รหัสของฉันพลังงานการทำงานและขดของการทำแผนที่ลดลงอย่างเหมาะสมสำหรับการทำซ้ำสองสามครั้ง (จากไม่กี่สิบถึงสองสามพันขึ้นอยู่กับขั้นตอนเวลา) แต่หลังจากนั้นการจำลองจะระเบิดขึ้น: พลังงานจะเพิ่มขึ้นถึง NAN ในการวนซ้ำน้อยมาก ฉันลองคำสั่งซื้อหลายรายการสำหรับความแตกต่างและการผสานรวม (สามารถแทนที่คำสั่งที่สูงกว่าเพื่อ cumptrapz ได้ที่นี่ ) และแผนการแก้ไขที่แตกต่างกัน แต่ฉันได้รับปัญหาเดียวกันเสมอ (แม้ในภาพที่ราบรื่นมาก ทุกคนจะสนใจดูรหัสและ / หรือปัญหาทางทฤษฎีที่ฉันกำลังเผชิญอยู่หรือไม่ รหัสค่อนข้างสั้น รหัสที่มีฟังก์ชั่นการดีบัก ฟังก์ชั่นการลงทะเบียน …

11 pde  finite-difference  matlab  fluid-dynamics  image-processing