Tôi cần tạo một danh sách các biến ngẫu nhiên theo các ràng buộc có thể được biểu thị dưới dạng trong đó là ma trận nếu có mục. Trong tất cả các trường hợp tôi đang xử lý, , ví dụ sẽ vào khoảng 14.000 và sẽ là 50. Tôi không chắc tôi sẽ sử dụng phương pháp nào để lấy mẫu ngẫu nhiên, bình thường hay thống nhất, không phải là rõ ràng cái nào là tốt nhất cho vấn đề tôi đang cố gắng giải quyết, nhưng tôi cần tất cả các biến được lấy mẫu từ các bản phân phối có cùng giá trị trung bình và phạm vi / phương sai.
Những gì tôi đã làm để giải quyết vấn đề này là giảm thành dạng phản hồi hàng, đặt tất cả các biến tương ứng với các cột ở bên phải của trục cuối cùng thành các giá trị ngẫu nhiên, sau đó giải phương trình ma trận vuông còn lại.
Tuy nhiên, có một vấn đề, để giải quyết đẳng thức ma trận vuông, tôi trừ các giá trị đã đặt từ phía bên phải. Thật không may, phương sai cũng thêm vào, vì vậy 50 giá trị cuối cùng của tôi có xu hướng thay đổi lớn, điều không may là không thể chấp nhận được trong vấn đề này.
Có cách nào tốt hơn để làm điều này? Tôi không thể nghĩ ra cách khắc phục phương pháp hiện tại tôi đang sử dụng. Tôi sử dụng R.