numerical-methods

    0

    1답변

    R에서 광범위한 계산을 시도했습니다. 18 시간이 지났지 만 RStudio는 계속 작동하는 것 같습니다. 스크립트를 다른 방식으로 작성하여 더 빠르게 만들 수 있는지 확신 할 수 없습니다. #defining the discretization of cells dt<-1 t<-50000 dz<-0.0075 z<-350*dz #velocity & d

    3

    1답변

    구체적으로 다음 두 종류의 코드는 동일한 S 및 V 아이디어를 얻을 수 있습니다. 그러나 두 번째 속도는 일반적으로 Matlab의 첫 번째 속도보다 빠릅니다. 누군가 이유를 말해 줄 수 있습니까? 게다가 어떤 방법이 수치 적으로 더 안정적입니까? 감사합니다. . [~,S,V] = svd(B,'econ'); [Qc,Rc] = qr(B',0); [U,S

    1

    1답변

    내 문제는 구문 문제라고 생각하기 때문에 여기에 게시 : 나는 2 차 ODE를 분석하려고합니다. y ''+ 2y '+ 2y = e^(- x) * sin (x) midpoint 메서드에 대한 MATLAB 코드를 사용하여. 먼저 ODE를 1 차 방정식 시스템으로 변환 한 후 아래에 적용하려고 시도했지만 이산 값 [m]이 증가함에 따라 출력은 .2718에서 멈추

    1

    1답변

    y=x^2+xe^(x)에 대한 MATLAB에서 라그랑주 보간법을 계산하려고했습니다. 그래서 x0=4.7에 대한 결과가 있었다 clc clear close all x0=4.7; n=10; x=linspace(0,5,n); y=x.^2+x.*exp(x); syms t L=sym(ones(1,n)); P_x=sym(0); for i=1:n

    0

    1답변

    에 나는 매트랩 라그랑주 보간 코드가 있습니다 function[r] = interpolation(x, y) r = 0; for i = 1:size(x,2) w = 1; for j = 1:size(x,2) if j ~= i w = conv(w, [1, -x(j)])/(x(i) - x(j));

    2

    1답변

    아주 적은 기능 평가를 사용하는 root-finding 알고리즘을 찾고 있습니다 (목표는 최소 임). 루트 - 발견 된 문제는 다음과 같은 특징이 있습니다 : f(x) = 0, R -> R 기능 (f(.)) 평가는 매우 비용이 많이 드는 *; 경계 간격 ([a,b])은 시작에 사용할 수 있습니다 (비교적 좋은 근사값, 야생 추정치는 아닙니다). f(.)은

    0

    1답변

    시스템 인덱스가 범위를 벗어났습니다. 예외 오류가 발생하여 어디서 실수를했는지 찾을 수 없었습니다. 아마도 인덱스에 오류가있었습니다. 저를 도와주세요? 당신의 친절에 감사드립니다. 내가 처음 여기에 오류가 발생합니다 : x1 = x1 - (hess_universe[1, 1] * grad[1] + hess_universe[1, 2] * grad[2] + he

    0

    2답변

    : log(a + b) = log(a * (1 + b/a)) = log a + log(1 + b/a) 나는이 방법이 더 효율적입니다 이유에 대한 혼란 스러워요 .. 사람이에 대한 아이디어가 있습니까 : To Compute log(a+b)이 때때로 log_sum는 다음과 같이 구현된다?

    0

    1답변

    최종 상태가 주어진 시스템의 초기 상태와 시스템의 다양한 측면의 변화율을 설명하는 비선형 미분 방정식 시스템을 풀 수 있습니까? 예를 들어 : V_x' = n * cos(b) V_y' = (n * sin(b)) - (g_0 * (r_0/r)^2) b' = ((g_0 * (r_0/r)^2) * cos(b))/||V|| b_final = 0 V_y_f

    -2

    1답변

    익명 함수를 저장하는 변수가 있습니다. temp_f = @(x)((1/h^3)*(((xi+2*i)-x)^3-4*((xi+i)-x)^3)) 나는 다음 단계에서, 그래서 내가 그 접근 할 수있는 방법 temp2_f = @(x) (temp_f(x) * varA) 를 얻을 수있는 'VARA'에 의해 곱 싶습니다 ? 당신이 varA 다음 새로운 익명 함수의