II. TỔNG QUAN VỀ AMÔN SUNPHÁT
2. CƠ SỞ THỰC NGHIỆM
3.3. Đánh giá về khả năng áp dụng vào sản xuất thực tế ở quy mô công nghiệp
3.3.8. Địa điểm áp dụng
Nếu dự án được triển khai ở quy mô công nghiệp, địa điểm dự kiến xây dựng nhà máy có thể lựa chọn ở khu vực lân cận hoặc ở trong mặt bằng nhà máy sản xuất
DAP Đình Vũ Hải Phịng vì các lý do sau:
- Có thể tận dụng được nguồn nguyên liệu tại chỗ, tiết kiệm chi phí vận chuyển.
có khả năng cơ giới hóa nguồn nguyên liệu trong sản xuất. + Gyps: Lấy ở bãi thải của Nhà máy DAP gần đó
+ Khí NH3: Cung cấp từ các nhà máy sản xuất NH3 được vận chuyển bằng đường bộ hoặc đường thủy về kho chứa trong nhà máy
+ Khí CO2: Nguyên liệu được cung cấp từ nhiều nguồn như: các nhà máy sản xuất đạm hoặc nhà máy tổng hợp etannol vận chuyển bằng đường bộ hoặc
đường thủy về kho chứa trong nhà máy. Khí CO2 cũng có thể xem xét tận dụng
từ nguồn khí thải của lị hơi đốt than đem xử lý tách lấy CO2 giống như dây chuyền đang hoạt động tại nhà máy đạm Phú Mỹ.
- Có thể sử dụng một số hệ thống phụ trợ của nhà máy DAP nhằm tiết kiệm chi
phí đầu tư
- Phù hợp với quy hoạch tập trung khu cơng nghiệp sản xuất hóa chất, phân bón. Dễ dàng quản lý trong vận hành sản xuất và mơi trường.
tại KCN Đình Vũ – Hải Phịng”
Học viên: Lưu Ngọc Vĩnh Trang 78
Trên thực tế việc sản xuất SA từ Gyps ở quy mơ cơng nghiệp có thể được lựa chọn với công suất 200.000 tấn/năm. Tuy nhiên, trước khi sản xuất ở quy mô công
nghiệp cần lập dự án sản xuất thử nghiệm (có thể lựa chọn công suất sản xuất SA là 2.000kg/h) để đánh giá tính khả thi của đề tài khi áp dụng trong sản xuất với các điều kiện trong môi trường thực tế.
Học viên: Lưu Ngọc Vĩnh Trang 79
KẾT LUẬN
***********
Đề tài nghiên cứu thăm dò khả năng thu hồi và sử dụng hiệu quả hợp chất lưu
huỳnh trong Gyps bằng phản ứng tạo thành Amơn Sunphát, từ đó xác định các
thông số công nghệ ảnh hưởng đến hiệu suất phản ứng tạo thành Amôn Sunphát. - Đề tài có tính thực tiễn cao do sản phẩm Amơn Sunphát hiện nay Việt Nam gần như hoàn toàn phải nhập khẩu, việc nghiên cứu công nghệ sử dụng nguyên liệu là
chất thải Gyps thải để sản xuất Amôn Sunphát nếu thành công sẽ là một bước đột
phá để xây dựng nhà máy sản xuất Amôn Sunphát cho tập đồn Hóa chất Việt Nam
- Mặt khác, khi chuyển sang sản xuất ở quy mơ cơng nghiệp có ý nghĩa cao về sử
dụng hiệu quả nguồn thải Gyps để sản xuất các sản phẩm có ích.. Tuy nhiên, khi tiến hành áp dụng từ nghiên cứu phịng thí nghiệm sang cơng nghiệp cần phải qua mơ hình pilot để kiểm nghiệm lại các kết quả nghiên cứu đã đạt được trước khi đưa vào sản xuất ở quy mơ cơng nghiệp. Do vậy đề tài cần có những bước nghiên cứu
tiếp ở quy mô pilot để đánh giá chính xác về các thơng số ảnh hưởng đến quá trình phản ứng tạo thành SA trước khi chuyển sang quy mô công nghiệp
Phạm vi nghiên cứu của đề tài:
- Trên cơ sở thu thập tài liệu tham khảo ngoài nước kết hợp với khảo sát sơ bộ trong quá trình nghiên cứu thực nghiệm, hướng nghiên cứu của đề tài chỉ tập vào các nội dung sau:
- Khảo sát nhu cầu sản xuất và tiêu thụ Gyps, SA trên thế giới và Việt Nam.
- Khảo sát, xác định thành phần Gyps thải của nhà máy DAP Đình Vũ – Hải Phịng,
đánh giá về thành phần nguyên liệu Gyps cho việc nghiên cứu sản xuất SA.
- Đưa ra quy trình sản xuất SA từ Gyps thải và Amôn Cácbonát.
- Xây dựng mơ hình thí nghiệm khảo sát các thơng số công nghệ nhiệt độ, tốc độ khuấy trộn, thời gian phản ứng ảnh hưởng đến hiệu suất tạo thành SA.
- Đưa ra phương trình hồi quy về các thông số công nghệ ảnh hưởng đến hiệu suất phản ứng, xác định được các thông số công nghệ cho hiệu suất cao, hiệu quả kinh tế cao khi áp dụng vào quy mơ cơng nghiệp.
tại KCN Đình Vũ – Hải Phòng”
Học viên: Lưu Ngọc Vĩnh Trang 80
Kết quả của đềtài đã thực hiện:
1. Lựa chọn và đưa ra phương pháp thích hợp, tiên tiến, cho hiệu suất cao, phù hợp với điều kiện nghiên cứu trong nước.
2. Đã xác định hàm lượng hợp chất lưu huỳnh hữu ích trong Gyps thải của nhà máy
DAP Hải Phòng ở dạng CaSO4.2H2O chiếm 75%, các tạp chất khác 25%, phù hợp với sản xuất Amôn Sunphát.
3. Đã xây dựng được phương trình hồi quy thực nghiệm ảnh hưởng của các thông
số công nghệđến hiệu suất phản ứng tạo SA.
4. Đã xác định được quy trình cơng nghệ trong phịng thí nghiệm để điều chế Amơn
Sunphát đi từ Gyps thải với 2 nguyên liệu chính là Gyps thải của Nhà máy DAP Đình Vũ được tinh chế qua các khâu : nghiền, rửa, sấy khô để tách nước, và nguyên
liệu thứ 2 là Amôn Cácbonát (25%); ở khâu phản ứng đã nghiên cứu hiệu suất phản
ứng theo một số thông số chính của q trình là nhiệt độ, tốc độ khuấy, thời gian
phản ứng, đã có kết quả. Hiệu suất phản ứng tạo thành Amôn Sunphát trên 93% ở
điều kiện:
+ Nhiệt độ phản ứng : 45 ± 1 0C
+ Tốc độ khuấy : 45 ± 2 vòng/phút
+ Thời gian phản ứng : 2,5 ± 0,1 (giờ)
Hạn chế của đề tài:
- Chưa tiến hành nghiên cứu ảnh hưởng của thành phần tạp chất đến hiệu suất phản
ứng.
Kiến nghịhướng mở rộng:
- Nghiên cứu q trình phản ứng tạo thành Amơn Sunphát chịu ảnh hưởng của
nhiều quá trình phức tạp. Từ quá trình phản ứng giữa Amơniắc với Cácboníc tạo
thành Amơn Cácbonát đến q trình phản ứng Amơn Cácbonát với bã thải Gyps tạo thành Amơn Sunphát.
- Do kinh phí cịn hạn hẹp, thời gian có hạn, việc nghiên cứu còn nhiều hạn chế, mới chỉ nghiên cứu ở quy mơ phịng thí nghiệm do đó cần cấp thêm kinh phí và thời gian để nghiên cứu các thông số ảnh hưởng rộng hơn, phản ánh đầy đủ hơn khi
Học viên: Lưu Ngọc Vĩnh Trang 81
TÀI LIỆU THAM KHẢO Tiếng Việt
[1] Nguyễn Minh Tuyển, Phạm Văn Thiêm. Kĩ thuật hệ thống công nghệ hóa học, tập 1. Nhà xuất bản khoa học và kĩ thuật, Hà Nội, 1997
[2] Nguyễn Minh Tuyển – Quy hoạch thực nghiệm- Nhà xuất bản khoa học kỹ thuật Hà Nội, 2005.
[3] Nguyễn Bin, Các quá trình, thiết bị trong cơng nghệ hóa chất và thực phẩm, Tập 5, Nhà xuất bản Khoa học & Kỹ thuật Hà nội
[4] Đào Trần Cao, Giáo trình Vật lý Chất rắn cơ sở, Viện Khoa Học Vật liệu, TTKHTN & CNQG, Hà nội, 10/2001, trang II-1 đến II-17.
[5] Chu Thị Thơm, Kỹ thuật sản xuất chế biến phân bón – Nhà xuất bản lao động - 2006
[6] Thuyết minh dự án đầu tư xây dựng nhà máy sản xuất DAP Đình Vũ - Hải
Phịng cơng suất 330.000 tấn/năm
Tiếng Anh
[7] Arab Gulf, Scient. Res. Math.Phys. A6 (l). PP77-83(1988),"Chemical Processing of Phốtphátghypsum for Production of Ammonium Sulfate Fertilizer by Selim F. Estefan.
[8] Wikkiam C. Burnett, Michael S., Carter H." Behavior of Radionuclides During Ammonocarbon of PhốtphátGyps", Florida State Univ.
[9] Gopinath, N.D., 1968 Disposal or use of Gyps in production of ammonium sulfate. In: Phosphoric Acid, vol. l, Part II, (ed. A.V. Slack), Marcel Dekker, New York, 541-566. Chen, Jing-Yu, (1990). The trial production results from demonstration plant for production of ammonium sulfate from phosphogypsum. Phosphogypsum Conference Proceedings, Orlando,
[10] Florida Department of Environmental Protection (FDEP) management requirements for phosphogypsum stack systems (17-673.100-900, Florida Administrative Code)
[11] CPCB under the publication series HAZWAMS/17/2000-01 of Hazardous Waste Management i.e criteria for hazardous waste landfills.
[12] Guidelines for Development of Basic Infrastructure of Gypsum Handling and Storage on Railway Siding and Stack Yard of State Pollution Control Board, Orissa issued in August 2010
[13] Mineral Fertilizer Production and the Environment Part I. The Fertilizer Industry's Manufacturing Processes and Environmental Issues of United Nations Environmental Programme: ISBN: 92-807-1640-9.
tại KCN Đình Vũ – Hải Phịng”
Học viên: Lưu Ngọc Vĩnh Trang 82
[15] Fertilizer Annual Report 2009 & 2010 published by the ‘Fertilizer Association of India’
[16] ‘Study on The Production of Ammonium Sulfate Fertilizer from Phosphogypsum’ published by Khalid K. Abbas, Material Engineering Department, University of Technology, Baghdad in ‘Eng. & Tech. Journal, Vol.29, No.4, 2011. [17] Eco- balance features and significance of Hemihydrate Phosphogypsum Reprocessing into Gypsum Binding Materials’ S. Gaiducis et al. published in JOURNAL OF CIVIL ENGINEERING AND MANAGEMENT, 2009, 15(2): 205– 213
[18] IFA Technical Conference held at Marrakech, Morocco (during September-1 October 1998); Beijing, China (during 20-23 April 2004)
[19] Potential groundwater contamination by fluoride from two South African phosphogypsum by MP Motalane* and CA Strydom, Department of Chemistry,
University of Pretoria, Pretoria 0002, South Africa published in Water SA Vol. 30
No. 4 October 2004, 465-468
[20] Phosphate Fertilizer Plant :Pollution Prevention and Abatement Handbook WORLD BANK GROUP Effective July 1998
[21] ‘Phosphogypsum Stack Closure’ by Dean Kleinschmidt, Senior Civil
Engineer/Land Manager, Cargil Fertilizer, Inc. (USA.)
[22] Atomic Energy Regulatory Board Directive No. 01/09-Use of Phosphogypsum in Building & Construction Materials & in Agriculture (Ref. No. No. CH/AERB/IPSD/78/2009 dated March 19, 2009) downloaded from AERB website i.e http://www.aerb.gov.in/cgi-bin/constitution/directives.asp#dir3
[23] Indian Standard: Specification for By-product Gypsum IS:10170-1982 and Indian Standard: Specification for gypsum for agricultural use (first revision) IS:6046-1982
[24] The Fertilizer (Control) Order, 1985 published in the Gazette of India vide No. GSR 758 (E) dated the 25 th September 1985 and latest amended vide No. S.O. 2886 (E) dated 3 rd December 2010.
[25] Study on The Production of Ammonium Sulfate Fertlizer From Phosphogypsum Khalid K. Material Engineering Department, University of Technology / Baghdad 2011, 814-821
Học viên: Lưu Ngọc Vĩnh Trang 83
PHỤ LỤC
Phần mềm quy hoạch tuyến tính (phi tuyến) tối ưu hóa hai mức tối ưu
clear;clc;clf
p=0.05;%Muc co nghia 0.01 0.05 0.2
% -----Quy hoach tuyen tinh--------------------------------% % y=[91 90.3 93.7 93.4 90.85 89.3 93.95
94.5]%ke hoach tai nhan 2^k y0=[95.2
95.5
95.4]%Thi nghiem lap tai tam
%---------Quy hoach phi tuyen---------------------------------------------------- % y=[...%ke hoach tai nhan 2^k
% ...% Canh tay don 2k % ]';%Tam n0=1
% y0=[ ]';%Thi nghiem lap tai tam %----------------------------------------------------------- [hang cot]=size(y); switch hang case 2^3 qh31(p,y,y0); shg; case 2^3+2*3+1 qh32(p,y,y0); shg; case 2^4 qh41(p,y,y0);shg; case 2^4+2*4+1 qh42(p,y,y0);shg; end function qh31(p,yy,yy0) % So cac yeu to k=3;N=2^k;shs=k+2+(k-1)*k/2; % -------------------------------------%
tại KCN Đình Vũ – Hải Phịng”
Học viên: Lưu Ngọc Vĩnh Trang 84
% Toa do khong thu nguyen xx0=[1 1 1 1 1 1 1 1 ]';
xx1=[-1 +1 -1 +1 -1 +1 -1 +1 ]'; xx2=[-1 -1 +1 +1 -1 -1 +1 +1 ]'; xx3=[-1 -1 -1 -1 1 1 1 1]';
X=[xx0 xx1 xx2 xx3 xx1.*xx2 xx1.*xx3 xx2.*xx3 xx1.*xx2.*xx3]; %--he so pt hoi quy------
B=(X'*X)^(-1)*X'*yy;
% Danh gia tinh co nghia he so Student f2=length(yy0)-1; S11binh=sum((yy0-sum(yy0)/length(yy0)).^2)/f2 Sb=sqrt(S11binh/N) t=abs(B)/Sb hsvn=0; for i=1:shs if t(i)<student(p,f2) B(i)=0; hsvn=hsvn+1; end end Bkiemdinh=B l=shs-hsvn; f1=N-l f2 studentpf2=student(p,f2) fisherpf21=fisher(p,f2,f1)
% Kiem dinh tinh tuong hop Fisher
y=B(1)+B(2).*xx1+B(3).*xx2+B(4).*xx3+B(5).*xx1.*xx2+B(6).*xx1.*xx3+B(7).* xx2.*xx3+B(8).*xx1.*xx2.*xx3;
Sdubinh=(sum((yy-y).^2))/f1 F=Sdubinh/S11binh
%--Viet ra man hinh-------------- if F<fisher(p,f2,f1)
hold on;
fill([0 10 10 0 0],[0 0 10 10 0],'y','EdgeColor','none'); axis off; text(0.1,9,['PTHQ bac 1 tuong hop voi buc tranh thuc nghiem:']) text(0.5,8,['y = b0 ']);
text(1,7,'+ b1.x1 + b2.x2 + b3.x3 ');
text(1,6,'+ b12.x1.x2 + b13.x1.x3 + b23.x2.x3') text(1,5,'+ b123.x1.x2.x3 ')
Học viên: Lưu Ngọc Vĩnh Trang 85
hold on;
plot([0,10],[4 4],'r');
text(1,3.5, ['b0= ' num2str(B(1))])
text(1,2.5, ['b1= ' num2str(B(2))]); text(4,2.5, ['b2= ' num2str(B(3))]); text(7,2.5, ['b3= ' num2str(B(4))]) ;
text(1,1.5, ['b12= ' num2str(B(5))]); text(4,1.5, ['b13= ' num2str(B(6))]); text(7,1.5, ['b23= ' num2str(B(7))]) ;
text(1,0.5, ['b123= ' num2str(B(8))]) ; else
hold on;
fill([0 10 10 0 0],[0 0 10 10 0],'y','EdgeColor','none'); axis off;
text(0.5,9.5,['PTHQ bac 1 khong tuong hop, can tien hanh hoi quy phi tuyen']) end end function qh32(p,yy,yy0) %----Cac he so-------- k=3;n0=1 ;anpha=1.215; N=2^k+2*k+n0 ;shs=(k+1)*(k+2)/2+1;
%------Toa do khong thu nguyen---------------------------------- xx0=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ]'; xx1=[-1 +1 -1 +1 -1 +1 -1 +1 -anpha anpha 0 0 0 0 0 ]'; xx2=[-1 -1 +1 +1 -1 -1 +1 +1 0 0 -anpha anpha 0 0 0 ]'; xx3=[-1 -1 -1 -1 1 1 1 1 0 0 0 0 -anpha anpha 0]'; x1p=xx1.^2-sum(xx1.^2)/N; x2p=xx2.^2-sum(xx2.^2)/N; x3p=xx3.^2-sum(xx3.^2)/N; X=[xx0 xx1 xx2 xx3 xx1.*xx2 xx1.*xx3 xx2.*xx3 xx1.*xx2.*xx3 x1p x2p x3p]; %-----He so phuong trinh hoi quy---------------------------------------
B=inv(X'*X)*X'*yy;
%---Kiem dinh tinh co ngia Student f2=length(yy0)-1; S11binh=sum((yy0-sum(yy0)/length(yy0)).^2)/f2 Sbop=sqrt(S11binh/(sum(xx0.^2))); Sbi=sqrt(S11binh/(sum(xx1.^2))); Sbju=sqrt(S11binh/(sum((xx2.*xx3).^2))); Sbjj=sqrt(S11binh/(sum(x1p.^2))); Sbiju=sqrt(S11binh/sum((xx1.*xx2.*xx3).*2)); % Sb=sqrt(Sbop^2+Sbjj^2*(sum(xx1.^2)/N))^2 t(1)=abs(B(1))/Sbop; t(8)=abs(B(8)/Sbiju);
tại KCN Đình Vũ – Hải Phòng”
Học viên: Lưu Ngọc Vĩnh Trang 86
for i=2:4 t(i)=abs(B(i))/Sbi; t(i+3)=abs(B(i+3))/Sbju; t(i+7)=abs(B(i+7))/Sbjj; end f2=length(yy0)-1; hsvn=0; for i=1:shs if t(i)<student(p,f2) B(i)=0; hsvn=hsvn+1; end end B(1)=B(1)-B(9)*(sum(xx1.^2)/N)-B(10)*(sum(xx2.^2)/N)-B(11)*(sum(xx3.^2)/N); B l=shs-hsvn; f1=N-l f2 studentpf2= student(p,f2) fisherpf2f1=fisher(p,f2,f1)
% Kiem dinh tinh tuong hop Fisher
y=B(1)+B(2).*xx1+B(3).*xx2+B(4).*xx3 + B(5).*xx1.*xx2 + B(6).*xx1.*xx3 + B(7).*xx2.*xx3... + B(8).*xx1.*xx2.*xx3... + B(9).*xx1.^2+ B(10)*xx2.^2+ B(11).*xx3.^2; Sdubinh=(sum((yy-y).^2))/(N-f1) F=Sdubinh/S11binh
%--Viet ra man hinh-------------- if F<fisher(p,f2,f1)
hold on;
fill([0 10 10 0 0],[0 0 10 10 0],'y','EdgeColor','none'); axis off; text(0.5,9.7,['PTHQ bac 2 tuong hop voi buc tranh thuc nghiem:']) text(0.8,9,['y = b0 ']) text(1,8.2,[' + b1.x1 + b2.x2 + b3.x3']) text(1,7.4,[' + b12.x1.x2 + b13.x1.x3 + b23.x2.x3 ']) text(1,6.6,[' + b123.x1.x2.x3']) text(1,5.8, [' + b11.x1^2 + b22.x2^2+b33.x3^2']) hold on; plot([0 10],[5.5 5.5],'r'); shg; text(1,5, ['b0= ' num2str(B(1))]);
Học viên: Lưu Ngọc Vĩnh Trang 87
text(1,4, ['b1= ' num2str(B(2))]); text(4,4, ['b2= ' num2str(B(3))]);text(7,4, ['b3= ' num2str(B(4))]);
text(1,3, ['b12= ' num2str(B(5))]);text(4,3, ['b13= ' num2str(B(6))]); text(7,3, ['b23= ' num2str(B(7))]);
text(1,2, ['b123= ' num2str(B(8))]);
text(1,1, ['b11= ' num2str(B(9))]); text(4,1, ['b22= ' num2str(B(10))]); text(7,1, ['b33= ' num2str(B(11))]);
else
hold on;
fill([0 10 10 0 0],[0 0 10 10 0],'y','EdgeColor','none'); axis off;
text(0.5,9.5,['PTHQ bac 2 khong tuong hop voi buc tranh thuc nghiem']) end end function qh41(p,yy,yy0) % So cac yeu to k=4;N=2^k;shs=k+1+(k-1)*k/2; % -------------------------------------% % Toa do khong thu nguyen
xx0=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]'; xx1=[-1 +1 -1 +1 -1 +1 -1 +1 -1 +1 -1 +1 -1 +1 -1 +1]'; xx2=[-1 -1 +1 +1 -1 -1 +1 +1 -1 -1 +1 +1 -1 -1 +1 +1 ]'; xx3=[-1 -1 -1 -1 1 1 1 1 -1 -1 -1 -1 1 1 1 1]'; xx4=[-1 -1 -1 -1 -1 -1 -1 -1 1 1 1 1 1 1 1 1]'; X=[xx0 xx1 xx2 xx3 xx4 xx1.*xx2 xx1.*xx3 xx1.*xx4 xx2.*xx3 xx2.*xx4 xx3.*xx4];
%--he so pt hoi quy------ B=(X'*X)^(-1)*X'*yy;
% Danh gia tinh co nghia he so Student f2=length(yy0)-1; S11binh=sum((yy0-sum(yy0)/length(yy0)).^2)/f2; Sb=sqrt(S11binh/(2^k)); t=abs(B)./Sb; hsvn=0; for i=1:shs if t(i)<student(p,f2) B(i)=0; hsvn=hsvn+1; end end Bkiemdinh=B
tại KCN Đình Vũ – Hải Phịng”
Học viên: Lưu Ngọc Vĩnh Trang 88
l=shs-hsvn f1=N-l f2 studentpf2=student(p,f2) fisherpf21=fisher(p,f2,f1) y=B(1)... + B(2).*xx1 + B(3).*xx2 + B(4).*xx3 +B(5).*xx4... + B(6).*xx1.*xx2 + B(7).*xx1.*xx3 + B(8).*xx1.*xx4 + B(9).*xx2.*xx3 + B(10).*xx2.*xx4 + B(11).*xx3.*xx4; Sdubinh=(sum((yy-y).^2))/f1; F=Sdubinh/S11binh
%--Viet ra man hinh-------------- if F<fisher(p,f2,f1)
fill([0 10 10 0 0],[0 0 10 10 0],'y','EdgeColor','none'); axis off; text(0.1,9,['PTHQ bac 1 tuong hop voi buc tranh thuc nghiem:']) text(0.5,8,['y = b0 ']); text(1,7,'+ b1.x1 + b2.x2 + b3.x3 + b4.x4'); text(1,6,'+ b12.x1.x2 + b13.x1.x3 + b14.x1.x4 '); text(1,5,'+ b23.x2.x3 + b24.x2.x4 + b34.x3.x4'); hold on; plot([0,10],[4 4],'r'); text(1,3.5, ['b0= ' num2str(B(1))])
text(1,2.5, ['b1= ' num2str(B(2))]); text(3,2.5, ['b2= ' num2str(B(3))]); text(5,2.5, ['b3= ' num2str(B(4))]) ; text(7,2.5, ['b4= ' num2str(B(5))]) ;
text(1,1.5, ['b12= ' num2str(B(6))]); text(3,1.5, ['b13= ' num2str(B(7))]); text(5,1.5, ['b14= ' num2str(B(8))]);
text(1,0.5, ['b23= ' num2str(B(9))]);text(3,0.5, ['b24= ' num2str(B(10))]);text(5,0.5, ['b34= ' num2str(B(11))]); else
fill([0 10 10 0 0],[0 0 10 10 0],'y','EdgeColor','none'); axis off;
text(0.5,9.5,['PTHQ bac 1 khong tuong hop, can tien hanh hoi quy phi tuyen']) end end function f=fc(x) % ------------------- B=[72.0125 -3.3625 2.2875 -3.0125 -0.1375
Học viên: Lưu Ngọc Vĩnh Trang 89 1.7625 -0.1875 -0.0625]; % -------------------------- [hang cot]=size(B); % [xx0 xx1 xx2 xx3 xx4 xx1.*xx2 xx1.*xx3 xx1.*xx4 xx2.*xx3 xx2.*xx4 xx3.*xx4] switch hang case 8 f=-sum([1 x(1) x(2) x(3) x(1)*x(2) x(1)*x(3) x(2)*x(3) x(1)*x(2)*x(3) ].*B'); case 11 f=-sum([1 x(1) x(2) x(3) x(1)*x(2) x(1)*x(3) x(2)*x(3) x(1)*x(2)*x(3) x(1)^2 x(2)^2 x(3)^2].*B'); end clear;clc; %___________________ %Toc do khuay Z1=[30 60];
%Thoi gian phan ung Z2=[1 4];
% Nhiet do phan ung Z3=[30 45]; %______________________ x0=[0; 0; 0]; lb=[-1 -1 -1]; ub=[1 1 1]; %________________________________ %------------- [x, fmin]=fmincon('fc',x0,[],[],[],[],lb,ub); Q1=x(1)*(Z1(2)-Z1(1))/2+(Z1(1)+Z1(2))/2; Q2=x(2)*(Z2(2)-Z2(1))/2+(Z2(1)+Z2(2))/2; T=x(3)*(Z3(2)-Z3(1))/2+(Z3(1)+Z3(2))/2; %_______________________________________
fill([0 10 10 0 0],[0 0 10 10 0],'y','EdgeColor','none'); axis off; shg;hold on;shg; text(3,9.3,'CHUONG TRINH TOI UU HOA');
plot([0 10],[8.8 8.8],'r')
text(3,0.5,'CNDT LUU NGOC VINH') plot([0 10],[0.9 0.9],'r')
tại KCN Đình Vũ – Hải Phịng”