Đến nội dung


Chú ý

Nếu các bạn đăng kí thành viên mà không nhận được email kích hoạt thì hãy kiểm tra thùng thư rác (spam). Nếu không biết cách truy cập vào thùng thư rác thì các bạn chịu khó Google hoặc đăng câu hỏi vào mục Hướng dẫn - Trợ giúp để thành viên khác có thể hỗ trợ.


Hình ảnh

Phân tích bình phương SOS cho bậc $3$ bằng Maple.

sos deg 3 stsos

  • Please log in to reply
Chưa có bài trả lời

#1 tthnew

tthnew

    Thượng sĩ

  • Thành viên
  • 239 Bài viết
  • Giới tính:Nam
  • Đến từ:Nơi cần đến.
  • Sở thích:Viết blog, viết SOS .v.v.. etc.

Đã gửi 09-04-2020 - 06:52

Có nhiều bạn hỏi mình về cách phân tích bình phương $S.$$O.$$S$ trên Maple như thế nào. Mình cũng chỉ mới nhập môn Maple thôi nên kiến thức còn hạn hẹp, nên mong mọi người đóng góp thêm những "code" khác để chạy các bậc cao hơn hoặc tối ưu hơn code mình chăng$?$

 

Code này mình viết cho bậc $3$, các bạn có thể tham khảo thêm các bậc khác trên Epsilon $10$, Epsilon $13$.

print("============================================================================");
print("Chuong trinh duoc viet boi tth_new");
print("version:1.0");
print("Copyright 2020-2021");
print("============================================================================");
sgm := proc(f) 
global f1, f2, f3; 
f1 := f; 
f2 := subs({a = b, b = c, c = a,x=y,y=z,z=x}, f1); 
f3 := subs({a = b, b = c, c = a,x=y,y=z,z=x}, f2); 
subs(f1 + f2 + f3); end proc:
pro := proc(f) 
global f1, f2, f3; 
f1 := f; 
f2 := subs({a = b, b = c, c = a}, f1); 
f3 := subs({a = b, b = c, c = a}, f2); 
subs(f1*f2*f3); end proc:
sos := proc(f) 
local A, B, C, J, temp; 
A := coeff(subs({a = a, b = 1, c = 0}, numer(f)), a^3); 
B := coeff(subs({a = a, b = 1, c = 0}, numer(f)), a^2); 
C := coeff(subs({a = a, b = 1, c = 0}, numer(f)), a); 
J := coeff(subs({a = a, b = 1, c = 1}, numer(f)), a) - B - C; 
temp := sgm((a*m[1] + b*m[2] + c*m[3])*(a - b)^2); 
if subs({a = 1, b = 1, c = 1}, f) = 0 then print(`SOS form is:`*subs({m[1] = 1/2*A - 1/6*B + 1/6*C, m[2] = 1/2*A + 1/6*B - 1/6*C, m[3] = 1/2*A + 1/2*B + 1/2*C}, temp)); 
else print(`False!`); end if; end proc:

Các bạn có thể nhập thẳng đoạn code trên vào Worksheet của Maple, hoặc có thể lưu lại vào Notepad.

Ở đây mình lưu file Notepad vào thư mục SOS ở Desktop với tên là sostool-deg3. Sau đây là một vài ví dụ:

 

is7m0J5.png

Chú ý rằng, các bất đẳng thức mà bạn muốn phân tích $SOS$ về dạng $\sum S_a (b-c)^2$ như trên phải có một trong những đẳng thức là xảy ra tại $a=b=c$. 

Nếu không phần mềm sẽ báo "False!"

Bài viết mình chỉ vọn vẹn nhiêu đó, mong mọi cùng đóng góp và chia sẻ các đoạn "code" trên Maple khác vào Topic này nhé!


Bài viết đã được chỉnh sửa nội dung bởi tthnew: 09-04-2020 - 06:56






Được gắn nhãn với một hoặc nhiều trong số những từ khóa sau: sos, deg 3, stsos

1 người đang xem chủ đề

0 thành viên, 1 khách, 0 thành viên ẩn danh