MATLAB và SciPy đưa ra các kết quả khác nhau cho hàm 'buttord'

In lại Tác giả: Vũ trụ không gian Thời gian cập nhật: 2023-11-03 20:25:02 29 4

Wp = 2*pi*10
Ws = 2 * pi * 100
Rp = 3
R = 80
[N, Wn] = mông(Wp, Ws, Rp, Rs, 's')

từ số pi nhập numpy
từ tín hiệu nhập scipy
Wp = 2*pi*10
Ws = 2 * pi * 100
Rp = 3
R = 80
(N, Wn) = signal.buttord(Wp, Ws, Rp, Rs, analog=True)

nhập numpy dưới dạng np
từ nhập khẩu scipy.signal buttord, bơ, tần số
nhập matplotlib.pyplot dưới dạng plt


# Kết quả thiết kế cho:
Wp = 2*np.pi*10
Ws = 2*np.pi*100
Rp = 1,5 # thay vì 3
R = 80

n_old = 5
wn_old = 99,581776302787929

n_new, wn_new = buttord(Wp, Ws, Rp, Rs, analog=True)

b_old, a_old = bơ(n_old, wn_old, analog=True)
w_old, h_old = tần số(b_old, a_old)

b_new, a_new = bơ(n_new, wn_new, analog=True)
w_new, h_new = freqs(b_new, a_new)


db_old = 20*np.log10(np.abs(h_old))
db_new = 20*np.log10(np.abs(h_new))

plt.semilogx(w_old, db_old, 'b--', label='old')
plt.axvline(wn_old, color='b', alpha=0,25)
plt.semilogx(w_new, db_new, 'g', label='new')
plt.axvline(wn_new, color='g', alpha=0,25)

plt.axhline(-3, color='k', ls=':', alpha=0.5, label='-3 dB')

plt.xlim(40, 1000)
plt.ylim(-100, 5)

xbounds = plt.xlim()
ybounds = plt.ylim()
trực tràng = plt.Rectangle((Wp, ybounds[0]), Ws - Wp, ybounds[1] - ybounds[0],
                     facecolor="#000000", edgecolor='none', alpha=0.1, Hatch='//')
plt.gca().add_patch(trực tràng)
trực tràng = plt.Rectangle((xbounds[0], -Rp), Wp - xbounds[0], 2*Rp,
                     facecolor="#FF0000", edgecolor='none', alpha=0,25)
plt.gca().add_patch(trực tràng)
trực tràng = plt.Rectangle((Ws, ybounds[0]), xbounds[1] - Ws, -Rs - ybounds[0],
                     facecolor="#FF0000", edgecolor='none', alpha=0,25)
plt.gca().add_patch(trực tràng)

plt.annotate("Đạt", (0,5*(xbounds[0] + Wp), Rp+0,5), ha='center')
plt.annotate("Dừng", (0,5*(Ws + xbounds[1]), -Rs+0,5), ha='center')
plt.annotate("Không quan tâm", (0,1*(8*Wp + 2*Ws), -Rs+10), ha='center')

plt.legend(loc='tốt nhất')
plt.xlabel('Tần số [rad/s]')
plt.ylabel('Tăng [dB]')
plt.show()

29 4 0

• matlab - Các tệp tập lệnh mẫu trong Matlab có khác với các tệp tập lệnh chuẩn trong Matlab không?

  Trong Matlab, bạn có tùy chọn tạo các tệp tập lệnh "ví dụ" mới cũng như các tập lệnh, hàm, lớp, v.v. Khi bạn tạo chúng, chúng sẽ có biểu tượng tập lệnh. Chúng được xử lý khác với các tệp tập lệnh tiêu chuẩn khác như thế nào? Có tài liệu nào về mục đích sử dụng của các loại tập lệnh mẫu này không?

• MATLAB - Có một số quy trình Matlab nhất định được sử dụng trong các tập lệnh Matlab không?

  Tôi đang chạy một tệp m lớn mà tôi không tự viết và nó dựa vào một số chức năng phụ. Tôi muốn biết liệu một hàm cụ thể (trong trường hợp của tôi là hàm eig.m (tính giá trị riêng)) có được sử dụng ở bất kỳ đâu trong tất cả các hàm lồng nhau hay không. Có cách nào nhanh chóng để làm điều này? liên quan đến

• matlab - Cách xem mã Matlab đằng sau hàm Matlab

  Có một hàm trong Matlab tên là copulafit. Làm cách nào tôi có thể xem mã đằng sau hàm này? Nhiều hàm numpy và scipy của Python dễ dàng là mã nguồn mở trên Github, nhưng vì lý do nào đó tôi

• matlab - Bài toán kế thừa OOP Matlab của bộ mã hóa Matlab (Khối hàm Simulink MATLAB)

  Tôi đã định nghĩa một hàm xử lý MeasureHandler < của lớp cơ sở trừu tượng để xác định giao diện cho tất cả các lớp được kế thừa. Hai lớp con của lớp này là a < MeasureHandler và b < me

• matlab - Làm cách nào để lặp qua các cột của ma trận trong matlab và thêm chúng riêng lẻ vào các cột cụ thể của ma trận tổng trong matlab?

  Giả sử có một ma trận A = 1 3 2 4 4 2 5 8 6 1 4 9 Ví dụ mình có một Vector chứa "class" của mỗi cột của ma trận

• matlab - Vô hiệu hóa việc đánh cắp tiêu điểm Matlab mà không sửa đổi tập lệnh Matlab

  Tôi có một tập lệnh Matlab đang chạy ở chế độ nền. Khi quá trình tính toán tiếp tục, nó sẽ tiếp tục xuất hiện trên cửa sổ thanh tiến trình. Điều này rất khó chịu. Vấn đề là tôi không tự viết kịch bản Matlab, nó là một đoạn mã dài và phức tạp và tôi không muốn làm nó rối tung lên. Vậy làm thế nào để làm điều đó mà không cần sửa đổi Matla

• matlab - Phát hiện các quy trình matlab từ bên trong matlab

  Có cách nào để phát hiện từ chương trình matlab có bao nhiêu tiến trình matlab đang chạy trên máy tính không? Tôi muốn có chính xác n tiến trình matlab đang chạy. Nếu tôi có quá ít tôi muốn tạo ra chúng, nếu có quá nhiều tôi muốn giết một số. Bạn nên

• matlab - Ứng dụng Matlab được triển khai sử dụng nhiều bộ nhớ hơn tập lệnh Matlab

  Tôi đang thử nghiệm một ứng dụng độc lập mà chúng tôi đang phát triển trong Matlab thì nhận thấy rằng mức sử dụng bộ nhớ của nó (theo Trình quản lý tác vụ Windows) đã đạt đỉnh nhiều lần trên 16gb. Tôi quyết định sử dụng hồ sơ trên tập lệnh đằng sau phiên bản đã biên dịch

• matlab - Chữ cái tiếng Nga trong cửa sổ lệnh Matlab

  Tôi đang phải đối mặt với một vấn đề khá khó khăn. Trong Matlab, lệnh S = char(1044) đặt chữ cái tiếng Nga д vào biến S. Nhưng disp(S) trả về các ký hiệu khoảng trắng, mặc dù nội dung thực sự đúng: >> S = char

• MATLAB có thực hiện các phép toán mảng cơ bản trong không gian không đổi không?

  Tôi đã gặp phải lỗi hết bộ nhớ trong dòng mã MATLAB này: result = (A(1:xmax,1:ymax,1:zmax) .* B(2:xmax+1,2:ymax+1,2:zmax+ 1) +

• Có kiểu dữ liệu giống như tập hợp trong MATLAB không?

  Tôi đang tìm cách so sánh dữ liệu thứ tự hữu hạn với thứ tự không xác định trong MATLAB. Về cơ bản, thứ tôi muốn là một mảng nhưng không áp đặt thứ tự cho các phần tử được chứa. Nếu tôi có đối tượng a = [xyz] và b = [xzy];

• matlab - Thay thế các giá trị trùng lặp dựa trên kích thước chuỗi - Matlab

  Tôi có ma trận 2D bao gồm 1 và 0. mat = [0 0 0 0 1 1 1 0 0 1 1 1 1 1 0 0 1 0 0 0 1 0 1 1 0 0 1];

• Đường hồi quy cho cụm điểm trong Matlab

  Tôi có thể vẽ đường hồi quy trong Matlab với tập hợp các điểm x, y. Tuy nhiên, nếu tôi có một tập hợp điểm (như hiển thị bên dưới), giả sử tôi có bốn tập hợp điểm và tôi muốn vẽ bốn đường hồi quy cho chúng...tôi phải làm thế nào? Tất cả các điểm được lưu trong x, y. Không có cách nào để tách chúng ra và đặt chúng

• c# — Khoảng trống trong sơ đồ bề mặt MATLAB

  Tôi đang cố vẽ một hình nón trong MATLAB bằng đoạn mã sau. Tuy nhiên, khi MATLAB tạo đồ thị, có một khoảng trống trên bề mặt, như minh họa trong hình bên dưới. Có ai có thể gợi ý cách tắt nó đi không? clearvars; đóng tất cả;

• Ánh xạ trong MATLAB là bao nhiêu byte?

  Tôi có một bản đồ tên là res_Map chứa một tập hợp các mảng có kích thước khác nhau. Tôi muốn tìm tổng bộ nhớ được sử dụng để lưu trữ res_Map. Như bạn có thể thấy bên dưới, có vẻ như res_Map hầu như không chiếm bộ nhớ và các phần tử riêng lẻ trong res_Map

• matlab - Cú pháp vectơ sang ma trận trong MATLAB

  Có cách nào để kết hợp 2 vectơ trong MATLAB như thế này không: mat = zeros(length(C),length(S)); for j=1:length(C) mat(j,:)=C(j)*S ;

• matlab - Phân rã cực trong MATLAB

  Đã kết thúc. Câu hỏi này không tuân thủ nguyên tắc Stack Overflow. Hiện tại nó không chấp nhận câu trả lời. Chúng tôi không cho phép các câu hỏi tìm kiếm đề xuất về sách, công cụ, thư viện phần mềm, v.v. Bạn có thể chỉnh sửa câu hỏi của mình để trả lời bằng sự kiện và trích dẫn. Đã đóng cửa 5 năm trước

• MatLab: Choline phải xác định dương

  Tôi đang cố gắng khớp một t copula trong MatLab với dữ liệu của mình và các hàm của tôi là: u = ksdense(range_1, range_1,'function','cdf');

• Giảm hình ảnh hình thái đa tỷ lệ trong Matlab

  Xin chào mọi người, tôi hiện đang cố gắng sử dụng toán tử SMMT Dorini từ bài báo "Giảm hình ảnh đa tầng". Vì trang này không thể truy cập được nếu không đăng ký nên tôi sẽ đăng các chi tiết có liên quan tại đây: Xin lưu ý rằng tôi đang đăng. bài viết liên quan Một số nội dung được xuất bản dưới dạng hình ảnh.

• MATLAB: Sử dụng các phiên bản khác nhau của hàm có cùng tên

  Tôi đang viết mã bằng MATLAB và cần sử dụng hàm có tên modwt tồn tại trong hai hộp công cụ tôi sử dụng cùng lúc (Wavelet và WMTSA). Vấn đề là một phiên bản chỉ trả về một đầu ra, trong khi phiên bản kia trả về ba đầu ra. tôi nên