จำนวน 194 ความเห็น, หน้า่ | 1| 2| 3| 4| 5| 6| 7| 8| -9- 10|

ความเห็นเพิ่มเติมที่ 158 12 ก.ค. 2550 (09:44)
สำหรับคำถามเกี่ยวกับ การใช้ภาษา C \ C++ ร่วมกับ Scilab
ผมลองค้นดูได้คำแนะนำดังนี้

(ANSI-) C

A simple C function meting out signature requirements has e.g. this shape:

#include
double
fgen(const double *x)
{
if (*x > 0.0)
return 1.0 / sqrt(*x);
else
return 0.0;
}
After compilation (e.g. cc -c fgen.c) the compiled code can be linked to Scilab and called with the integration routine.

link('fgen.o', 'fgen', 'c');
[ires, ierr, neval, info] = intgen(0.0, 1.0, 'fgen')


C++
A C++ source for a function could look like this:

#include

extern "C" {
double fgk(const double *x);
}

double
fgk(const double *x)
{
return 2.0 / (2.0 + sin(10.0 * M_PI * (*x)));
}
After compilation (e.g. c++ -c fgen.c) the compiled code can be linked to Scilab and called with the integration routine.

link('fgk.o', 'fgk', 'c');
[ires, ierr, neval, info] = intgk(0.0, 1.0, 'fgk', 0, %eps, '15-31')

MathGuy
ร่วมแบ่งปันความรู้และความเห็นแล้ว 1988 ครั้ง - ได้รับดาวแล้ว 323 ดวง - โหวตเพิ่มดาว

ความเห็นเพิ่มเติมที่ 159 12 ก.ค. 2550 (15:58)
ขอบคุณครับคุณ MathGuy สำหรับ function plot2d3 ผมได้ลองพิจารณาอยู่ก่อนแล้วครับ
จุดประสงค์ของการใช้คำสั่ง stem ของ MATLAB เพื่อไว้ plot graph of discreate value in the DSP ในหนังสือตัวอย่างเกี่ยวกับ DSP จะใช้คำสั่งนี้เกือบทั้งหมด ผมเลยอยากที่จะสามารถใช้คำสั่งที่สามารถทำแบบเดียวกันได้ใน SCILAB( ผมใช้เครื่องบริษัทหนะครับ เลยลงได้แต่ฟรีโปรแกรม)
ข้างล่างเป็นตัวอย่างจากหนังสือ ที่ผมกำลังศึกษาอยู่

% Chapter 02: Example 02.01: Signal Synthesis
%
figure(1); clf
% a) x(n) = 2*delta(n+2) - delta(n-4), -5<=n<=5
n = [-5:5];
x = 2*impseq(-2,-5,5)-impseq(4,-5,5);
subplot(2,2,1); stem(n,x); title('Sequence in Example 2.1a')
xlabel('n'); ylabel('x(n)'); axis([-5,5,-2,3])
%
% b) x(n) = n[u(n)-u(n-10)]+10*exp(-0.3(n-10))(u(n-10)-u(n-20)); 0<=n<=20
n = [0:20];
x1 = n.*(stepseq(0,0,20)-stepseq(10,0,20));
x2 = 10*exp(-0.3*(n-10)).*(stepseq(10,0,20)-stepseq(20,0,20));
x = x1+x2;
subplot(2,2,2);stem(n,x);
title('Sequence in Example 2.1b')
xlabel('n');ylabel('x(n)');axis([0,20,-1,11])
%
% c) x(n) = cos(0.04*pi*n) + 0.2*w(n); 0<=n<=50, w(n): Gaussian (0,1)
n = [0:50];
x = cos(0.04*pi*n)+0.2*randn(size(n));
subplot(2,2,3);stem(n,x);title('Sequence in Example 2.1c')
xlabel('n');ylabel('x(n)');axis([0,50,-1.4,1.4])
%
% d) x(n) = {...,5,4,3,2,1,5,4,3,2,1,...}; -10<=n<=9
n=[-10:9];
x=[5,4,3,2,1];
xtilde=x' * ones(1,4);
xtilde=(xtilde(:))';
subplot(2,2,4);stem(n,xtilde);title('Sequence in Example 2.1d')
xlabel('n');ylabel('xtilde(n)');axis([-10,9,-1,6])

ผลลัพธ์ครับ

ความเห็นเพิ่มเติมที่ 160 12 ก.ค. 2550 (17:17)
จากคำสั่ง plot2d3 ของ scilab

ที่เราต้องปรับ คือ ที่จุดปลายของต้นไม้ หรือกิ่งเล็กๆแต่ละอัน เราต้องหาวิธีให้เขียนวงกลมรูปเล็กๆไว้ที่ตรงปลาย

เท่าที่ผมค้นๆดูในคำสั่งทำ graphic ของ scilab
โดยความรู้สึก(ของคนที่ยังเป็นมือใหม่ ยังไม่ใช่ expert) ผมว่า ตัวอย่าง คำอธิบายต่างๆ ยังดูเข้าใจยาก

บางทีเราต้องให้เวลาที่จะเข้าไปเรียนรู้เองพอสมควรเหมือนกันครับ

ความเห็นเพิ่มเติมที่ 161 13 ก.ค. 2550 (15:32)
จาก คห. 158 ไม่เข้าใจเลยครับ

รบกวนอธิบายทีนะคับ

ความเห็นเพิ่มเติมที่ 162 13 ก.ค. 2550 (16:10)
ผมพยายามแปลง code Matlab ในการสร้างรูป Fractol ให้เป็น SCilab กรุณาช่วย
ดูและให้คำแนะนำด้วยครับ


function M=Mandelbrot(re,im,m,tol,ma);
// MANDELBROT computes and displays the Mandelbrot set
// M=Mandelbrot(re,im,m,tol,ma); computes the Mandelbrot set defined
// by the behavior of the iteration x(n+1)=x(n)^2+nu whre nu and x(n
// are complex numbers. Starting with x(0)=0 we compute all values
// nu on a grid [re(1),re(2)]x[im(1),im(2)] with grid points in each
// direction, the iteration stops if |x(n)|>tol or n=ma
// and then n defines a color of the pixel at location x(0) in the result M.

regrid(m)=0;
imgrid(m)=0;
M(m)=0;
M(m,m)=0;

// construct grid
for j=1:m,
regrid(j)= re(1) + (j-1)*(re(2)-re(1))/(m-1);
imgrid(j)= im(1) + (j-1)*(im(2)-im(1))/(m-1);
end

realpart=0;
imagpart=0;

for j=1:m,
for k=1:m,
realpart=regrid(j);
imagpart=imgrid(k);

// initialize variables for iteration
nu(1)=realpart;
nu(2)=imagpart;
re_xn=0;
im_xn=0;
re_xn1=0;
im_xn1=0;
n=0;

while sqrt(re_xn*re_xn+im_xn*im_xn)<=tol | n<ma;
// splitted into real and imaginary parts
re_xn1=(re_xn*re_xn)-(im_xn*im_xn)+nu(1);
im_xn1=(2*re_xn*im_xn)+nu(2);
re_xn=re_xn1;
im_xn=im_xn1;
n=n+1;
end;

M(j,k)=n;

end;
end;

// show picture

endfunction

ความเห็นเพิ่มเติมที่ 163 19 ก.ค. 2550 (14:50)
นำไฟล์เสียง มา วาดเป็นกราฟได้รึป่าวคับ

ความเห็นเพิ่มเติมที่ 164 20 ก.ค. 2550 (16:05)
code fractal จาก MatlaB ไปเป็น scilab

function Mandelbrot(iter,pixel)

// mandelbrot uses iter = 23, pixel = 400.
// MANDELBROT(ITER,PIXEL) draws the Mandelbrot set;
// ITER gives the number of iterations and
// PIXEL indicates the number of points on the x-axis.

// MANDELBROT uses iter = 23, pixel = 400.

// This file was generated by students as a partial fulfillment
// for the requirements of the course "Fractals", Winter term
// 2004/2005, Stuttgart University.

// Author : Jochen Schurr
// Date : Nov 2004
// Version: 1.0

// default setting


// decompose x- and y-axis according to the chosen
// proportion of r = 3:4.
// ใน scilab ต้องกำหนด i ด้วย มันไม่รู้จัก
i=(-1)^(0.5)
r = 3/4;
x = linspace(-2.5,1.5,pixel);
y = linspace(-1.5,1.5,round(pixel*r))';

// create the matrix C containing all pixel
// to be analysed.
[Re,Im] = meshgrid(x,y);
C = Re + i * Im;

// The matrix B assigns a number to every pixel which is
// characteristic with respect to the velocity of divergence.
// Therefore B and C have the same dimension.
B = zeros(round(pixel*r),pixel);

Cn = B; // C0 = 0+0i
for l = 1:iter
Cn = Cn.*Cn + C; // performing the map
B = B + (abs(Cn)<2); // if |cn| > 2 the sequence diverges
end;

// plot settings
Matplot(B);

endfunction

ความเห็นเพิ่มเติมที่ 165 20 ก.ค. 2550 (16:12)
นี่เป็นอีก code fractal หนึ่งที่แปลงได้ครับ ตอน run ต้องกำหนด i ด้วยนะครับ

function Julia(c,k,v)

// Julia uses c = 0.2+0.65i, k = 14, v = 500.
// JULIA(C,K,V) draws the Julia set with the following parameters:
// c is a complex number used in the map f(z) = zฒ + c.
// k gives the number of iterations
// v determines the number of points on the x-axis.

// JULIA uses c = 0.2+0.65i, k = 14, v = 500.

// This file was generated by students as a partial fulfillment
// for the requirements of the course "Fractals", Winter term
// 2004/2005, Stuttgart University.

// Author : Sylvia Frey
// Date : Nov 2004
// Version: 1.1

// default settings
if rand() < 3
c = 0.2+0.65*i;
k = 14;
v = 500;
end

// radius of the circle beyond which every point diverges
r = max(abs(c),2);

// divide the x-axis
d = linspace(-r,r,v);

// create the matrix A containing complex numbers
A = ones(v,1)*d+i*(ones(v,1)*d)';

// create the point matrix
B = zeros(v,v);

// iteration
for s = 1:k
B = B+(abs(A)<=r);
// the map
A = A.*A+ones(v,v).*c;
end;

// plot settings
Matplot(B);
xset("colormap",jetcolormap(64))


endfunction

ความเห็นเพิ่มเติมที่ 166 20 ก.ค. 2550 (17:41)
code ใน คหพ#164 เป็น code ที่เขียนได้กระทัดรัดมาก

ผมปรับให้เป็นดังนี้

function MDB()

x=linspace(-2.5,1.5,400);
y=linspace(-1.5,1.5,300);

[Re,Im]=meshgrid(x,y);
C=Re+%i*Im;

B=zeros(300,400);

Cn=B;

for j=1:23
Cn=Cn.*Cn+C;
B=B+(abs(Cn)<2);
end

Matplot(B)

endfunction

เมื่อ run แล้วได้ผลดังนี้

ความเห็นเพิ่มเติมที่ 167 25 ก.ค. 2550 (21:05)
คัดลอกมาจากกระทู้ : http://www.vcharkarn.com/include/vcafe/showkratoo.php?Pid=101548

โดยคุณ gte723s
...................................................................................................................
SCILAB เป็นโปรแกรมภาษาขั้นสูงที่ถูกพัฒนาขึ้นโดยความร่วมมือกันระหว่างนักวิจัยจากสถาบัน Institut National De Recherche En Informatique Et En Automatique (INRIA) และ École nationale des ponts et chausses (ENPC) ประเทศฝรั่งเศส ตั้งแต่ปี ค.ศ. 1990 โดยมีจุดมุ่งหมายเพื่อใช้ในการคำนวณเชิงตัวเลขและแสดงผลกราฟริกที่ซับซ้อน นอกจากนี้โปรแกรม SCILAB ยังเป็นโปรแกรมที่ให้ฟรี (ไม่ต้องเสียเงินค่าลิขสิทธิ์ซอฟต์แวร์) ผู้อ่านสามารถดาวน์โหลดได้จาก http://www.scilab.org

วิธีการใช้งานโปรแกรม SCILAB ไม่ยากครับ เมื่อเข้าใจแล้วก็สามารถที่จะนำมาใช้ในการแก้ไขปัญหาทางคณิตศาสตร์ ฟิสิกส์ และวิศวกรรม ได้

ตัวอย่างเช่น การหาตัวคูณร่วมน้อยและตัวหารร่วมมาก

ถ้ากำหนดให้ p คือเวกเตอร์ของเลขจำนวนเต็มบวก (positive integer) หรือเวกเตอร์ของพหุนาม (polynomial) โปรแกรม SCILAB สามารถคำนวณหาตัวคูณร่วมน้อยและตัวหารร่วมมากได้โดยใช้คำสั่ง lcm และ gcd ตามลำดับ ตัวอย่างเช่น

-->p = int16([12 15 30]) //ทำให้เป็นเลขจำนวนเต็มขนาด 16 บิต

-->pp = lcm(p) //นั่นคือ ตัวคูณร่วมน้อยของ 12, 15, และ 30 คือ 60
pp =
60

-->qq = gcd(p) //นั่นคือ ตัวหารร่วมมากของ 12, 15, และ 30 คือ 3
qq =
3

นอกจากนี้คำสั่ง lcm และ gcd ยังสามารถใช้หาตัวคูณร่วมน้อยและตัวหารร่วมมากพหุนามได้เช่นกัน ตัวอย่างเช่น

-->x = poly(0, 'x'); //กำหนดให้ x เป็นตัวแปรพหุนาม

-->p = [x^2 – 5 * x - 14, x^2 – 12 * x + 35] //พหุนาม และ
p =
2 2
- 14 - 5x + x 35 - 12x + x

-->lcm(p)
ans =
2 3
70 + 11x - 10x + x

-->gcd(p)
ans =
- 7 + x

ความเห็นเพิ่มเติมที่ 168 23 ส.ค. 2550 (09:37)
Website แนะนำการใช้งานโปรแกรม SCILAB เบื้องต้น http://home.npru.ac.th/~t3058/Scilab.html

ความเห็นเพิ่มเติมที่ 169 23 ส.ค. 2550 (21:05)
ขอบคุณ คุณ gte723s ที่ให้ link แนะนำการใช้งานโปรแกรม ที่ดีมากครับ

ความเห็นเพิ่มเติมที่ 170 26 ส.ค. 2550 (15:49)
ขอโปรกแรม CAI graph หน่อย

ความเห็นเพิ่มเติมที่ 171 31 ส.ค. 2550 (21:13)
ช่วยหนูทีนะคะ..พอดีหนูทำสัมมนาเรื่องการวาดกราฟด้วย SCILAB ค่ะที่มหาวิยาลัย
โดยหนูต้องตอบคำถามอาจารย์เมื่อพูดเสร็จค่ะ..แต่หนูตอบคำถามอาจารย์ไม่ได้อยู่หลายคำถามทีเดียวค่ะ...จึงขอรบกวนอาจารย์ช่วยตอบคำถามหน่อยนะคะ...ขอบคุณล่วงหน้าค่ะ (หนูต้องตอบคำถามอาจารย์ในวันที่ 7/9/50 ค่ะ)
- SCILAB ย่อมาจากอะไรคะ..(หนูหาไม่เจอค่ะ..เลยต้องรบกวนอาจารย์ค่ะ)
- plot(y) ที่กำหนดพารามิเตอร์ x = 1 จะวาดกราฟแล้วเหมือนกับ plot (1,y) มั้ยคะ
คือหนูทดลองปรากฏว่ามันต่างกัน หนูเกรงว่าหนูจะใช้คำสั่งผิดค่ะ เลยลองถามอาจารย์ดูค่ะ
- อาจารย์คะ ใน SCILAB มีสีทั้งหมด 32 สีใช่มั้ยคะ
(กลัวจะตอบอาจารย์ผิดอ่ะค่ะเลยลองถามอาจารย์ปิยะดูค่ะเพื่อความมั่นใจ)
- อาจารย์คะ...nb_colors มีค่าสูงสุดที่ 30 ใช่มั้ยคะ เพราะว่าหนูใส่เยอะกว่านี้จะ error ค่ะ

ความเห็นเพิ่มเติมที่ 172 3 ก.ย. 2550 (17:40)
ถ้าต้องการเขียนภาษา c++ โจทย์ y = x^n เขียนยังไงค่ะ รบกวนท่านผู้รู้ช่วยด้วยค่ะ

ความเห็นเพิ่มเติมที่ 173 4 ก.ย. 2550 (20:20)
มันแสดงผลยังไงอ่ะ

ความเห็นเพิ่มเติมที่ 174 31 ต.ค. 2550 (15:18)
เขียน scilab ไม่เป็นทำไงดีค่ะ หนู่ได้ทำโปรเจ็กคะ Egyptian fraction คะ
ต้องเริ่มยังไงคะ

ความเห็นเพิ่มเติมที่ 175 31 ต.ค. 2550 (20:08)
ตัวอย่างการหา 1/x + 1/y +1/z + 1/w =1


function EFraction(n)
for i=2:n
for j=(i+1):n
for k=(j+1):n
for m=(k+1):n
if ((j*k*m+i*k*m+i*j*m+i*j*k)==(i*j*k*m)) then
x=[i j k m]; disp(x);
end
end
end
end
end
endfunction

ผลการ run-->EFraction(100)

2. 3. 7. 42.

2. 3. 8. 24.

2. 3. 9. 18.

2. 3. 10. 15.

2. 4. 5. 20.

2. 4. 6. 12.

-->

ความเห็นเพิ่มเติมที่ 176 1 พ.ย. 2550 (12:43)
splitting method
สูตร 1/x = 1/x+1 + 1/x(x+1)
ให้ p/q เป็นเศษส่วนอย่างตำ โดยที่ p และ q และ p/q น้อยกว่า 1
ถ้าเขียน scilab ทำยังไงคะ
ขอบคุณคะ

ความเห็นเพิ่มเติมที่ 177 2 พ.ย. 2550 (14:42)
ต้องเริ่มต้นวิเคราะห์ว่า

1) เรามีอะไร เรากำหนดอะไรก่อนบ้าง

2) เราต้องการทำให้ได้อะไร จากข้อ (1)


ตอบคำถาม 2 ข้อนี้ให้ได้ก่อนนะครับ ถึงจะเริ่มเขียนโปรแกรมได้