Maplevn2008's Blog

Việc lập trình tính từng bước đòi hỏi viết code một cách công phu và cũng mất nhiều thời gian.Bạn có thể tự mình lập trình, nếu có vấn đề cần hỏi thì bạn có thể post lên ở đây để chúng ta cùng trao đổi.
Chú ý: Khi lập trình tính từng bước với việc sử dung Maplet thì nhất thiết bạn phải sử dụng đến 2 lệnh Get và Set trong gói Maplets[Tools]. Bạn có thể xem hướng dẫn 2 lệnh này trong help.

Mình quên nhắc là lệnh
> dk:=ImportVector(”C:/data.txt”,skiplines=1);
ở phần trả lời trước thì dk ở đây có kiểu Vector. Chú ý: Trong lệnh QPSolve thì điều kiện phải có dạng tập hợp (set) hay dạng danh sách (list), do đó bạn phải thay đổi kiểu của dk, chẳng hạn sang kiểu list:
>dk:=convert(dk,list);
Với ví dụ của bạn: file C:\dk.txt chứa các điều kiện nằm trên các hàng khác nhau
>with(Optimization):
>f:=x1^2+x2^2;
>constr:=ImportVector(”C:/dk.txt”);
>constr:=convert(constr,list);
>K:=QPSolve(f,constr);

Mình chỉ đưa ra một lời giải tham khảo cho câu 1 và câu 3. Câu 2 sẽ dễ dàng hơn nếu bạn cho hàm số cụ thể.
Câu 1:
> y:=(m-2)*x^3-m*x+2;
> pt:=factor(y);
> pt1:=op(1,pt);
pt2:=op(2,pt);
Trường hợp 1: x=1 là nghiệm của pt2.
> subs(x=1,pt2);
solve(%,{m});
> plot(subs(m=3,y),x=-1..2);
Trường hợp 2: pt2 có nghiệm kép.
> a:=coeff(pt2,x,2);
b:=coeff(pt2,x,1);
c:=coeff(pt2,x,0);
> Delta:=b^2-4*a*c;
> solve(Delta,{m});
Tất nhiên ta chỉ xét m khác 2.
> plot(subs(m=-6,y),x=-1..0);
Câu 3:
> for i from 1159 to 1200 do
if isprime(i) then print(i); fi;
end;

Bạn có thể tham khảo ở đây:
http://www.adeptscience.co.uk/products/mathsim/maple/Demos/java/JavaVBOpenMaple1.html
http://www.nd.edu/~hpcc/solaris8_usr_local/src/maple12/samples/OpenMaple/msvb/
Ngoài ra bạn nên xem thêm các ví dụ trong Help của Maple
>?OpenMaple/VB
Nói thêm : Những câu hỏi liên quan đến việc kết nối Maple với các ngôn ngữ lập trình khác các bạn có thể comment ở trang Maple/Java để mọi người tiện theo dõi.

Mình không hiểu ý bạn hỏi cho lắm. Vấn đề biểu thị kí hiệu tập rỗng cũng như dấu hỏi bạn có thê xử lý như ví dụ sau:
>S:=[``,``,``,``,``,``];# dấu nháy đơn nghiêng.
>G:=[`?`,`?`,`?`,`?`,`?`];
Tuy nhiên theo mình vấn đề chính không phải ở đây, mà ở chỗ phép toán so sánh >=. Để có thể sử dụng phéo toán so sánh >= cho trường hợp so sánh 2 đối tượng đặc biệt mà không làm mất “tính năng so sánh ban đầu” của nó, bạn có thể quan sát ví dụ sau. Ví dụ này cho phép sio sánh 2 list dữ liệu theo giá trị đầu tiên của nó.
> Sosanh:=module()
option package;
export `>=`;
`>=`:=proc(x::list,y::list)
option overload;
if evalb(x[1]>=y[1]) then true;
else false;
end if;
end;
end:
> with(Sosanh);
> [4,3,4]>=[3,4];
> [2,3,4]>=[3,4];
> 2>=3;evalb(%);
Đối với trường hợp của bạn bạn phải định nghĩa phép toán so sánh >= chặt chẽ hơn, chẳng hạn như trường hợp có cùng dấu chấm hỏi , độ dài 2 list dữ liệu khác nhau, có phải luôn so sánh được không,…Sau đó bạn có thể tự mình tiến hành lập trình là xong.

Chú ý:
- dấu nháy kép thẳng đứng (đây là lỗi do WordPress tạo ra nên khi copy lệnh bạn phải sửa lại).
- khi tạo ra file data.txt trong ổ C: thì giữa các hàng dữ liệu không có dòng trống.
Trong Maple có gói Statistics làm việc với dữ liệu thống kê cũng như tính toán đối với đại lượng ngẫu nhiên 1 chiều.
> with(Statistics):
Ví dụ với 2 mẫu dữ liệu:
> X:=[2,3,1,2,3,4,5,3,4,5,6,7,8]:
> Y:=[7,8,9,6,5,7,8,0,3,2,1,4,6]:
- giá trị trung bình:
> mX:=Mean(X);mY:=Mean(Y);
-phương sai mẫu(ko chệch)
> Variance(X);
- ma trận tương quan ( X,Y phải cùng số phần tử)
> CovarianceMatrix(Matrix([Vector(X),Vector(Y)]));
- ma trận hệ số tương quan ( X,Y phải cùng số phần tử)
> CorrelationMatrix(Matrix([Vector(X),Vector(Y)]));

Rất tiếc là mình không có code. Tuy nhiên chúng ta có thể giải nhờ kết quả đơn giản sau: Giả sử cần tìm ma trận nghích đảo của ma trân A, ta viết ma trận mở rộng [A,I], với I là ma trận đơn vị. Thực hiện các phép biến đổi sơ cấp hàng đưa ma trận mở rộng này về ma trân dạng [I,B]. Lúc này B chính là ma trận nghích đảo của ma trận A.
Các phép biến đổi sơ cấp bạn có thể xem ở đây:
http://maplevn2008.wordpress.com/2009/01/25/phep-toan-tren-ma-tr%E1%BA%ADn-t4/

Bạn có thể xem gợi ý ở bài viết “Sử dụng file “.m”" :
http://maplevn2008.wordpress.com/2009/05/12/s%E1%BB%AD-d%E1%BB%A5ng-file-m/

Maple chỉ cho xem code các lệnh ở “phần chính” của nó, tức là các thủ tục con bên trong theo mình được biết là không xem được. Do đó đối với các hàm cơ bản thì ta hầu như xem được tất cả code của nó. Bạn quan sát ví dụ sau (chú ý bổ sung thêm tham số):
> interface(verboseproc=2);
> print(cos);
> showstat(cos);
Đối với ví dụ của bạn, bạ có thể viết:
> interface(verboseproc=2);
> print(Optimization[LPSolve]);
Bạn có thể thấy rằng bên trong thủ tục này có chứa vài thủ tục con như Convert:-MatrixForm:-ProcessLP, LP:-SolveInteger , …, nhưng các thủ tục này lại không nằm trong các gói chuyên ngành mà Maple cung cấp.

Đối với một số miền đơn giản thì Maple có thể giúp ta tính trực tiếp. Bạn có thể xem thêm trong Help của Maple.
> with(Student[VectorCalculus]);
- đường loại 1:
> ?PathInt
- đường loại 2:
> ?LineInt
- mặt ( không thấy tham số gì cho phép phân biệt mặt loại nào cả)
> ?SurfaceInt
Ví dụ về tích phân đường:
http://maplevn2008.wordpress.com/2008/11/23/tich-phan-d%C6%B0%E1%BB%9Dng/
Khi có tham số “inert” thì kết quả là công thức dạng tích phân bội. Tuy nhiên trong một số trường hợp mà ta không thêm tham số này nhưng do tích phân bội thu được không thể tính ra kết quả giải tích thì kết quả trả về cũng là một biểu thức tính phân. Lúc này bạn có thể dùng lệnh evalf để lấy giá trị xấp xỉ của nó, ví dụ:
> kq:=int(exp(-x^3)*sin(x^2),x=0..1);
> evalf(kq);

Theo mình được biết có 1 cách giấu code đi là bạn tạo package nhờ module. Các hàm ở phần export bạn không tiến hành định nghĩa trực tiếp mà dùng một biến local để xây dựng hàm có chức năng tương tự, rồi sau đó chỉ việc gọi ra trong thủ tục chính. Đây là một vấn đề hay, cuối tuần này mình sẽ post một ví dụ minh họa.

Bạn có thể dùng hàm time, ví dụ tính thời gian tính giai thừa cho bởi hàm fac:
> fac:=proc(n)
option remember;
if n*(n-1)=0 then return 1;
else return n*procname(n-1); fi;
end;
> time(fac(1000));
Nếu bạn muốn tính 1 lần cho nhiều thủ tục và hàm thì có thể làm như sau:
> t:=time(): # lưu thời điểm đầu
fac(2000);
evalf(int(x^x*sin(x),x=1..3));
time()-t; # thời gian tính bằng thời điểm kết thúc trừ đi thời điểm ban đầu

Điều này là hoàn toàn bình thường. Tốc độ tính toán phụ thuộc vào nhiều yếu tố như tài nguyên của máy ở thời điểm hiện tại (số lượng quá trình đang hoạt động, dung lượng RAM còn rãnh,…), cấu trúc lệnh sử dụng( một số option có thể giúp ta trả về cho hàm một bảng nhớ các giá trị, điều này cho phép rút ngắn thời gian ở các bước tiếp theo,…),…
Bạn xài Maple 12 thì chạy bình thường các file viết trên Maple 9. Tuy nhiên các file viết trên Maple 12 thì có khả năng không chạy được trên Maple 9 do có nhiều hàm ,thủ tục và bố cục giao diện mà chỉ có trong phiên bản mới.

Bạn có thể xem hướng dẫn ở đây để biết thêm cách thiết lập sự kiện khi kích vào Button:
http://maplevn2008.wordpress.com/2009/06/18/maplet-m%E1%BB%99t-vai-s%E1%BB%B1-ki%E1%BB%87n-khi-kich-button/

Các package trong Maple đều có phần giới thiệu tổng quan ở phần mở đầu( trong Help), ví dụ:
>?Student
>?LinearAlgebra
>?Maplets
…………………………………

Các ví dụ minh họa sẽ được trình bày trong từng lệnh cụ thể của package.
Một vài thuật ngữ toán các bạn có thể vào Dictionary trong Maple để tra cứu
Help-> Manuals,Dictionary and more->Dictionary.

Để tính vi phân bạn có thể xem ví dụ trong Help:
>?difforms/d

Chào bạn. Vấn đề lý thuyết đồ thị thì mình chỉ biết sơ sài thôi. Một vài gợi ý có tính chất chủ quan của mình bạn có thể tham khảo:
- Nếu bạn có kiến thức tốt về Lý thuyết đồ thị và có kĩ năng lập trình thì có thể viết lại một vài lệnh như : tìm đường đi ngắn nhất giữa hai đỉnh với ma trận liền kề, có thể in ra từng bước và minh họa đồ thị, hay sử dụng các định lý kiểm tra đồ thị đã cho có chu trình Haminton, Euler,… và minh họa đồ thị trong trường hợp tìm thấy,….
- Nếu khả năng lập trình hạn chế thì bạn có thể sử dụng các hàm dựng sẵn, lúc này chủ yếu là bạn cho người dùng biết những thông tin có thể về đồ thị cũng như minh họa hình học của nó ( có thể xây dựng giao diện minh họa cho đẹp mắt,…),…

Chào bạn.
Bạn có thể tham khảo cách tổ chức dữ liệu về đỉnh, cạnh trong gói GraphTheory. Chẳng hạn đối với đồ thị vô hướng ta sử dụng kiểu dữ liệu tập hợp (set) để khai báo tập đỉnh, các cạnh. Còn đối với đồ thị có hướng thì ta sử dụng kiểu dữ liệu list. Việc xét tính liên thông, tìm đường đi,… là vấn đề về thuật toán, bạn phải tự tìm hiểu thôi.

Chào bạn.
Bạn có thể sử dụng gói Optimization.
> with(Optimization);